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Syllabus 

Fluid Mechanics for Master Students 
https://fluidmech.ninja/ 

Summer semester 2018 
Olivier Cleynen 

Welcome to the Fluid Mechanics course of the Chemical and Energy Engineering Tprograxnl 


Objectives 


This one-semester course is designed for students working towards a Master in engineer¬ 
ing, who have had little or no previous experience with fluid dynamics. My objectives 
are simultaneously: 

• to provide you with a solid understanding of what can, an cannot, be done with 
fluid mechanics in engineering; 

• to enable you to solve selected engineering fluid dynamics problems with confi¬ 
dence; 

• to request the minimum feasible amount of your time and energy for the above. 

If all goes well, at the end of the semester, you should be well-prepared to begin a course 
in Computational Fluid Dynamics, where the knowledge and skills you acquire here can 
be used to solve applied problems in great detail. 


Course homepage 


The course homepage is located at a hard-to-unremember location: 

https://fluidmech.ninja/ 

From that page, you can download the latest version of the present document, as well as 
lecture slides and archived examinations. 

The course script is the only document you need to work through and succeed in this 
course. I highly recommend that you print it out (at 200 pages, a copy at the University 
print shop costs less than 9 €). If you can’t afford this, print at least the exercise sheets 
and the first page of each chapter. 


Lectures 


We will have lectures on Thursdays from 11:15 to 12:45 in room G05-H4. In class, we will 
cover a lot of content, at a rather brisk pace. The most efficient way to proceed is to read 
the notes before the lecture. Try doing that: even a brief skim over the relevant chapter 
helps. We will also experiment with new tools: if you have one at hand, try bringing an 
Internet-connected device with you. 
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Exercise sessions 


Exercise sessions take place on Fridays in room G22A-209. The schedule is slightly more 
complicated: 

• For first-semester students, on odd weeks, from 13:00 to 15:00, and on even weeks, 
from 15:00 to 17:00; 

• For second (or later) semester students, just the opposite: on even weeks, from 
13:00 to 15:00, and on odd weeks, from 15:00 to 17:00. 

I am available during those sessions to help you work through the exercises, and will 
strive to answer all of your questions. The objective should be for you to work as fast as 
possible in order to cover each sheet entirely; attempting the first exercises on the prior 
day helps greatly in this respect. Although I occasionally use the board and projector, 
these sessions are hands-on and you should expect to work autonomously, alone or in 
groups. Bring your own exercise sheets, a calculator, a sugary drink, and start right 
away! 


Assessment 


The end-of-semester examination will take place on July 12*, from 15:00 to 17:00 in 
room G26-HS1. This examination is the only assessment for the course. It lasts exactly 
2 hours and is closed-book (a formula sheet is provided, containing the formulas and data 
contained in the preamble of every exercise sheet). A briefing regarding the examination 
will take place in class immediately after we cover Chapter 5, on May 24*. 

The examination consists exclusively of lightly-modified exercise sheet problems. As 
reported by students in the previous years, “the exam is very hard if you haven’t done the 
exercises, and very easy if you have done them”. It rewards your ability to methodically, 
resolutely and reliably solve the problems studied in class. To succeed, make sure you 
can solve all of them without outside help! A handful of archived examinations, together 
with their full solutions, is provided on the course homepage to help you practice. 


Books 


The present document is, I hope, enough to allow you to work and succeed through the 
entire course. If you need a second source, I recommend the books of Crowe et al. [10] 
(particularly since a dozen copies are available in our library). White [13], Qengel 
et al. [16], or Munson et al. [18] (a few personal copies of which can be borrowed from 
my office). These are solid, beautifully-illustrated, well-established textbooks, with plenty 
of worked-out exercises. You will see that this course matches closely their content and 
style. 
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Time plan 


Our time plan, which is affected by the holidays in Saxony-Anhalt, should be as follows: 


Week 14: Ch. 0 (Important concepts) 

Week 15: Ch. 1 (Effects of pressure), No exercise session 
Week 16: Ch. 2 (Effects of shear) 

Week 17: Ch. 3 (Analysis of existing flows) part 1/2 
Week 18: Ch. 3 (Analysis of existing flows) part 2/2 
Week 19: No lecture 
Week 20: Ch. 4 (Predicting fluid flow) 

Week 21: Ch. 5 (Elow in pipes). Exam briefing 

Week 22: No lecture 

Week 23: Ch. 6 (Scale effects) 

Week 24: Ch. 7 (Elow near walls) 

Week 25: Ch. 8 (Large- and small-scale flows). Feedback 
Week 26: Ch. 9 (Introduction to compressible flow) 

Week 27: revision & practice 

Week 28: Examination: July 12* from 15:00 to 17:00 in room G26-HS1 


Copyright and remixing 


This document still contains several figures from cited, fully-copyrighted works. The 
source files for all of the remaining figures (which stem from many different authors, not 
all associated with this course) can be accessed and re-used by clicking the links in the 
captions. The text of this document is licensed under a Creative Commons CC-by-sa 
(attribution, share-alike) license. The Latex sources can be downloaded from the git 
repository accessed from the course homepage. I am grateful for receiving all types of 
feedback, including reports of inaccuracies and errors. 


Contact 


I can always be reached by email (olivier.cleynen@ovgu.de). It is far more efficient to 
ask questions during exercise sessions. If you cannot make it there or for urgent matters, 
contact me directly by mail (please be precise and succinct) or afternoons at my office 
G14-111. 

It’s a pleasure to join you for this course this semester! Fluid mechanics is one of the 
most exciting disciplines out there. Now, let us begin! 


Olivier Cleynen 
April 2018 
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Fluid Mechanics 
Chapter 0 - Important concepts 
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These lecture notes are based on textbooks by White [13], Qengel & al.[16], and Munson & al.[18]. 


0.1 Concept of a fluid 


We call fluid a type of matter which is continuously deformable, and which 
spontaneously tends to adapt its shape to its container by occupying all of 
the space made available to it. 


0.2 Fluid mechanics 


0.2.1 Solution of a flow 

Fluid mechanics is the study of fluids subjected to given constraints. The 
most common type of problem in this discipline is the search for a complete 
description of the fluid flow around or through a solid object. This problem is 
solved when the entire set of velocities of fluid particles has been described. 
This set of velocities, which is often a function of time, can be described either 
as a set of discrete values (“pixelized” data) or as a mathematical function; it 
is called the flow solution. 
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If the solution is known, the shear and pressure efforts generated on the 
surface of the object can be calculated. Other quantities, such as the force 
and moments applying on the object, or the fluid’s energy gains or losses, 
can also be calculated. 

0.2.2 Modeling of fluids 

Like all matter, fluids are made of discrete, solid molecules. However, within 
the scope of fluid mechanics we work at the macroscopic scale, in which 
matter can be treated like a continuum within which all physical properties 
of interest can be continuously differentiable. 

There are about 2 • 10^^ molecules in the air within an “empty” 1-liter bottle 
at ambient temperature and pressure. Even when the air within the bottle is 
completely still, these molecules are constantly colliding with each other and 
with the bottle walls; on average, their speed is equal to the speed of sound: 
approximately 1 000 km/h. 

Despite the complexity of individual molecule movements, even the most 
turbulent flows can be accurately described and solved by considering the 
velocities of groups of several millions of molecules collectively, which we 
name fluid particles. By doing so, we never find out the velocity of individual 
molecules: instead, those are averaged in space and time and result in much 
simpler and smoother trajectories, which are those we can observe with 
macroscopic instruments such as video cameras and pressure probes. 

Our selection of an appropriate fluid particle size (in effect defining the lower 
boundary of the macroscopic scale), is illustrated in fig. 0.1. We choose to 
reduce our volume of study to the smallest possible size before the effect of 
individual molecules becomes meaningful. 

Adopting this point of view, which is named the continuum abstraction, is not 
a trivial decision, because the physical laws which determine the behavior 
of molecules are very different from those which determine the behavior of 
elements of fluid. For example, in fluid mechanics we never consider any 
inter-element attraction or repulsion forces; while new forces “appear” due 
to pressure or shear effects that do not exist at a molecular level. 

A direct benefit of the continuum abstraction is that the mathematical com¬ 
plexity of our problems is greatly simplified. Finding the solution for the 



Figure 0.1 - Measurement of the average value of a property (here, velocity V; but it 
could be pressure, or temperature) inside a given volume. As the volume shrinks 
towards zero, the fluid can no longer be treated as a continuum; and property 
measurements will oscillate wildly. 

Figure CC-by-sa o.c. 
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bottle of “still air” mentioned above, for example, requires only a single 
equation (section 1.3.3 p.33) instead of a system of 2 • 10^^ equations with 
2 • 10^^ unknowns (all leading up to 14 verage,x,t;,z,t = 0!)- 

Another consequence is that we cannot treat a fluid as if it were a mere set of 
marbles with no interaction which would hit objects as they move by. Instead 
we must think of a fluid -even a low-density fluid such as atmospheric air- 
as an infinitely-flexible medium able to fill in almost instantly all of the space 
made available to it. 

0.2.3 Theory, numerics, and experiment 

Today, fluid dynamicists are typically specializing in any one of three sub¬ 
disciplines: 

Analytical fluid mechanics which is the main focus of these lectures and 
which consists in predicting fluid flows mathematically. As we shall 
see, it is only able to provide (exact) solutions for very simple flows. 
In fluid mechanics, theory nevertheless allows us to understand the 
mechanisms of complex fluid phenomena, describe scale effects, and 
predict forces associated with given fluid flows; 

Numerical fluid mechanics also called Computational Fluid Dynamics or 
CFD, which consists in solving problems using very large numbers of 
discrete values. Initiated as a research topic in the 1970s, cfd is now 
omnipresent in the industry; it allows for excellent visualization and 
parametric studies of very complex flows. Nevertheless, computational 
solutions obtained within practical time frames are inherently approxi¬ 
mate: they need to be challenged using analysis, and calibrated using 
experimental measurements; 

Experimental fluid mechanics which consists in reproducing phenom¬ 
ena of interest within laboratory conditions and observing them using 
experimental techniques. A very mature branch (it first provided useful 
results at the end of the 19* century), it is unfortunately associated 
with high human, equipment and financial costs. Experimental mea¬ 
surements are indispensable for the validation of computational simu¬ 
lations; meanwhile, the design of meaningful experiments necessitates 
a good understanding of scale effects. 

Our study of analytical fluid mechanics should therefore be a useful tool to 
approach the other two sub-disciplines of fluid mechanics. 


0.3 Important concepts in mechanics 


Mechanics in general deals with the study of forces and motion of bodies. A 
few concepts relevant for us are recalled here. 

0.3.1 Position, velocity, acceleration 

The description of the movement of bodies (without reference to the causes 
and effects of that movement) is called kinematics. 



The position in space of an object can be fully expressed using three com¬ 
ponents (one for each dimension) This is usually done with a vector (here 

written x). If the object moves, then this vector varies with time. 

—> 

The velocity V of the object is the rate of change in time of its position: 




dx 

dt 


( 0 / 1 ) 


The magnitude V of velocity V is called speed and measured in m s“^ Con¬ 
trary to speed, in order to express velocity completely, three distinct values 
(also each in m s“^) must be expressed. Generally, this is done by using one 
of the following notations: 





I Ux \ 

Uy 

\ j 

(Ui) = {Ui,U2,U3) 
( Ur \ 

ue 

\ I 


The acceleration a of the object is the rate of change in time of its velocity: 


a = 


dV 

dt 


( 0 / 2 ) 


Acceleration is especially important in mechanics because it can be deduced 
from Newton’s second law (see eq. 0/23 p.l9 below) if the forces applying 
on the object are known. Acceleration can then be integrated with respect 
to time to obtain velocity, which can be integrated with respect to time to 
obtain position. 

Like velocity, acceleration has three components (each measured in m s“^). 
It shows at which rate each component of velocity is changing. It may not 
always point in the same direction as velocity (fig. 0.2). 
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Figure 0.2 - Acceleration vectors of a body following a curved trajectory when its 
speed (the magnitude of its velocity) remains constant (top) and when it changes 
continuously (bottom) 

Figure CC-0 o.c. 


0.3.2 Forces and moments 

A force expresses an effort exerted on a body. Force has a direction as well 

—> 

as a magnitude, and so it is expressed with a vector (typically noted F). 
The three main types of forces relevant to fluid mechanics are those due 
to pressure, those due to shear, and those due to gravity (the force due to 
gravity is called weight). 

When a force exerts at a position r away from a reference point, it exerts a 
torsion (or “twisting”) effort named moment. Like a force, a moment has a 
direction as well as a magnitude, and so is best expressed with a vector. This 
vector M is expressed as the cross-product of the arm r and the force F: 

M^rAF (0/3) 


0.3.3 Energy 

Energy, measured in joules 0), is in most general terms the ability of a body 
to set other bodies in motion. It can be accumulated or spent by bodies in a 
large number of different ways. The most relevant forms of energy in fluid 
mechanics are: 

Kinetic energy £k, accumulated as motion 

£k = (0/4) 

Work W, which is energy spent on displacing an object over a distance I 
with a force F: 


W = F-l 


(0/5) 
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Internal energy I stored as heat within the body itself. As long as no phase 
changes occurs, the internal energy I of fluids is roughly proportional 
to their absolute temperature T. 


0.4 Properties of fluids 


Beyond velocity, which is the primary unknown for us in fluid mechanics, 
there are a few other important fluid properties. 


0.4.1 Density 


The density p is defined as the inverse of mass-specific volume v, 

1 m 

where p the density (kgm“^); 

V is the specific volume (m^ kg“^); 
m is the considered mass (kg); 
and is the considered volume (m^). 


( 0 / 6 ) 


At ambient conditions air has a density of approximately pair = 1,2 kg m ^ 
that of water is almost a thousand times greater, at Pwater = 1 000kgm“^. 


0.4.2 Phase 

Fluids can be broadly classified into phases, which are loosely-defined sets 
of physical behaviors. Typically one distinguishes liquids which are fluids 
with large densities on which surface tension effects play an important role, 
from gases or vapors which have low densities and no surface tension effects. 
Phase changes are often brutal (but under specific conditions can be blurred 
or smeared-out); they usually involve large energy transfers. The presence of 
multiple phases in a flow is an added layer of complexity in the description 
of fluid phenomena. 


0.4.3 Temperature 

Temperature is a scalar property measured in Kelvins (an absolute scale). It 
represents a body’s potential for receiving or providing heat and is defined, 
in thermodynamics, based on the transformation of heat and work. 

A conversion between Kelvins and degrees Celsius by subtracting 273,15 
units: 


T{°C) = T(K) - 273,15 


(0/7) 


Although we can “feel” temperature in daily life, it must be noted that the 
human body is a very poor thermometer in practice. This is because we 
constantly produce heat, and we infer temperature by the power our body 
loses or gains as heat. This power not only depends on our own body 
temperature (hot water “feels” hotter when we are cold), but also on the 
heat capacity of the fluid (cold water “feels” colder than air at the same 
temperature) and on the amount of convection taking place (ambient air 
“feels” colder on a windy day). 
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In spite of these impressions, the fact is that the heat capacity of fluids 
is extremely high (cair » 1 kjkg“^ K, crater ~ 4kjkg“^ K). Unless very high 
velocities are attained, the temperature changes associated with fluid flow 
are much too small to be measurable in practice. 

0.4.4 Perfect gas model 

Under a specific set of conditions (most particularly at high temperature 
and low pressure), several properties of a gas can be related easily to one 
another. Their absolute temperature T is then modeled as a function of their 
pressure p with a single, approximately constant parameter R = pv/T: 

pv = RT 

- =RT 
P 

where R depends on the state and nature of the gas 0K“' kg“h; 
and p is the pressure (Pa), defined later on. 

This type of model (relating temperature to pressure and density) is called an 
equation of state. Where R remains constant, the fluid is said to behave as a 
perfect gas. The properties of air can satisfactorily be predicted using this 
model. Other models exist which predict the properties of gases over larger 
property ranges, at the cost of increased mathematical complexity. 

Many fluids, especially liquids, do not follow this equation and their tem¬ 
perature must be determined in another way, most often with the help of 
laboratory measurements. 

0.4.5 Speed of sound 

An important property of fluids is the speed at which information (in particu¬ 
lar, pressure changes due to the movement of an object) can travel within the 
fluid. This is equal to the average speed of molecules within the fluid, and it 
is called the speed of sound, noted c. 

It is important to quantify how fast the fluid is flowing relative to e speed 
of sound. For this, we define the Mach number [Ma] as the ratio of the local 
fluid speed V to the local speed of sound c: 

V 

[Ma] = - (0/10) 

c 

Since both V and c can be functions of space in a given flow, [Ma] may not be 
uniform (e.g. the Mach number around an airliner in flight is different at the 
nose than above the wings). Nevertheless, a single value is typically chosen 
to identify “the” Mach number of any given flow. 

It is observed that providing no heat or work transfer occurs, when fluids flow 
at [Ma] < 0,3, their density p stays constant. Density variations in practice 
can be safely neglected below [Ma] = 0,6. These flows are termed incom¬ 
pressible. Above these Mach numbers, it is observed that when subjected to 
pressure variations, fluids exert work upon themselves, which translates into 
measurable density and temperature changes: these are called compressibility 
effects, and we will study them in chapter 9. 

In most flows, the density of liquids is almost invariant - so that water flows 
are generally entirely incompressible. 


( 0 / 8 ) 

(0/9) 
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When the fluid is air (and generally within a perfect gas), we will show in 
chapter 9 (§9.3 p.l96) that c depends only on the absolute temperature: 

c = ^|YRT (0/11) 

in all cases for a perfect gas, 

where c is the local speed of sound (m s“^); 

p is a gas property, approx, constant (dimensionless); 
and T is the local temperature (K). 


0.4.6 Viscosity 

We have accepted that a fluid element can deform continuously under pres¬ 
sure and shear efforts; however this deformation may not always be “for 
free”, i.e. it may require force and energy inputs which are not reversible. 
Resistance to strain in a fluid is measured with a property named viscosity. 

In informal terms, viscosity is the “stickiness” of fluids, i.e. their resistance to 
being sheared: honey and sugar syrups are very viscous while water has low 
viscosity. 

Formally, viscosity is defined as the ratio between the shear effort applied 
to a fluid particle, and the time rate of its deformation (its strain rate). For 
example, if a fluid element is strained uniformly in the x-direction, the 
velocity Ux will be different everywhere inside the fluid element, and shear 
Tyx (the shear in the x-direction in the plane perpendicular to the y-direction) 
will occur. The viscosity p is then defined as: 



where is the viscosity (N s m~^ or Pa s); 

Tyx is the shear in the x-direction along the 
plane perpendicular to the y-direction (Pa); 
and Ux is the velocity in the x-direction (m s“^). 


This equation 0/12 can be generalized for all directions (i, j and k) as: 


k = 



(0/13) 


For most fluids, viscosity is the same in all directions and for all strain rates. 
This characteristic makes them Newtonian fluids. Viscosity is affected by 
temperature. We will study the effects of shear and viscosity in chapter 2. 
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0.5 Forces on fluids 


Fluids are subjected to, and subject their surroundings and themselves to 
forces. Identifying and quantifying those forces allows us to determine how 
they will flow. Three types of forces are relevant in fluid dynamics: gravity, 
pressure and shear. 

0.5.1 Gravity 

Gravity (here, the attraction effort exerted on fluids at a distance by the 
Earth) is expressed with a vector g pointing downwards. Within the Earth’s 
atmosphere, the magnitude g of this vector varies extremely slightly with 
altitude, and it may be considered constant at ^ = 9,81 ms“^. The weight 
force exerted on an object of mass m is then simply quantified as 

■^weight ~ (0/14) 


0.5.2 Pressure 

The concept of pressure can be approached with the following conceptual 
experiment: if a flat solid surface is placed in a fluid at zero relative velocity, 
the pressure p will be the ratio of the perpendicular force F± to the surface 
area A: 

p = ^ (0/15) 

where p is the pressure (N m~^ or pascals, 1 Pa = INm“^); 

Fj_ is the component of force perpendicular to the surface (N); 
and A is the surface area (m^). 

Although in si units pressure is measured in Pascals (1 Pa = INm“^), in 
practice it is often measured in bars (1 bar =1-10^ Pa). Ambient atmo¬ 
spheric pressure at normal altitude varies with the weather and is approxi¬ 
mately 1 bar. 

The pressure distribution within a fluid is related to their velocity distribution 
according to relations that we will study later. We shall then be looking for a 
pressure field p(^x,y,z,t), a function of space an time. 

The pressure p exerting on an element of fluid can be thought of as the 
time- and space-average of the perpendicular component of impact force 
of its molecules on its neighbors. It is strictly a macroscopic property, that 
is, it cannot be defined at a microscopic level (there is no such thing as the 
“pressure of a molecule”). In subsonic flows, it is a scalar property, meaning 
that for a given particle, it is the same in all directions. Pressure effects are 
explored in further detail in chapter 1. 

0.5.3 Shear 

In the same thought experiment as above, shear t expresses the efforts of a 
force parallel to a surface of interest: 



where t is the shear (N m ^ or Pa); 

F|| is the component of force parallel to the surface (N); 


(0/16) 
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and A is the surface area (m^). 


Contrary to pressure, shear is not a scalar, i.e. it can (and often does) take 
different values in different directions: on the flat plate above, it would have 
two components and could be represented by a vector r = Tyj. We will 
explore shear in further detail in chapter 2. 


0.6 Basic flow quantities 


A few fluid-flow related quantities can be quantified easily and are worth 
listing here. 

Mass flow is noted m and represents the amount of mass flowing through a 
chosen surface per unit time. When the velocity across the surface is 
uniform, it can be quantified as: 


m = pV±A (0/17) 

where m is the mass flow (kg s“^); 

p is the fluid density (kg s“h; 

A is the area of the considered surface (m^); 
and Vj_ is the component of velocity perpendicular to the surface (m s“*). 

Alternatively, one might consider instead the full fluid velocity and the 
area of a surface perpendicular to that velocity: 

m = pVAx. (0/18) 

where V is the flow speed (m s“^); 

and Ax is the area of a surface perpendicular to the flow velocity (m^). 

Volume flow is noted W and represents the volume of the fluid flowing 
through a chosen surface per unit time. Much like mass flow, when the 
velocity is uniform, it is quantified as: 


W = VxA = VAx = - (0/19) 

P 

where 'V is the volume flow (m^ s“h- 

Mechanical power is a time rate of energy transfer. Compression and 
expansion of fluids involves significant power. When the flow is in¬ 
compressible (see §0.8 below), the mechanical power W necessary to 
force a fluid through a fixed volume where pressure losses occur is: 

^ ~ f'net, pressure ‘ kfluid, average ~ ^ |/^plloss (0/20) 

where W is the power spent as work (W); 

and |Ap| is the pressure loss occurring due to fluid flow (Pa). 

Power as heat is also a time rate of energy transfer. When fluid flows 
uniformly and steadily through a fixed volume, its temperature T 
increases according to the net power as heat Q: 

Q = m CfluidAT (0/21) 

where Q is the power spent as heat (W); 

AT is the temperature change occurring in the fluid (K); 

and Cfiuid is the specific heat capacity of the fluid (JK“^ kg“^). 
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The heat capacity of fluids varies strongly according to the amount of 
work that they are performing. When no work is performed in a steady 
flow, the heat capacity is termed Cp. In fluids such as liquid water and 
air, this capacity is almost independent from temperature. 


0.7 Five equations 


The analytical and numerical resolution of problems in fluid mechanics is 
arched upon five important physical principles, which are sometimes called 
conservation principles. 


1. Mass conservation: 

Without nuclear reactions, the total amount of matter at hand in a 
given phenomenon must remain constant. This statement can be 
expressed as: 


^system ~ CSt 


dm 


system 


dt 


= 0 


( 0 / 22 ) 


2. Conservation of linear momentum: 

This is a form of Newton’s second law, which states that the net force 
Tnet = applying to any given system of mass m is equal to the time 
change of its linear momentum mV: 

Fnet = (mV) (0/23) 


3. Conservation of angular momentum: 

This is a different form of Newton’s second law, which states that 

—^ ^ 

the net moment Mnet, x = SMx applied on a system about a point X 
is equal to the time change of its total angular momentum about this 
same point r A mV: 

Mnet, X = -^ (7 A mV) (0/24) 


4. Conservation of energy: 

This equation, also known as the “first principle of thermodynamics”, 
states that the total amount of energy within an isolated system must 
remain constant: 


•iEisolated system 

dt 


= 0 


(0/25) 


5. Second principle of thermodynamics: 

By contrast with all four other principles, this is not a conservation 
statement. Instead this principle states that the total amount of en¬ 
tropy S within a constant-energy system can only increase as time 
passes; in other words, the entropy change of any system must be 
greater than the heat transfer Q divided by the temperature T at which 
it occurs: 

^^system ^ j, (0/26) 
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These are the only five important equations written in fluid mechanics. We 
usually apply these statements to our problem in either one of two ways: 

• We may wish to calculate the net (overall) change in fluid properties 
over a zone of interest so as to relate them to the net effect (forces, 
moments, energy transfers) of the fluid flow through the zone’s bound¬ 
aries. This method is called integral analysis and we will develop it 
in chapter 3. 

• We may instead wish to describe the flow through our zone of in¬ 
terest in an extensive manner, aiming to obtain vector fields for the 
velocity and pressure everywhere, at all times. This method is called 
differential analysis and we will develop it in chapter 4. 


0.8 Classification of fluid flows 


As we will see progressively, it is extremely difficult to obtain general so¬ 
lutions for fluid flow. Thus, whenever possible or reasonable, simplifying 
hypothesis are made about the behavior of any particular flow, that allow 
us to proceed with the analysis and obtain a reasonable, if inexact, specific 
solution. It is therefore a habit of fluid dynamicists to classify fluid flows in 
various categories, which are not necessarily incompatible. When approach¬ 
ing a given problem, we typically look out for the following characteristics: 

Time dependence Flows which do not vary with time are called steady. 
Steadiness is dependent on the chosen point of view: for example, 
the air flow around an airliner in cruise flight is seen to be steady 
from within the aircraft, but is obviously highly unsteady from the 
point of view of the air particles directly in the path of the airliner. 
Steadiness is also dependent on the selection of a suitable time frame: 
time-variations in a flow may be negligible if the time window used in 
the study is short enough. 

Mathematically, flows are steady when either the partial time derivative 
d/dt or total time derivative D/Dt (concepts we will study in chapter 4) 
of properties is zero: this considerably eases the search for solutions. 
It is important not to confuse steadiness (property fields constant in 
time) with uniformity (properties constant in space). 

To find out whether or not a flow is steady, compare instantaneous 
representations (photos, measurement readings, vector field represen¬ 
tations, etc.) taken at different times. If they are all identical, the flow 
is steady. There does not exist an analytical method, however, that 
allows us to predict whether a flow will be steady. 

Compressibility When the density of the fluid is uniform, the flow is said to 
be incompressible. The term is treacherous, because it refers to density, 
not pressure (incompressible flows almost always feature non-uniform 
pressure fields). 

To find out whether a gas flow is compressible, compute the Mach 
number (eq. 0/10 p.l5). Below 0,3 the flow is always incompressible. 
Compressibility effects can be reasonably neglected below [Ma] = 0,6. 
Unless very specific phenomena such as phase changes or extreme 
speeds occur, the flow of liquids is always incompressible. 
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Temperature distribution In a fluid flow, temperature changes can occur 
due to three phenomena: 

• Heat transfer from solid bodies; 

• Changes in pressure and density due to work being performed on 
the fluid (by moving solid bodies, or by the fluid itself); 

• Heat created through internal friction within the fluid. 

To find out whether a flow has uniform temperature (i.e. whether it is 
isothermal), three steps must be taken: 

1. Quantify the heat transfer from external bodies. If that is zero, 
the flow is adiabatic, meaning there is no heat transfer. 

2. Find out whether the flow is compressible. If not, the density 
changes are negligibly small and temperature changes do not 
occur due to compression or expansion. 

3. Quantify the mechanical power lost by the fluid as it flows through 
the domain of interest. This power is integrally lost to heat. 

Because the heat capacity of fluids is generally very high, temperature 
changes due to this third factor are usually negligibly small. This is 
assessed with an example in exercise 0.6 p.25. 

Turbulence One last characteristic that we systematically attempt to iden¬ 
tify in fluid flows is turbulence (or its opposite, laminarity). While 
laminar flows are generally very smooth and steady, turbulent flows 
feature multiple, chaotic velocity field variations in time and space. 

We shall first approach the concept of turbulence in chapter 5, and 
study it more formally in chapter 6. In the meantime, we may assess 
whether a flow will become turbulent or not using a non-dimensional 
parameter named the Reynolds number, noted [Re]: 

[Re] = ^ (0/27) 

b 

To find out whether a flow will become turbulent, quantify [Re] by 
using the fluid properties p and p, a representative fluid speed V, and a 
length L which is representative of the flow domain (for example, the 
stream-wise length of an obstacle). If the result is on the order of 10^ 
or more, the flow is very likely to become turbulent over the length L. 
By contrast, with [Re] of the order of 10^ or less, the flow is very likely 
to remain laminar over this length. 

This crude quantification has more meaning than it first appears. In par¬ 
ticular, it assesses the influence of viscosity, which acts to to attenuate 
non-uniformities and unsteadiness in fluid flow, effectively “damping” 
its kinematics. We will explore the meaning behind this procedure, 
and apply it to concrete examples, in chapters 5 and following. 


0.9 Limits of fluid mechanics 


Fluid mechanics is a complex discipline. 

It is easy to observe that flows as ordinary as sea waves crashing on a reef, 
water flowing down a river stream, or air blown into one’s hands, display 
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tremendous geometrical complexity. Even after choosing to describe only the 
movement of macroscopic fluid particles instead of individual molecules (and 
thereby avoiding studying thousands of billions of individual movements), 
we still need to describe a three-dimensional field of properties (pressure, 
temperature, etc.), one of which, velocity, is itself three-dimensional. 

Thus, even before we begin describing the exact problem and a procedure to 
obtain its solution, we know that the mere description of a solution can have 
tremendous complexity. 

Additionally, we have to admit that in practice much progress can be made 
in the field of fluid mechanics. For example, our weather forecasts have 
almost no value beyond one week, and aircraft manufacturers with bud¬ 
gets measured in billions of dollars still make extensive use of wind tunnel 
models - this despite our staggering continuous rate of progress in com¬ 
puting technology, and many decades of efforts dedicated to analytical fluid 
mechanics. 

In our present study of fluid mechanics, therefore, we shall proceed modestly, 
and will always take care to show the limits of our analysis. 
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Exercise sheet 0 (Important concepts) 

last edited April 1, 2018 


Except otherwise indicated, we assume that fluids are Newtonian, and that: 

Pwater = 1 000 kg m“^ patm. = 1 bar; Patm. = 1,225 kgm“^ Tatm. = ll,3°C;Patm. = 1,5 • lO'^Nsm”^; 
g = 9,81 ms“^. Air is modeled as a perfect gas (i?air = 287 JK“' kg“^; Yaa = 1,4; Cpair = 1 005Jkg“' K“'). 


0.1 Compressibility effects 

An aircraft is flying in air with density 0,9 kg and temperature -5 °C. Above which 
flight speed would you expect the air flow over the wings to become compressible? 


0.2 Pressure-induced force 

A 2 m by 2 m flat panel is used as the wall of a swimming pool (fig. 0.3). On the left side, 
the pressure is uniform at 1 bar. 



Figure 0.3 - Pressure distribution on a flat plate 


Figure CC-0 o.c. 


1. What is the pressure force exerted on the left side of the plate? 

On the right side of the plate, the water exerts a pressure which is not uniform: it 
increases with depth. The relation, expressed in pascals, is: 

pwater = 1,3 • 10^ - 9,81 • 10^ X Z (0/28) 


2. What is the pressure force exerted on the right side of the plate? 

[Hint: we will explore the required expression in chapter 1 as eq. 1/14 p.35] 
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0.3 Shear-induced force 

A fluid flows over a 3 m by 3 m flat horizontal plate, in the x-direction as shown in tig. 0.4. 
Because of this flow, the plate is subjected to uniform shear = 1,65 Pa. 

Z 

A 



Figure 0.4 - Shear force exerting on a plate 

1. What is the shear force applying on the plate? 


Figure CC-0 o.c. 


2. What would be the shear force if the shear was not uniform, but instead was a 
function of x expressed (in pascals) as Tzx = 1,65 - 0,01 x x^? 

[Hint: we will explore the required expression in chapter 2 as eq. 2/20 p. 53] 


0.4 Speed of sound 

White [13] Pl.87 

Isaac Newton measured the speed of sound by timing the interval between observing 
smoke produced by a cannon blast and the hearing of the detonation. The cannon is 
shot 8,4 km away from Newton. What is the air temperature if the measured interval 
is 24,2 s? What is the temperature if the interval is 25,1 s? 


0.5 Power lost to drag 

—^ .. 

A truck moves with constant speed T on a road, with V = 100kmh“k Because it 
experiences cross-wind, it is subjected to a drag Fd with Fd = 5 kN at an angle 9 = 20°, 
as shown in fig. 0.5. 

1. What is the power required to overcome drag? 

The drag force Fd is applying at a distance 0,8 m behind the center of gravity of the truck. 

—> 

2. What are the magnitude and the direction of the moment exerted by the drag Fd 
about the center of gravity? 






Figure 0.5 - Top view of a truck traveling at velocity V and subject to a drag force Fd 

Figure CC-0 o.c. 
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0.6 Go-faster exhaust pipe 

The engine exhaust gases of a student’s hot-rod car are flowing quasi-steadily in a 
cylindrical outlet pipe, whose outlet is slanted at an angle 9 = 25° to improve the good 
looks of the car and provide the opportunity for an exercise. 




e 


V 


top view 


Figure 0.6 - Exhaust gas pipe of a car. The outlet cross-section is at an angle 9 relative to the axis 
of the pipe. 

Figure CC-0 o.c. 

Photo cropped, mirrored and edited from an original CC-by-sa by kazandrew2 


The outlet velocity is measured at 15 ms and the exhaust gas density is 1,1 kg m 
The slanted outlet section area A is 420 cm^. 

1. What is the mass flow m through the pipe? 

2. What is the volume flow V of exhaust gases? 

Because of the shear within the exhaust gases, the flow through the pipe induces a 
pressure loss of 21 Pa (we will learn to quantify this in chapter 5). In these conditions, 
the specific heat capacity of the exhaust gases is Cpgases = 1 lOOJkg”^ K“k 

3. What is the power required to carry the exhaust gases through the pipe? 

4. What is the gas temperature increase due to the shear in the flow? 


0.7 Acceleration of a particle 

Inside a complex, turbulent water flow, we are studying the trajectory of a cubic fluid 
particle of width 0,1 mm. The particle is accelerating at a rate of 2,5 m s“^. 
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1. What is the net force applying to the particle? 

2. In practice, which types of forces could cause it to accelerate? 

0.8 Flow classifications 

1. Can an incompressible flow also be unsteady? 

2. Can a very viscous fluid flow in a turbulent manner? 

3. [more difficult] Can a compressible flow also be isothermal? 

4. Give an example of an isothermal flow, of an unsteady flow, of a compressible flow, 
and of an incompressible flow. 
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Answers 


0.1 If you adopt [Ma] = 0,6 as an upper limit, you will obtain V^ax = 709kmh“^ 
(eqs. 0/10 & 0/11 p.l6). Note that propellers, fan blades etc. will meet compressiblity 
effects far sooner. 


0.2 1) Fieft = 400kN (eq. 0/15 p.l7); 

0.3 1) Fi = 14,85 N (eq. 0/16 p.l7); 

0.4 26,7 °C & 5,6 °C. 

0.5 l)W = /drag • /truck = 130,5 kW; 

2) M = 11 r A /drag II = 1 368 N m, M = 


2) Fright = 480 kN (eq. 1/14 p.35). 
2) F 2 = 14,58 N (eq. 2/20 p.53). 


0 (points vertically upwards). 
-1368 / 


0.6 l)m = 0,2929kg s“^ (eq. 0/17 p.l8); 2) = 266,2Ls“^ (eq. 0/19 p.l8); 

3)W = 5,59W (eq. 0/20 p.l8); 4) AT = +0,0174K (eq. 0/21 p.l8), an 

illustration of remarks made in §0.8 p.21 regarding temperature distribution. 


0.7 1) Tnet = 2,5 • 10“^ N (eq 0/23 p.l9), such are the orders of magnitude involved in 

CFD calculations! 2) Only three kinds: forces due pres¬ 

sure, shear, and gravity. 


0.8 1) yes, 2) yes if [Re] is high enough, 3) yes (in very specific cases such as high 

pressure changes combined with high heat transfer or high irreversibility, therefore 
generally no), 4) open the cap of a water bottle and turn it upside down: you have 
an isothermal, unsteady, incompressible flow. An example of compressible flow 
could be the expansion in a jet engine nozzle. 
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Chapter 1 - Effects of pressure 
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These lecture notes are based on textbooks by White [13], Qengel & al.[16], and Munson & al.[18]. 


1.1 Motivation 


In fluid mechanics, only three types of forces apply to fluid particles: forces 
due to gravity, pressure, and shear. This chapter focuses on pressure (we will 
address shear in chapter 2), and should allow us to answer two questions: 

• How is the effect of pressure described and quantified? 

• What are the pressure forces generated on walls by static fluids? 



Video: pre-lecture briefing for 
this chapter 


by o.c. (CC-by) 
https://youtu.be/4MVnlPoxGqY 


1.2 Concept of pressure 


We approached the concept of pressure in the introductory chapter with the 
notion that it represented force perpendicular to a given flat surface (eq. 0/15), 
for example a flat plate of area A: 

P = ^ ( 1 / 1 ) 

To appreciate the concept of pressure in fluid mechanics, we need to go 
beyond this equation. 

1.2.1 The direction of pressure 

An important concept is that in continuum mechanics, the flat surface is 
imaginary. More precisely, a fluid is able to exert pressure on not only on 
solid surfaces, but also upon and within itself. In this context, we need to 
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rework eq. 1/1 so that now pressure is defined as perpendicular force per 
area on an infinitesimally small surface of fluid: 

p 

p ^ lira (1/2) 
A^O A 

Equation 1/2 may appear unsettling at first sight, because as A tends to zero, 
Fj_ also tends to zero; nevertheless, in any continuous medium, the ratio of 
these two terms tends to a single non-zero value. 

This brings us to the second particularity of pressure in fluids: the pressure 
on either side of the infinitesimal flat surface is the same regardless of its 
orientation. In other words, pressure has no direction: there is only one (scalar) 
value for pressure at any one point in space. 

Thus, in a fluid, pressure applies not merely on the solid surfaces of its 
container, but also everywhere within itself. We need to think of pressure as 

a scalar property Held p(^x,y,z,t)- 

1.2.2 Pressure on an infinitesimal volume 

While pressure has no direction, it may not have the same value everywhere 
in a fluid, and so the gradient (the space rate of change) of pressure may not 
be null. For example, in a static water pool, pressure is uniform in the two 
horizontal directions, but it increases along with depth. 

Instead of a flat plate, let us now consider an infinitesimally small cube within 
the fluid (fig. 1.1). Because the cube is placed in a scalar field, the pressure 
exerting on each of its six faces may be different. The net effect of pressure 
will therefore have three components: one for each pair of opposing faces. 
Further down in this chapter, we will express this effect with the help of the 
gradient of pressure Vp, a vector field. For now, it is enough to notice that 
pressure itself is quantified as a (one-dimensional) scalar field, and the effect 
of pressure on a submerged volume as a (three-dimensional) vector field. 



Figure 1.1 - The pressure on each face of an infinitesimal volume may have a 
different value. The net effect of pressure will depend on how the pressure varies in 
space. These changes are labeled dp|, in each of the i = x,y,z directions. 

figure CC-0 o.c. 









1.3 Pressure in static fluids 


1.3.1 Fluid statics 


—^ 

Fluid statics is the study of fluids at rest, i.e. whose velocity field V is 
everywhere null and constant: 


V = 0 



(1/3) 


We choose to study this type of problem now, because it makes for a con¬ 
ceptually and mathematically simple case with which we can start analyzing 
fluid flow and calculate fluid-induced forces. 


1.3.2 Shear in static fluids 

We will see in chapter 2 that shear efforts can be expressed as a function 
of viscosity and velocity: ? = (2/16). We will also see that for most 

ordinary fluids (termed Newtonian fluids) the viscosity p is independent from 
the velocity distribution. 

For now, it is enough to accept that in a static fluid (one in which ^ = 0 

since V[x,y^,t) = 0); there cannot be any shear f. Static fluids are unable to 
generate shear effort. 

This is a critical difference between solids and fluids. For example: 

• A static horizontal beam extruding from a wall carries its own weight 
because it is able to support internal shear; but this is impossible for a 
static fluid. 

• A solid brick resting on the ground exerts pressure on its bottom face, 
but because of internal shear, the pressure on its sides is zero. On the 
contrary, a static body of water of the same size and weight cannot 
exert internal shear. Because pressure has no direction, the body exerts 
pressure not only on the bottom face, but also laterally on the sides of 
its container. 


1.3.3 Pressure and depth 

We now wish to express pressure within a fluid as a function of depth. Let us 
consider, inside a static fluid reservoir, a small element of fluid. We choose 
an element with finite horizontal surface 5 but infinitesimal height dz, as 
shown in fig. 1.2. 

The element’s geometry is such that all of the lateral forces (which are due to 
pressure only) compensate each other. In the vertical direction, we are left 
with three forces, which cancel out because the element is not accelerating: 

Fp bottom Fp top + — 0 

Pbottom‘5 ~ Ptop‘5 ~ — 0 
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Figure 1.2 - A fluid element of (finite) horizontal surface S and infinitesimal height dz 
within a static fluid. 

figure CC-0 o.c. 


If we now rewrite the bottom pressure as pbottom = ptop + dp, we obtain: 


(ptop “t dp)5 ptop‘5 — b 

dp = jg (1/4) 


This equation 1/4 expresses the fact that the infinitesimal variation of pres¬ 
sure dp between the top and bottom surfaces of the element is due only to its 
own weight. If we evaluate this variation along the infinitesimal height dz, 
we obtain: 


dp m m 
dz " " d^^ 


dp 

dz 


= pg 


(1/5) 


Thus, equation 1/5, which is obtained from Newton’s second law applied to 
a static fluid, links the vertical variation in pressure to the fluid’s density p 
and the local gravitational acceleration g. This equation is so useful that it is 
sometimes called “the fundamental equation of fluid statics”. 

If our coordinate system is not aligned with the vertical direction (and thus 
with gravity), we will have to adapt equation 1/5 and obtain a more general 
relationship: 


dp 

dp 

Tz 


dp^ dp^ 


P9x 


P9y 

(1/6) 

P9z 


II 

^1 ^ 

(1/7) 


We can simplify the writing of this last equation by using the mathematical 
operator gradient (see also Appendix A2 p.217), defined as so: 


g _^d d 

dx dy dz 


( 1 / 8 ) 
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And thus, equation 1/7 is more elegantly re-written as: 


Vp = pg 


(1/9) 


This equation is read “the gradient of pressure is equal to the density times 


the gravitational acceleration vector”. It expresses concisely a very important 
idea: in a static fluid, the spatial variation of pressure is due solely to 
the fluid’s own weight. 

In chapter 4, we shall see that eq. 1/9 is a special case in a much more potent 
and complex relation called the Navier-Stokes equation (eq. 4/37 p.94). For 
now, it suffice for us to work with just pressure and gravity. 

1.3.4 Pressure distribution 

An immediate consequence stemming from equation 1/9 is that in a static 
fluid (e.g. in a glass of water, in a swimming pool, in a calm atmosphere), 
pressure depends solely on height. Within a static fluid, at a certain altitude, 
we will measure the same pressure regardless of the surroundings (fig. 1.3). 



Figure 1.3 - Pressure at a given depth (or height) in a static fluid does not depend on 
the environment. Here, as long as the fluid remains static, pA = pB = pc = pD- 


Figure CC-0 o.c. 


1.3.5 Hydrostatic pressure 

How is pressure distributed within static liquid water bodies? The density 
of liquid water is approximately constant: Pwater = 1 OOOkgrn”^. In a water 
reservoir, if we align the z-direction vertically downwards, equation 1/9 
becomes: 


Vp = pg 




1 000 X 9,81 = 9,81 • lO^Pam”^ = 9,81 • 10“^barm“^ 


( 1 / 10 ) 
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Therefore, in static water, pressure increases by approximately 0,1 bar/m 
as depth increases. For example, at a depth of 3 m, the pressure will be 
approximately 1,3 bar (which is the atmospheric pressure plus Az x 


1.3.6 Atmospheric pressure 


The density pair of atmospheric air is not uniform, since —unlike for water— 
it depends strongly on temperature and pressure. If we model atmospheric 
air as a perfect gas, once again orienting z vertically downwards, we can 
express the pressure gradient as: 


dz) 


— Pair 9 ~ P 


atm. 


TR 


( 1 / 11 ) 


This time, the space variation of pressure depends on pressure itself (and it is 
proportional to it). A quick numerical investigation for ambient temperature 
and pressure (1 bar, 15 °C) yields: 


dzj 


atm. ambient 


1 • 10^ X 


1 

288,15 


11,86Pam ^ = 


9,81 

X- 

287 

1,186 • 10“^ barm“^ 


This rate (approximately 0,1 mbar/m) is almost a thousand times smaller than 
that of water (fig. 1.4). 

Since the rate of pressure change depends on pressure, it also varies with 
altitude, and the calculation of pressure differences in the atmosphere is a 
little more complicated than for water. 

If we focus on a moderate height change, it may be reasonable to consider 
that temperature T, the gravitational acceleration g and the parameter R (gas 
constant) are uniform. In this admittedly restrictive case, equation 1/11 can 



Figure 1.4 - Variation of pressure as a function of altitude for water and air at the 
surface of a water reservoir. The gradient of pressure with respect to altitude is 
almost a thousand times larger in water than in air. 

Figure CC-0 o.c. 
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be integrated as so: 




RTcst 


dz 


exp 


RTcst. 


( 1 / 12 ) 


Of course, air temperature varies significantly within the atmosphere (at 
moderate altitudes the change with altitude is approximately -6Kkm“^). 
Adapting equation 1/12 for a uniform temperature gradient (instead of uni¬ 
form temperature) is the subject of exercise 1.8 p.43. 

In practice, the atmosphere also sees large lateral pressure gradients (which 
are strongly related to the wind) and its internal fluid mechanics are com¬ 
plex and fascinating. Equation 1/12 is a useful and convenient model, but 
refinements must be made if precise results are to be obtained. 


1.4 Wall pressure forces and buoyancy 


1.4.1 Pressure forces on plane surfaces 


When the pressure p exerted on a flat surface of area S is uniform, the 
resulting force F is easily easily calculated (f = pcst.5). In the more general 
case of a three-dimensional object immersed in a fluid with non-uniform 
pressure, the force must be expressed as a vector and obtained by integration: 

-^pressure ~ X dF = £pndS (1/13) 

where the 5-integral denotes an integration over the entire surface; 

and n is a unit vector describing, on each infinitesimal surface element d5, the 

direction normal to the surface. 



Video: when pres sure-induced 
forces in static fluids matter: 24 
hours of heavy tonnage transit 
through the Miraflores locks in 
Panama. Can you quantify the 
force applying on a single lock 
door? 

by Y:XyliboxFrance (styl) 
https ://youtu.be/ke4Rf9FrHnw 


Equation 1/13 is easily implemented in software algorithms to obtain numeri¬ 
cally, for example, the pressure force resulting from fluid flow around a body 
such as an aircraft wing. In our academic study of fluid mechanics, however, 
we will restrict ourselves to the simpler cases where the surface is perfectly 
flat (and thus n is everywhere the same). Equation 1/13 then becomes: 

^pressure ~ | dF = j p dS (1/14) 

Js Js 

for a flat surface. 


What is required to calculate the scalar E in eq. 1/14 is an expression ofp 
as a function of S. In a static fluid, this expression will be given by eq. 1/9. 
Typically in two dimensions x and y we re-write d5 as dS = dx dy and we 
may then proceed with the calculation starting from 


^ pressure 


= JJ P(x,y) dx dy 


(1/15) 
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Such a calculation gives us the magnitude of the pressure force, but not its 
position. This position can be evaluated by calculating the magnitude of 
moment generated by the pressure forces about any chosen point X. This 
moment Afx, using notation shown in fig. 1.5, is expressed as: 


Mx = 



rxF A p n dS 


(1/16) 


where ocf is a vector expressing the position of each infinitesimal surface relative to 
point X. 


Much like eq. 1/13 above, this eq. 1/16 is easily implemented in software 
algorithms but not very approachable on paper. In the simple case where we 
study only a flat surface, and where the reference point X is in the same plane 
as the surface, eq. 1/16 simplifies greatly and we can calculate the magnitude 
Afx as: 


Afx = 




rxF dF - 



rxF p dS 


for a flat surface, with X in the plane of the surface. 


(1/17) 


Once both fpressure and Afx pressure have been quantified, the distance Rxf 
between point X and the application point of the net pressure force is easily 
computed: 

(1/18) 

-^pressure 

While we restrict ourselves to simple cases where surfaces are flat, the 
process above outlines the method used by computational fluid dynamics 
(cfd) software packages to determine the position of forces resulting from 
fluid pressure on any arbitrary solid object. 




Figure 1.5 - Moment generated about an arbitrary point X by the pressure exerted 
on an arbitrary surface (left: perspective view; right: side view). The vector n is a 
convention unit vector everywhere perpendicular to the infinitesimal surface dS 
considered. 


Figure CC-0 o.c. 




1.4.2 Buoyancy and Archimedes’ principle 


Any solid body immersed within a fluid is subjected to pressure on its walls. 
When the pressure is not uniform (for example because the fluid is subjected 
to gravity, although this may not be the only cause), then the net force due 
to fluid pressure on the body walls will be non-zero. 

When the fluid is purely static, this net pressure force is called buoyancy. 
Since in this case, the only cause for the pressure gradient is gravity, the 
net pressure force is oriented upwards. The buoyancy force is completely 
independent from (and may or may not compensate) the object’s weight. 



Video: playing around with an 
air pump and a vacuum cham¬ 
ber 

by Y;Roobert33 (styl) 
https://youtu.be/ViuQKqUQlU8 


Since it comes from equation 1/9 that the variation of pressure within a fluid 
is caused solely by the fluid’s weight, we can see that the force exerted on 
an immersed body is equal to the weight of the fluid it replaces (that is to 
say, the weight of the fluid that would occupy its own volume were it not 
there). This relationship is sometimes named Archimedes’ principle. The 
force which results from the static pressure gradient applies to all immersed 
bodies: a submarine in an ocean, an object in a pressurized container, and 
of course, the reader of this document as presently immersed in the earth’s 
atmosphere. 





Figure 1.6 - Immersion in a static fluid results in forces that depend on the body’s vol¬ 
ume. They can evidenced by the removal of the fluid (for example in a depressurized 
semi-spherical vessel). 

Figure CC-0 o.c. 
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Exercise sheet 1 (Effects of pressure) 

last edited April 18, 2018 


These lecture notes are based on textbooks by White [13], Qengel & al.[16], and Munson & al.[18]. 


Except otherwise indicated, we assume that fluids are Newtonian, and that: 

Pwater = 1 000 kg m“^ patm. = 1 bar; Patm. = 1,225 kgm“^ Tatm. = ll,3°C;Patm. = 1,5 • 10“5Nsm“2; 
g = 9,81 ms“^. Air is modeled as aperfect gas (i?air = 287JK“^ kg“'; pair = 1,4; Cpair = 1005Jkg“^ K“^). 


1.1 Pressure in a fluid 

A small water container whose geometry is described in fig. 1.7 is filled with water. What 
is the pressure at the bottom of the container? 



Figure 1.7 - A small water container. 


Figure CC-0 o.c. 


1.2 Pressure measurement 


A tube is connected to a pressurized vessel as shown in fig. 1.8. The U-tube is filled with 
water. What is the pressure pint, in the vessel? 



13 cm 
4 cm 

12 cm 
19 cm 


Figure 1.8 - Working principle of a simple liquid tube manometer. The outlet is at atmospheric 
pressure patm.- 

Figure CC-0 o.c. 


what would be the height difference shown for the same internal pressure if mercury 
(Pmercury = 13 600 kg m“^) was used instead of water? 
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1.3 Straight water tank door 
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A water tank has a window on one of its straight walls, as shown in figure 1.9. The 
window is 3 m high, 4 m wide, is positioned 0,4 m above the bottom of the tank. 



Figure 1.9 - An aquarium tank with a window installed on one of its walls 


Figure CC-0 o.c. 


1. what is the force applying on the window due to the pressure exerted by the 
water? 

2. At which height does this force apply? 


1.4 Water canal access door 

A water canal used in a laboratory is filled with stationary water (fig. 1.10). An obser¬ 
vation window is installed on one of the walls of the canal, to enable observation and 
measurements. The window is hinged on its bottom face. 

The hinge stands 1,5 m below the water surface. The window has a length of 0,9 m and 
a width of 2 m. The walls of the canal are inclined with an angle 9 = 60° relative to 
horizontal. 



Figure 1.10 - A door installed on the wall of a water canal. The water in the canal is perfectly 
still. 

Figure CC-0 o.c. 

1. Represent graphically the pressure of the water and atmosphere on each side of 
the window. 

2. What is the moment exerted by the pressure of the water about the axis of the 
window hinge? 






























3. If the same door was positioned at the same depth, but the angle 6 was decreased, 
would the moment be modified? (briefly justify your answer, e.g. in 30 words or 
less) 


1.5 Water lock 


A system of lock doors is set up to allow boats to travel up the side of a hill (figure 1.11). 

door A 



— 

top view 

— 

door A 2 5 m 

— 

— 


— 

— 

ill 

^ - 

- ^ - 



Figure 1.11 - Outline schematic of a water canal lock. 


Figure 1 CC-0 o.c.; figure 2 CC-by-sa o.c. 


We are studying the force and moment exerted by the water on the lock door labeled A. 

At this instant, the water levels are as shown in figure 1.11: 2 m in the lower canal, 5 m 
in the lock, and 9 m in the upper canal. The door is 2,5 metres wide. 
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1. Sketch the distribution of the pressure exerted by the water and the atmosphere 
on both sides of door A. 

2. What is the net force exerted by the water on door A? 

3. At which height is this force exerting? 


1.6 Buoyancy force on a tin can 


A student contemplates a tin can of height 10 cm and diameter 7 cm. 

1. If one considers that the atmospheric density is uniform, what is the buoyancy 
force generated by the atmosphere on the can when it is positioned vertically? 

2. What is the force generated when the can is immersed in water at a depth of 20 cm? 
At a depth of 10 m? 

3. What is the buoyancy force generated when the can is immersed in the water in a 
horizontal position? 


Instead of uniform density, we now wish to calculate the atmospheric buoyancy under 
the hypothesis of uniform temperature (room temperature 20 °C). 


4. Starting from equation 1/9: Vp = pg, show that when the temperature Test, is 
assumed to be uniform, the atmospheric pressure p at two points 1 and 2 separated 
by a height difference Az is such that: 


p2 

— = exp 

Pi 


RZ 


CSt. 


( 1 / 12 ) 


5. What is the buoyancy generated by the room atmosphere? 


1.7 Buoyancy of a barge 

A barge of very simple geometry is moored in a water reservoir (fig. 1.12). 



Figure 1.12 - Basic layout of a barge floating in water. 

Figure CC-0 o.c. 

1. Sketch the distribution of pressure on each of the immersed walls of the barge (left 
and right sides, rear , bottom and slanted front). 

2. What is the force resulting from pressure efforts on each of these walls? 

3. What is the weight of the barge? 
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1.8 Atmospheric pressure distribution 


The integration we carried out in §1.3.6 p.34 to model the pressure distribution in the 
atmosphere was based on the hypothesis that the temperature was uniform and constant 
(T = Test.). In practice, this may not always be the case. 

1. If the atmospheric temperature decreases with altitude at a constant rate (e.g. 
of -7 Kkm“b> how can the pressure distribution be expressed analytically? 

A successful fluid dynamics lecturer purchases an apartment at the top of the Burj 
Khalifa tower (800 m above the ground). Inside the tower, the temperature is controlled 
everywhere at 18,5 °C. Outside, the ground temperature is 30 °C and it decreases linearly 
with altitude (gradient: -7Kkm“b. 

A door is opened at the bottom of the tower, so that at zero altitude the air pressure 
(1 bar) is identical inside and outside of the tower. 

For the purpose of the exercise, we pretend the tower is entirely hermetic (meaning air is 
prevented from flowing in or out of its windows). 

2. What is the pressure difference between each side of the windows in the apartment 
at the top of the tower? 


1.9 Reservoir door 

Munson & al. [18] 2.87 

A water reservoir has a door of width 3 m which is held in place with a horizontal cable, 
as shown in fig. 1.13. The door has a mass of 200 kg and the friction in the hinge is 
negligible. 



Figure 1.13 - A sealed, hinged door in a water reservoir. The width across the drawing (towards 
the reader) is 3 m 

Figure CC-0 o.c. 


1. What is the force in the cable? 

2. If the water height was decreased, how would this force be modified? (briefly 
justify your answer, e.g. in 30 words or less) 
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non-examinable. CC-0 o.c. 


1.10 Atmospheric buoyancy force 

If the atmospheric density is considered uniform, estimate the buoyancy force exerted 
on an Airbus A380, both on the ground and in cruise flight (pcmise = 0,4kgm“^,Tcruise = 
-40 °C). 



Figure 1.14 - Airbus A380-800 

Drawing CC-by-sa by Julien Scavini 


Length overall 

72,73 m 

Wingspan 

79,75 m 

Height 

24,45 m 

Wing area 

845 m^ 

Aspect ratio 

7,5 

Wing sweep 

33,5° 

Maximum take-off weight 

560 000 kg 

Typical operating empty weight 

276 800 kg 


Table 1.1 - Characteristics of the Airbus A380-800. 


1.11 Lock with diagonally-mounted doors 

^ non-examinable. CC-0 o.c. 

The doors of the lock studied in exercise 1.5 p.41 are replaced with diagonally-mounted 
doors at an angle a = 20°, mounted such that no bending moment is sustained by the 
hinges (fig. 1.15). 

What is the force exerted by door “A” on the other door? 



top view 


Figure 1.15 - Lock doors mounted with an angle relative one to another. The angle relative to the 
case in exercise 1.5 is a = 20°. The width of the canal is still 5 m. 

Figure CC-0 o.c. 
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Answers 


1.1 pA = patm. + 0,039 bar « 1,039 bar. 

1.2 l)pinside = patm. + 0,0157 bar « 1,0157 bar; 2) Az 2 = 1,1765 cm. 

1.3 1) Pnet ~ P9^ (-^max-bmax ~ 2’^rnax) ~ 176,6 kN; 

2) Mnet, bottom hinge = ppb (|Zmaxbmax “ jb^^x) = 176,6 kN m SO the force exerts at 

R = = 1 m above the bottom hinge. 

1.4 2) Mnet = 7,79 kN m; 3) observe the equation used to calcu¬ 

late Mnet to answer this question. If needed, ask for help in class! 

1.5 2) Pnet, door A = ppb = +686,7 kN (positive in right direction) 

3) Mnet, door A, about bottom axis ~ Pnet 2 ~ 2,4689 AIN m (positive clockwise): Pnet, door A 
exerts at R 2 = 3,595 m from the bottom. 

1.6 1) Pvertical = 4,625 mN upwards; 2) Pyertical = 3,775 N in both cases; 

3) There is no change; 4) See §1.3.6; 

5) Pvertical = 4,487 mN upwards (and so in question 1 we were off by 3 %). 

1.7 2) Prear = 0,2453 MN, Pside = 2,1714 MN, Pbottom = 13,734 MN, Pfront = 0,3468 MN; 

3) Pbuoyancy = 13,979 MN (1 425 t). 

1.8 l)fe = (l + ^)“. 

1.9 Mw.,er = I (-J[f = 0,4186 MNm; so, feable = 81,13 kN. 

1.10 Assuming a volume of approx. 2 600 m^, we obtain approx. 30,9 kN on the ground, 
10,4 kN during cruise. 

1.11 The moment is brought to zero by an inter-door force Psideways = 1,068 MN (per¬ 
pendicular to the canal axis). 
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These lecture notes are based on textbooks by White [13], Qengel & al.[16], and Munson & al.[18]. 


2.1 Motivation 


In fluid mechanics, only three types of forces apply to fluid particles: forces 
due to gravity, pressure, and shear. This chapter focuses on shear, and should 
allow us to answer two questions: 

• How is the effect of shear described and quantified? 

• What are the shear forces generated on walls by simple flows? 


2.2 Concept of shear 


We approached the concept of shear in the introductory chapter with the 
notion that it represented force parallel to a given flat surface (eq. 0/16), for 
example a flat plate of area A: 



( 2 / 1 ) 


Like we did with pressure, to appreciate the concept of shear in fluid mechan¬ 
ics, we need to go beyond this equation. 


2.2.1 The direction of shear 


Already from the definition in eq. 2/1 we can appreciate that “parallel to a 
flat plate” can mean a multitude of different directions, and so that we need 
more than one dimension to represent shear. Furthermore, much in the same 
way as we did for pressure, we do away with the flat plate and accept that 
shear is a Geld, i.e. it is an effort applying not only upon solid objects but also 
upon and within fluids themselves. We replace eq. 2/1 with a more general 
definition: 



A-^O A 


( 2 / 2 ) 



Video: pre-lecture briefing for 
this chapter 

by o.c. (CC-by) 
https://youtu.be/BKKXJWgLwJg 
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Video: cloud movements in a 
time-lapse video on an interest¬ 
ing day are evidence of a highly- 
strained atmosphere: pilots and 
meteorologists refer to this as 
wind shear. 

by Y:StormsFishingNMore (styl) 
https://youtu.be/LjWeYPEmCk8 



Contrary to pressure, shear is not a scalar, i.e. it can (and often does) take 
different values in different directions. At a given point in space we represent 
it as a vector r = (r^, r^, Tj, j, and in a fluid, there is a shear vector Geld: 


I Tx ^ 

^(x,y,z,t) = I 

' ' (x,y,z,t) 


(2/3) 


2.2.2 Shear on an infinitesimal volume 


Describing the changes in space of the shear vector field requires another 
mathematical dimension (called order). 

Instead of a flat plate, let us consider an infinitesimally small cube within the 
fluid (fig. 2.1). Because the cube is immersed inside a vector field, the shear 
vector exerting on each of its six faces may be different. 

In order to express the efforts on any given face, we express a component 
of shear with two subscripts, the first indicating the direction normal to the 
surface of interest, and the second indicating the direction of the effort. For 
example, Txy represents the shear in the y-direction on a surface perpendicular 
to the x-direction. On this face, the shear vector would be: 

'^xj ~ '^xx '^xy L)cz (2/4) 

~ ^xx^ '^xyj ^xz^ (2/5) 

where the subscript xj indicates all of the directions (J = x,y,z) on a face perpendicu¬ 
lar to the x-direction. 


In eq. 2/4, the reader may be surprised to see the term Txx appear — a shear 
effort perpendicular to the surface of interest. This is because the faces of 
the infinitesimal cube studied here (shown in fig. 2.1) are not solid. They are 
permeable, and the local velocity on each one may (in fact, must, if there is 
to be any flow) include a component of velocity through the face of the cube. 
Thus, there is no reason for the shear effort, which is three-dimensional, to be 
aligned along each flat surface. As the fluid travels across any face, it can be 
sheared (which results in strain) in any arbitrary direction, regardless of the 
local pressure — and thus shear can and most often does have a component 
(tu) perpendicular to an arbitrary surface inside a fluid. 



Figure 2.1 - Shear efforts on a cubic fluid particle (with only the efforts on the visible 
faces 1 to 3 represented). The shear tensor has six members of three components 
each. 

Figure CC-0 o.c. 
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Now, the net shear effect on the cube will have eighteen components: one 
tree-dimensional vector for each of the six faces. Each of those components 
may take a different value. The net shear could perhaps be represented as 
en entity —a tensor— containing six vectors Ti,T 2 ,Ts. .. t^. By convention, 
however, shear is notated using only three vector components: one for each 
pair of faces. Shear efforts on a volume are thus represented with a tensor 
Geld Tif. 



( 2 / 6 ) 


In this last equation 2 / 6 , each of the nine components of the tensor acts as the 
container for two contributions: one for each of the two faces perpendicular 
to the direction expressed in its first subscript. 

So much for the shear effort on an element of fluid. What about the net foree 
due to shear on the fluid element? Not every element counts: part of the 
shear will accelerate (change the velocity vector) the particle, while part of it 
will merely strain (deform) the particle. Quantifying this force thus requires 
making a careful selection within the eighteen components of Tij. We may 
start with the x-direction, which consists of the sum of the component of 
shear in the x-direction on each of the six cube faces: 

—^ 

fshear x ~ 3 ~ ^6^zx 6 

+S2'^yx 2 - SsTyx 5 

'^^I'^xx 1 ~ ^i'^xx 4 (2/7) 

Given that S 3 = Se = dx dy, that S 2 = S 5 = dx dz and that Si = S 4 = dz dy, 
this is re-written as: 


f^shear x — dx dy(Tzx 3 ~ '^zx &) 

+ dx dz{Tyx 2 ~ '^yx 5) 

-h dz dy{Txx 1 ~ '^xx 4) (2/8) 


In the same way we did with pressure in chapter 1 (§1.3.3 p.31), we express 
each pair of values as derivative with respect to space multiplied by an 
infinitesimal distance: 


dr. 


fshear x = dx dy ( dz-^ ) -h dx dz ( dy-^ 1 - 1 - dz dy I dx 


= d^ 


dz 

dUx dr^ 


dz 


dr,, 


dy 


dzx 


\ dx 




yx ^ dr XX 


dy 


dx 


(2/9) 
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If we make use of the operator divergent (see also Appendix A2 p.217), 
written V • : 


V- = 
V ■! = 

V-Aij = 


dx 

dA, 

dx 

dAx 


d - 
dA 

+ 


y 


dx 

dAxi 

dx 

dAx-. 


+ 


dy 

dA^x 


dy 


+ 


dA 


dy 


dx 


+ 


dA, 


dy 


dA, 


+ 


+ 

+ 

+ 


dz 

dA^x 

dz 

dz 
dA,, 
dz / 



f y-^\ 

= 

v-4 

\ V • Ai^ / 


( 2 / 10 ) 

( 2 / 11 ) 


( 2 / 12 ) 


we can re-write eq. 2/9 and see that the net shear force in the x-direction 
is equal to the particle volume times the divergent of the shear in the x- 
direction: 


f shear x ~ d'V V • Tjx 


(2/13) 


The y- and z-direction are taken care of in the same fashion, so that we can 
gather up our puzzle pieces and express the force per volume due to shear as 
the divergent of the shear tensor. 



^ Pshear x \ 

IV •^-.i \ 

'f'shear “ 

Pshear i/ 1 “ 

IV 


\ -f'shear z / 

IV-41 / 


dA' V • fij 


(2/14) 


This exploration goes beyond the theory required to to through this chapter 
(we will come back to it when we start concerning ourselves with the dynam¬ 
ics of fluid particles in chapter 4, where the divergent of shear will make part 
of the glorious Cauchy equation). It suffices for now to sum up our findings 
as follows: 

• Shear at a point in space has three components — it is a vector field; 

• The effect of shear on a volume of fluid has eighteen components - it 
is a second-order tensor field; 

• The net force due to shear on a volume of fluid, expressed using the 
divergent of the shear tensor, has three components — it is a vector 
field. 


2.3 Slip and viscosity 


2.3.1 The no-slip condition 

We observe that whenever we measure the velocity of a fluid flow along 
a solid wall, the speed tends to zero as we approach the wall surface. In 
other words, the fluid “sticks’ to the surface regardless of the overall faraway 
flow velocity. This phenomenon is called the no-slip condition and is of 
paramount importance in fluid mechanics. One consequence of this is that 
fluid flows near walls are dominated by viscous effects (internal friction) due 
to the large velocity gradients. 
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2.3.2 Viscosity 


When a solid wall is moved longitudinally within a fluid, the fluid generates 
an opposing friction force through viscous effects (fig. 2.2). We call viscosity ^ 
the ratio between the fluid velocity gradient and the shear effort. For example, 
the norm of the shear Txy on a surface perpendicular to the x-direction, in 
the y-direction, can be expressed as: 


dun 


\Txy\ \ = A'- 


dx 



In the general case, viscosity is defined as the (scalar) ratio between the norm 
of shear and the corresponding strain rate. The strain rate corresponding 
to the shear in the j-direction is the rate of change in the /-direction of the 
velocity in the j-direction (dujldi): 


m 

( 2 / 15 ) 

dui 


( 2 / 16 ) 


in which the subscript i is an arbitrary direction (x, y or z) and j is the direction 

following it in order (e.g. j = z when i = y); 

and where y is the viscosity (or dynamic viscosity) (Pas). 

Viscosity p is measured in Pa s, which is the same as N s m“^ or kg m“^ s“k 
It has historically been measured in poise (1 poise = 0,1 Pa s). 






Figure 2.2 - Any velocity gradient in the flow results in a shear force F|| in the 
direction y. The ratio between the shear and the velocity gradient is called viscosity. 

Figure CC-0 o.c. 


















The concept of kinematic viscosity v (Greek letter nu) is sometimes used 
instead of the dynamic viscosity; it is defined as 



(2/17) 


where v is measured in s ^ 


2.3.3 Newtonian Fluid 


Fluids for which p is independent from ^ are called Newtonian fluids. 

Most fluids of interest in engineering fluid mechanics (air, water, exhaust 
gases, pure gases) can be safely modeled as Newtonian fluids. Their viscos¬ 
ity p varies slightly with pressure and mildly with temperature - in our study 
of fluid mechanics, we will not take these dependencies into account. 

The values of viscosity vary very strongly from one fluid to another: for ex¬ 
ample, honey is roughly ten thousand times more viscous than water, which 
is roughly a hundred times more viscous than ambient air. The viscosities of 
various fluids are quantified in fig. 5.11 p.ll9. 

Oil-based paint, blood and jelly-based fluids are strongly non-Newtonian; 
they require more complex viscosity models (fig. 2.3). 



Figure 2.3 - Various possible viscosity characteristics of fluids. Those for which g is 
independent of duj/di are called Newtonian. 


Figure CC-0 o.c. 




2.4 Wall shear forces 


The calculation of forces due to shear is very similar to that which we used 
in the previous chapter with pressure. 

When the shear ? exerted on a flat surface of area S is uniform, the resulting 
force F is easily easily calculated {F = rcst.5). In the more general case of a 
three-dimensional object immersed in a fluid with non-uniform shear, the 
force must be expressed as a vector and obtained by integration: 

-fshear ~ I dF = fsurface (2/18) 

Much like equation 1/13 in the previous chapter, eq. 2/18 is easily imple¬ 
mented in software algorithms to obtain numerically, for example, the force 
resulting from shear due to fluid flow around a body such as an aircraft wing. 
In our academic study of fluid mechanics, however, we will restrict ourselves 
to the simpler cases where the surface is perfectly flat and the shear has 
uniform direction. The force in eq. 2/18, exterting in the /-direction due to 
shear on a surface perpendicular to the j-direction, then becomes: 

f shear 11 ^5 (2/19) 


What is required to calculate the scalar F in eq. 2/19 is an expression of r as 
a function of S. We proceed like we did in the previous chapter, splitting d5 
as dS = di dk before proceeding with the calculation starting from; 


(?■ 

II 

(2/20) 


(2/21) 

where i is the direction of the force; 

j is the direction perpendicular to the flat surface; 
and k is the third (orthonormal) direction. 



The above expression 2/21 is perhaps easier to read when it is developed. 
For example, the shear force fshear^i exerting on a plate perpendicular to the 
y-direction is: 


Fshearyi 


^ PshesLTyx 

0 

\ -Fshearj/z 


with 


^sheart/x ~ F 


Ilf. 


cbcdz 


^ shearyz 


-// 


dw 

dy 


dxdz 


Just like for pressure, the moment Mxj generated about a point X by the shear 
forces in the /-direction on a plane surface perpendicular to the J-direction 
can be calculated as 


M- 


X/ 


— I dAfx — I Txf A dF — I rxF A Zy dS 
Us Us Us 


( 2 / 22 ) 
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If the point X is in the plane of the surface of interest, this can be expressed 
as: 

^Xj = P JJ' fXF-^ didk (2/23) 

Mxy = 


54 



Exercise sheet 2 (Effects of shear) 

last edited April 30, 2017 


These lecture notes are based on textbooks by White [13], Qengel & al.[16], and Munson & al.[18]. 


Except otherwise indicated, we assume that fluids are Newtonian, and that: 

Pwater = 1 000 kg m“^ patm. = 1 bar; Patm. = 1,225 kgm“^ Tatm. = ll,3°C;patm. = 1,5 • 10‘^Nsm“2; 
g = 9,81 ms“^. Air is modeled as aperfect gas (i?air = 287JK“^ kg“'; pair = 1,4; Cpair = 1005Jkg“^ K“^). 


2.1 Elow in between two plates 

^ Munson & al [18] Ex 1.5 

A fluid is forced to flow between two stationary plates (fig. 2.4). We observe that the flow 
is laminar (smooth and fully steady), with a univorm velocity profile u = f{y) which is 
linked to the average fluid velocity i4verage by the relationship: 


(2/25) 

where y is measured from the middle of the gap; 
and H is half of the gap length. 


3 

/ y V 

^ “ “Average 

'-(s) 



Figure 2.4 - Velocity distribution for laminar flow in between two plates, also known as Couette 
flow. 

Figure CC-0 o.c. 


The fluid viscosity is 2Nsm the average velocity is 0,6ms ^ and the two plates 
are 10 mm apart. 

1. What is the shear effort plate generated on the lower plate? 

2. What is the shear effort Tyx in the middle plane of the flow? 
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2.2 Friction on a plate 

A plate the size of an A4 sheet of paper (210 mm X 297 mm) is moved horizontally at 
constant speed above a large flat surface (fig. 2.5). We assume that the velocity profile of 
the fluid betweeen the plate and the flat surface is entirely uniform, smooth, and steady. 



Figure 2.5 - A plate moved horizontally across a flat surface. 


Figure CC-0 o.c. 


1. Express the force Fy^ due to shear on the plate as a function of its velocity {7piate> 
the gap height H, and the properties of the fluid. 

2. The plate speed is f/piate = 1 m s“^ and the gap height is H = 5 mm. What is the 
shear force if the fluid is air (patm. = 1,5 • 10“^ N s m“^), and if the fluid is honey 
(/^honey ~ 40 N S m )? 

3. If a very long and thin plate with the same surface area was used instead of the 
A4-shaped plate, would the shear force be different? (briefly justify your answer, 
e.g. in 30 words or less) 


2.3 Viscometer , , 

Qengel & al. [16] 2-78 

An instrument designed to measure the viscosity of fluids (named viscometer) is made of 
two coaxial cylinders (fig. 2.6). The inner cylinder is immersed in a liquid, and it rotates 
within the stationary outer cylinder. 


The two cylinders are 75 cm tall. The inner cylinder diameter is 15 cm and the spacing 
is 1 mm. 


When the inner cylinder is rotated at 300 rpm, a friction-generated moment of 0,8 N m is 
measured. 


1. If the flow in between the cylinders corresponds to the simplest possible flow case 
(steady, uniform, fully-laminar), what is the viscosity of the fluid? 

2. Would a non-Newtonian fluid induce a higher moment? (briefly justify your 
answer, e.g. in 30 words or less) 
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[Note: in practice, when the inner cylinder is turned at high speed, the flow displays 
mesmerizing patterns called Taylor—Couette vortices, the description of which is much more 
complex!] 





15 cm 1 mm 


Figure 2.6 - Sketch of a cylinder viscometer. The width of the gap has been greatly exaggerated 
for clarity. 

Figure CC-0 o.c. 


2.4 Boundary layer 


White [13J PI.56 


A laminar fluid flow occurs along a wall (fig. 2.7). Close to the wall (y < 5), we observe 
that viscous effects dominate the mechanics of the flow. This zone is designated boundary 
layer. The speed can then be modeled with the relation: 


u = U sin 



in which U is the flow speed far away from the wall. 


(2/26) 



u 




Figure 2.7 - Velocity profile across the boundary layer. 


Figure CC-0 o.c. 


The fluid is helium at 20 °C and 1 bar (9-10 ^ N s m ^); measurements yield U = 10,8 m s ^ 
and (5=3 cm. 

1. What is the shear effort Tyx on the wall? 

2. At which height yi above the surface will the shear effort be half of this value? 

3. What would be the wall shear if helium was replaced with water (1 • 10“^ kg m“^ s“^)? 
































2.5 Clutch 


Qengel & al. [16] 2-74 

Two aligned metal shafts are linked by a clutch, which is made of two disks very close 
one to another, rotating in the same direction at similar (but not identical) speeds. The 
disk diameters are both 30 cm and the gap between them is 2 mm ; they are submerged 
in Sae30w oil with viscosity 0,38 kg m“^ s“k 



Figure 2.8 - Sketch of the two disks constituting the clutch. The gap width has been exaggerated 
for clarity. 

Figure CC-0 o.c. 


The power shaft rotates at 1450 rpm, while the powered shaft rotates at 1 398 rpm. We 
consider the simplest possible flow case (steady, laminar) in between the two disks. 

1. What is the moment imparted by one disk to the other? 

2. How would the moment change if the radius of each disk was doubled? 

3. What is the transmitted power and the clutch efficiency? 

4. Briefly (e.g. in 30 words or less) propose one reason why in practice the flow in 
between the two disks may be different from the simplest-case flow used in this 
exercise. 
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Answers 

2.1 1) Tyx\y=—H ~ 720Nm ; 2) Tyx\y=o = ONm 

2.2 1) Fyz = LiL 2 /i^ = 1,87 • 10““^N for air, and SOON for honey(!). 

2.3 u =—^^ = 1,281 • 10“^ Ns 

^ 2n<j)R\H 

2.4 1) ryx,wall = yufy = 5,09 • 10“^ N^ for helium; 

2) yi = |(5 = 2 cm; 3) = 0,565 N m~^ for water. 

2.5 1) M = f = 0,8228 N m; 3) = 120 W; //dutch = | = 

96,4 % (remember this is a very low-power, low-relative speed, laminar-flow case). 
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These lecture notes are based on textbooks by White [13], Qengel & al.[16], and Munson & al.[18]. 


3.1 Motivation 


Our objective for this chapter is to answer the question “what is the net effect 
of a given fluid flow through a given volume?”. 

Here, we develop a mass, momentum and energy accounting methodology to 
analyze the flow of continuous medium. This method is not powerful enough 
to allow us to describe extensively the nature of fluid flow around bodies; 
nevertheless, it is extremely useful to quantify forces, moments, and energy 
transfers associated with fluid flow. 

3.2 The Reynolds transport theorem 


3.2.1 Control volume 

Let us begin by describing the flow which interests us a a generic velocity 

—> —> 

field V = {u,v,-w) which is a function of space and time (V = f{x,y,z,t)). 

Within this flow, we are interested in an arbitrary volume named control 
volume (CV) which is free to move and change shape (fig. 3.1). We are going 
to measure the properties of the fluid at the borders of this volume, which 
we call the control surface (CS), in order to compute the net effect of the flow 
through the volume. 

At a given time, the control volume contains a certain amount of mass which 
we call the system (sys). Thus the system is a fixed amount of mass transiting 
through the control volume at the time of our study, and its properties 
(volume, pressure, velocity etc.) may change in the process. 

All along the chapter, we are focusing on the question: based on measured 
fluid properties at some point in space and time (the properties at the control 
surface), how can we quantify what is happening to the system (the mass 
inside the control volume)? 



Video: pre-lecture briefing for 
this chapter, part 1/2 

by o.c. (CC-by) 
https://youtu.be/lLXlFVtPoCY 
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control volume 



Figure 3.1 - A control volume within a flow. The system is the amount of mass 
included within the control volume at a given time. At a later time, it may have left 
the control volume, and its shape and properties may have changed. The control 
volume may also change shape with time, although this is not represented here. 

Figure CC-0 o.c. 


3.2.2 Rate of change of an additive property 

In order to proceed with our calculations, we need a robust accounting 
methodology. We start with a “dummy” fluid property B, which we will later 
replace with physical variables of interest. 

Let us thus consider an arbitrary additive property B of the fluid. By the term 
additive property, we mean that the total amount of property is divided if the 
fluid is divided. For instance, this is true of mass, volume, energy, entropy, 
but not pressure or temperature. 

The specific (i.e. per unit mass) value of B is designated b = B/m. 

We now wish to compute the variation of a system’s property B based on 
measurements made at the borders of the control volume. We will achieve 
this with an equation containing three terms: 

• The time variation of the quantity B within the system is measured 
with the term This may represent, for example, the rate of change 
of the fluid’s internal energy as it travels through a jet engine. 

• Within the control volume, the enclosed quantity Bqv can vary by 
accumulation (for example, mass may be increasing in an air tank fed 
with compressed air): we measure this with the term 

• Finally, a mass flux may be flowing through the boundaries of the 
control volume, carrying with it some amount of B every second: we 
name that net flow out of the system Bnet = 5out ~ Bin. 
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We can now link these three terms with the simple equation: 


dB 


sys 


dt 


dBcv 

dt 


^net 


(3/1) 


the rate of change 
of B for the system 


the rate of change 
of B within the + 
control volume 


the net flow of B 
through the boundaries 
of the control volume 


Since B may not be uniformly distributed within the control volume, we like 
to express the term as the integral of the volume density ^ with respect 
to volume: 



Obtaining a value for this integral may be difficult, especially if the volume 
of the control volume CV is itself a function of time. 

The term Bnet can be evaluated by quantifying, for each area element dA of 
the control volume’s surface, the surface flow rate pbVj_ of property B that 
flows through it (fig. 3.2). The integral over the entire control volume surface 
CS of this term is: 


^net ~ 



L 


pb (Vrel • n) dA 


(3/3) 


where flows and velocities are positive outwards and negative inwards by convention, 
CV is the control volume, 

CS is the the control surface (enclosing the control volume), 
n is a unit vector on each surface element dA pointing outwards, 

Vrel is Jhe local velocity of fluid relative to the control surface, 
and V_L = Vrei • n is the local cross-surface speed. 


By inserting equations 3/2 and 3/3 into equation 3/1, we obtain: 




N 

U 



Figure 3.2 - Part of the system may be flowing through an arbitrary piece of the 
control surface with area dA. The n vector defines the orientation of dA surface, 
and by convention is always pointed outwards. 

Figure CC-0 o.c. 


63 









(3/4) 


Equation 3/4 is named the Reynolds’ transport theorem-, it stands now as a 
general, abstract accounting tool, but as we soon replace B by meaningful 
variables, it will prove extremely useful, allowing us to quantify the net effect 
of the flow of a system through a volume for which border properties are 
known. 

In the following sections we are going to use this equation to assert four key 
physical principles (§0.7) in order to analyze the flow of fluids: 

• mass conservation; 

• change of linear momentum; 

• change of angular momentum; 

• energy conservation. 


3.3 Mass conservation 



Video: with sufficient skills (and 
lots of practice!), it is possible 
for a musician to produce an un¬ 
interrupted stream of air into 
an instrument while still con¬ 
tinuing to breathe, a technique 
called circular breathing. Can 
you identify the different terms 
of eq. 3/5 as they apply to the 
clarinetist’s mouth? 

by David Hernando Vitores (CC-by-sa) 
https://frama.link/yVesHaSk 


In this section, we focus on simply asserting that mass is conserved (eq. 0/22 
p.l9). Our study of the fluid’s properties at the borders of the control volume 
is made by replacing variable B by mass m. Thus becomes which 
by definition is zero. 

In a similar fashion, b = B/m = mim = 1 and consequently the Reynolds 
transport theorem (3/4) becomes: 


d/TZgyg 

dt 


^ Iff + JJ^ P (hel • n) dA (3/5) 


the time change 
of the system’s mass 


the rate of change 
of mass inside + 
the control volume 


the net mass flow 
at the borders 
of the control volume 


This equation 3/5 is often called continuity equation. It allows us to compare 
the incoming and outgoing mass flows through the borders of the control 
volume. 

When the control volume has well-defined inlets and outlets through which 
the term p(Vrel • n) can be considered uniform (as for example in fig. 3.3), this 
equation reduces to: 

0=+Z +Z 


(3/6) 


(3/7) 
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In equation 3/6, the term pVj^A at each inlet or outlet corresponds to the local 
mass flow ±m (positive inwards, negative outwards) through the boundary. 







Kel 


I A2I 
I I 



Figure 3.3 - A control volume for which the system’s properties are uniform at 
each inlet and outlet. Here eq. 3/5 translates as 0 = -4 fff pd'V + p 3 |V'j_ 3 |A 3 + 
P2lV^2lA2 - 

Figure CC-0 o.c. 


With equation 3/7 we can see that when the flow is steady (§0.8), the last two 
terms amount to zero, and the integral p dV (the total amount of mass 
in the control volume) does not change with time. 


3.4 Change of linear momentum 


In this section, we apply Newton’s second law: we assert that the variation 
of system’s linear momentum is equal to the net force being applied to it 
(eq. 0/23 p.l9). Our study of the fluid’s properties at the control surface is 
carried out by replacing variable B by the quantity mV: momentum. Thus, 

becomes which is equal to Fnet, the vector sum of forces applied 

on the system as it transits the control volume. 

In a similar fashion, b = B/m = V and equation 3/4, the Reynolds transport 
theorem, becomes: 



net 


the vector sum of 
forces on the system 


iffL 


pvdrv + 


If. 


pV (Vrel • n) dA 
(3/8) 


the rate of change 
of linear momentum + 

within the control volume 


the net flow of linear momen¬ 
tum through the boundaries 
of the control volume 


When the control volume has well-defined inlets and outlets through which 

—^ ^ 

the term pV (t^el • can be considered uniform (fig. 3.3), this equation reduces 

to: 


^net 


iffl pRdV + g {{p\VJA)v] - ^ {{p\V^\A)v} (3/9) 
fffcv ^ ~ 
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I A2I 
I I 



Figure 3.4 - The same control volume as in fig. 3.3. Here, since the system’s proper¬ 
ties are uniform at each inlet and outlet, eq. 3/8 translates as Fnet = JJfcv 

+ ^21^X21^2^2 “ Pl|Vj_i|AiVi. 

Figure CC-0 o.c. 


Let us observe the four terms of equation 3/10 for a moment, for they are full 
of subtleties. 

First, we notice that even if the flow is steady (an therefore that Xnet (^) = 
0 kg s“^ since = 0), the last two terms do not necessarily cancel each other 

(i.e. it is possible that Xnet 0). 

The reverse also applies: it is quite possible that Fnet ^ 0 even if the net 
momentum flow through the boundaries is null (that is, even if ^^net {jnV^ = 

0). This would be the case, if fJJ^y p dV (the total amount of momentum 
within the borders of the control volume) varies with time. Walking forwards 
and backwards within a rowboat would cause such an effect. 



Video: as a person walks, the 
deflection of the air passing 
around their body can be used 
to sustain the flight of a paper 
airplane (a walkalong glider). 
Can you figure out the momen¬ 
tum flow entering and leaving a 
control volume surrounding the 
glider? 

by Y:sciencetoymaker (styl) 
https://youtu.be/S6JKwzK37_8 


From equation 3/10 therefore, we read that two distinct phenomena can 
result in a net force on the system: 

—> 

• A difference between the values of mV of the fluid at the entrance and 
exit of the control volume (caused, for example, by a deviation of the 
flow or by a mass flow imbalance); 

• A change in time of the momentum mV within the control volume 
(for example, with the acceleration or the variation of the mass of the 
control volume). 

These two factors may cancel each other, so that the system may well be able 
to travel through the control volume without any net force being applied 
to it. 
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3.5 Change of angular momentum 


In this third spin on the Reynolds transport theorem, we assert that the 
change of the angular momentum of a system about a point X is equal to 
the net moment applied on the system about this point (eq. 0/24 p.l9). Our 
study of the fluid’s properties at the borders of the control volume is made 
by replacing the variable B by the angular momentum fxm A mV. Thus, 

becomes , which is equal to Mnet) the vector sum of moments 

applied on the system about point X as it transits through the control volume. 

In a similar fashion, b = B/m = rAV and equation 3/4, the Reynolds transport 
theorem, becomes: 

Mnet,X = JfJ A pV d^ + JJ Txm ^ P (t^rel ' n)V dA 

(3/11) 

the sum of the rate of change of the net flow of angular 

moments applied = the angular momentum + momentum through the 
to the system in the control volume control volume’s boundaries 


d(^Xm ^ ^^)sys 
dt 



Video: pre-lecture briefing for 
this chapter, part 2/2 

by o.c. (CC-by) 
https://youtu.be/nmEe7Dq0lAU 


in which Am is a vector giving the position of any mass m relative to point X. 


When the control volume has well-defined inlets and outlets through which 
the term fxm A p (Vrel • n)V can be considered uniform (fig. 3.5), this equation 
reduces to: 


MnetX = JJj 0(m A pr dV + ^ |o<m A inV) - ^ A mV} 



□ 



Video: rocket landing gone 

wrong. Can you compute the 
moment exerted by the top 
thruster around the base of the 
rocket as it (unsuccessfully) at¬ 
tempts to compensate for the 
collapsed landing leg? 

by Y:SciNews (styl) 
https://youtu.be/4cvGGxTsQxO 


Figure 3.5 - A control volume for which the properties of the system are uniform 
at each inlet or outlet. Here the moment about point X is Mnet,x ~ Jlfcv ^ 
pV &V + r2 A \m 2 \V 2 - Ti A 

Figure CC-0 o.c. 


Equation 3/12 allows us to quantify, with relative ease, the moment exerted 
on a system based on inlet and outlet velocities of a control volume. 
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3.6 Energy conservation 


We conclude our frantic exploration of control volume analysis with the first 
principle of thermodynamics. We now simply assert that the change in the 
energy of a system can only be due to well-identified transfers (eq. 0/25 p.l9). 
Our study of the fluid’s properties at the borders of the control volume is 
made by replacing variable B by an amount of energy Egys- Now, can be 
attributed to three contributors: 


dE. 


'sys 


_ 1 . 1 _ —I— _ Cl 


(3/13) 


where Qnet in is the net power transferee! as heat; 

!^shaft, net in IS the net power added as work with a shaft; 
and Vl^ressure, net in IS the net power required to enter and leave the control volume. 


It follows that b = B/m = E/m = e; and e is broken down into 


6 — i 6p (3/14) 

where i is the specific internal energy (Jkg“^); 

ej^ the specific kinetic energy (Jkg“^); 
and ep the specific potential energy (Jkg“^). 

Now, the Reynolds transport theorem (equation 3/4) becomes: 


d£, 


sys 


dt 


— Qnet 


in “f ^^haft, net in “f ^^ressure, net in 


d 

dt 


fff ped^-h rr pe(Vrern)dA 

JJJcv JJes 


(3/15) 


When the control volume has well-defined inlets and outlets through which 
the term p e • fi) can be considered uniform, this equation reduces to: 


Qnet in “f l^haft, net in “f l^ressure, net in 


Qnetin + l^haft, net ir 


- Ej + ^k + ep ) I i^ressure, net in 

in 


Making use of the concept of enthaipy defined as h = i + pip, we obtain: 


Cnet in + kl^haft, net in = JJJ^ p e dV -t ^ |m |/l -t -t pz j| - ^ |m |/l -t 1pz j| 


(3/17) 


the net power 
received by the system 
68 


the rate of change 
of energy inside -h 
the control volume 


the net energy flow rate 
exiting the control volume boundaries 






This equation 3/17 is particularly attractive, but it necessitates the input of 
a large amount of experimental data to provide useful results. It is indeed 
very difficult to predict how the terms i and pip will change for a given flow 
process. For example, a pump with given power Qnet in and Wshaft, net in will 
generate large increases of terms p, V and z if it is efficient, or a large increase 
of terms i and 1/p if it is inefficient. This equation 3/17, sadly, does not allow 
us to quantify the net effect of shear and the extent of irreversibilities in a 
fluid flow. 


3.7 The Bernoulli equation 


The Bernoulli equation has very little practical use for us; nevertheless it is 
so widely used that we have to dedicate a brief section to examining it. We 
will start from equation 3/17 and add five constraints: 

1. Steady flow. 

Thus^]j(jrc^pedW = 0. 

In addition, m has the same value at inlet and outlet; 

2. Incompressible flow. 

Thus, p stays constant; 

3. No heat or work transfer. 

Thus, both Qnet in and Wshaft, net in are zero; 

4. No friction. 

Thus, the fluid energy i cannot increase due to an input from the 
control volume; 

5. One-dimensional flow. 

Thus, our control volume has only one known entry and one known 
exit, all fluid particles move together with the same transit time, and 
the overall trajectory is already known. 


With these five restrictions, equation 3/17 simply becomes: 

0 = ^ |m(i + ^ + Ck + ep)| - ^ ■|m(i -l- ^ -l- -t ep)| 


out 

rh 


(j + — + -V 2 + gzz) - (i + — + -V^ + gzi) 
p 2, P ^ 


and we here obtain the Bernoulli equation: 


^ + ^V,^+gz,=^ + ^Vi+gz2 

p 2 P 2 


(3/18) 


This equation describes the properties of a fluid particle in a steady, incom¬ 
pressible, friction-less flow with no energy transfer. 

The Bernoulli equation can also be obtained starting from the linear momen¬ 
tum equation (eq. 3/8 p.65): 

F„et = pV dW + pV (Vrel ’ «) dA 
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When considering a fixed, infinitely short control volume along a known 
streamline s of the flow, this equation becomes: 

fi-Fpressure + dfshear + dFgj-avity ~ IfL pVd^ + pVAdV 

along a streamline, where the velocity V is aligned (by definition) with the streamline. 

Now, adding the restrictions of steady flow (d/ dt = 0) and no friction 
(dFshear = 0)> we already obtain: 

dFpi-essure d" dFgi-avity ~ pVA dV^ 

—^ 

The projection of the net force due to gravity dfgravity on the streamline 
segment ds has norm dFgravity • d? = -gpA dz, while the net force due to 
pressure is aligned with the streamline and has norm dfpressure,s = -A dp. 
Along this streamline, we thus have the following scalar equation, which we 
integrate from points 1 to 2: 


-A dp - pgA dz = pVA dV 
--dp-gdz = VdV 

p r ^ 




The last obstacle is removed when we consider flows without heat or work 
transfer, where, therefore, the density p is constant. In this way, we arrive to 
equation. 3/18 again: 


h + lv,^ + jz, = ^ + ly? + 9Z2 

p 2 P 2 


Let us insist on the incredibly frustrating restrictions brought by the five 
conditions above: 

1. Steady flow. 

This constrains us to continuous flows with no transition effects, which 
is a reasonable limit; 

2. Incompressible flow. 

We cannot use this equation to describe flow in compressors, turbines, 
diffusers, nozzles, nor in flows where M > 0,6. 

3. No heat or work transfer. 

We cannot use this equation in a machine (e.g. in pumps, turbines, 
combustion chambers, coolers). 

4. No friction. 

This is a tragic restriction! We cannot use this equation to describe a 
turbulent or viscous flow, e.g. near a wall or in a wake. 

5. One-dimensional flow. 

This equation is only valid if we know precisely the trajectory of the 
fluid whose properties are being calculated. 

Among these, the last is the most severe (and the most often forgotten): 
the Bernoulli equation does not allow us to predict the trajectory of 
fluid particles. Just like all of the other equations in this chapter, it requires 
a control volume with a known inlet and a known outlet. 
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3.8 Limits of integral analysis 


Integral analysis is an incredibly useful tool in fluid dynamics: in any given 
problem, it allows us to rapidly describe and calculate the main fluid phenom¬ 
ena at hand. The net force exerted on the fluid as it is deflected downwards 
by a helicopter, for example, can be calculated using just a loosely-drawn 
control volume and a single vector equation. 

As we progress through exercise sheet 3, however, the limits of this method 
slowly become apparent. They are twofold. 

• First, we are confined to calculating the net effect of fluid flow. The net 
force, for example, encompasses the integral effect of all forces —due to 
pressure, shear, and gravity— applied on the fluid as it transits through 
the control volume. Integral analysis gives us absolutely no way of 
distinguishing between those sub-components. In order to do that (for 
example, to calculate which part of a pump’s mechanical power is lost 
to internal viscous effects), we would need to look within the control 
volume. 

• Second, all four of our equations in this chapter only work in one 
direction. The value dflsys/ dt of any finite integral cannot be used to 
find which function pbV_i dA was integrated over the control surface to 
obtain it. For example, there are an infinite number of velocity profiles 
which will result in a net force of -12 N. Knowing the net value of an 
integral, we cannot deduce the conditions which lead to it. 

In practice, this is a major limitation on the use of integral analysis, 
because it confines us to working with large swaths of experimental 
data gathered at the borders of our control volumes. From the wake 
below the helicopter, we deduce the net force; but the net force tells us 
nothing about the shape of the wake. 

Clearly, in order to overcome these limitations, we are going to need to open 
up the control volume, and look at the details of the flow within — perhaps 
by dividing it into a myriad of sub-control volumes. This is what we set 
ourselves to in chapter 4, with a thundering and formidable methodology we 
shall call derivative analysis. 
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Exercise sheet 3 (Analysis of existing flows) 

last edited September 4, 2017 


These lecture notes are based on textbooks by White [13], Qengel & al.[16], and Munson & al.[18]. 


Except otherwise indicated, we assume that fluids are Newtonian, and that: 

Pwater = 1 000 kg m“^ patm. = 1 bar; Patm. = 1,225 kgm“^ Tatm. = ll,3°C;patm. = 1,5 • 10‘^Nsm“2; 
g = 9,81 ms“^. Air is modeled as aperfect gas (i?air = 287JK“^ kg“'; pair = 1,4; Cpair = 1005Jkg“^ K“^). 


Reynolds Transport Theorem: 

dBsys d 

dt dt 


fff pbd'y+ ff pb(Vrei-n)dA (3/4) 

JJJcv JJcs 


Mass conservation; 

dm 


sys ^ u 

—^ = 0 = — 

dt dt 


fff pdrV+ ff p(Vrel-n)dA (3/5) 

JJJcv JJcs 


Change in linear momentum: 
d(mVsys) 


dt 


= Fr 


net 




Change in angular momentum: 


d(rXm ^ ^^)sys d 

- 1 , -= ""««=d7 


pVd’V+ It pV(V„i-ii)iA (3/8) 
lev JJcs 


fff rxmApVd^+ ff rxmAp(Vrei-n)VdA 

JJJcv JJcs 

(3/11) 


3.1 Pipe bend 

^ CC-0 o.c. 

A pipe with diameter 30 mm has a bend with angle 0 = 130°, as shown in fig. 3.6. Wa¬ 
ter enters and leaves the pipe with the same speed Vi = R 2 = 1,5 ms“^ The velocity 
distribution at both inlet and outlet is uniform. 



Figure 3.6 - A pipe bend, through which water is flowing. We assume that the velocity distribution 
at inlet and outlet is identical. 

Figure CC-0 o.c. 
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1. What is the mass flow traveling through the pipe? 

2. What is the force exerted by the pipe bend on the water? 

3. Represent the force vector qualitatively. 

4. What would be the new force if all of the speeds were doubled? 


3.2 Water jet 


White [13] Ex3.9 


A horizontal water jet hits a vertical wall and is split in two symmetrical vertical flows 
(fig. 3.7). The water nozzle has a 3 cm^ cross-sectional area, and the water speed at the 
nozzle outlet is = 20ms“h 

In this exercise, we consider the simplest possible case, with the effects of gravity and 
viscosity neglected. In that case, the water is perfectly split into two vertical jets of 
identical mass flow, the speed is identical everywhere, and the velocity distributions are 
uniform. 


1. What is the velocity of each water flow as it leaves the wall? 

2. What is the force exerted by the water on the wall? 

Now, the wall moves longitudinally in the same direction as the water jet, with a speed 

Vwall = 15ms“h 

3. What is the new force exerted by the water on the wall? 

4. What is the velocity of each water jet as it leaves the wall, relative to the nozzle? 

5. How would the force be modified if the volume flow was kept constant, but the 

diameter of the nozzle was reduced? (briefly justify your answer, e.g. in 30 words 
or less) 
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Figure 3.7 - A water jet flowing out of a nozzle (left), and impacting a vertical wall on the right. 

Figure CC-0 o.c. 












3.3 Exhaust gas deflector 


A deflector is used behind a stationary aircraft during ground testing of a jet engine 
(fig. 3.8). 



Figure 3.8 - A mobile exhaust gas deflector, used to deflect hot jet engine exhaust gases upwards 
during ground tests. 

Figure CC-0 o.c. 


The deflector is fed with a horizontal air jet with a quasi-uniform velocity profile; the 
speed is fjet = 600 kmh“\ temperature 400 °C and the pressure is atmospheric. As the 
exhaust gases travel through the pipe, their heat losses are negligible. Gases are rejected 
with a 40° angle relative to the horizontal. 

The inlet diameter is 1 m and the horizontal outlet surface is 6 m^. 

1. What is the force exerted on the ground by the deflection of the exhaust gases? 

2. Describe qualitatively a modification to the deflector that would reduce the hori¬ 
zontal component of force. 

3. What would the force be if the deflector traveled rearwards (positive x-direction) 
with a velocity of 10 m s“^? 


3.4 Felton water turbine 

White [13] P3.56 

A water turbine is modeled as the following system: a water jet exiting a stationary 
nozzle hits a blade which is mounted on a rotor (fig. 3.9). In the ideal case, viscous effects 
can be neglected, and the water jet is deflected entirely with a 180° angle. 



Figure 3.9 - Schematic drawing of a water turbine blade. This type of turbine is called Felton 
turbine. 

Figure CC-0 o.c. 
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The nozzle has a cross-section diameter of 5 cm and produces a water jet with a speed 
t^et = 15 m s“^ The rotor diameter is 2 m and the blade height is negligibly small. 

We first study the case in which the rotor is prevented from rotating, so that the blade is 
stationary (Vjjiade = Oms“^). 

1. What is the force exerted by the water on the blade? 

2. What is the moment exerted by the blade around the rotor axis? 

3. What is the power transmitted to the rotor? 

We now let the rotor rotate freely. Friction losses are negligible, and it accelerates until it 
reaches maximum velocity. 

4. What is the rotor rotation speed? 

5. What is the power transmitted to the rotor? 

The rotor is now coupled to an electrical generator. 

6. Show that the maximum power that can be transmitted to the generator occurs for 
Mjlade ~ jl^ater- 

7. What is the maximum power that can be transmitted to the generator? 

8. How would the above result change if viscous effects were taken into account? 
(briefly justify your answer, e.g. in 30 words or less) 


3.5 Snow plow 


derived from Gerhart & Gross [6] Ex5.9 


A road-based snow plow (fig. 3.10) is clearing up the snow on a flat surface. We wish to 
quantify the power required for its operation. 

The snow plow is advancing at 25 kmh“^; its blade has a frontal-view width of 4 m. 


The snow on the ground is 30 cm deep and has density 300 kg m 


The snow is pushed along the blade and is rejected horizontally with a 30° angle to the 
left of the plow. Its density has then risen to 450kgm“^. The cross-section area Aoutlet of 
the outflowing snow in the x-y plane is 1,1 m^. 


1. What is the force exerted on the blade by the deflection of the snow? 

(Indicate its magnitude and coordinates) 

2. What is the power required for the operation of the snow plow? 

3. If the plow velocity was increased by 10 %, what would be the increase in power? 
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25 km/h 





Figure 3.10 - Outline schematic of a blade snow plow. 

Figure CC-0 o.c. 

3.6 Drag on a cylindrical profile 

In order to measure the drag on a cylindrical profile, a cylindrical tube is positioned 
perpendicular to the air flow in a wind tunnel (fig. 3.11), and the longitudinal component 
of velocity is measured across the tunnel section. 

Upstream of the cylinder, the air flow velocity is uniform (ui = 17 = 30 m s“^). 

Downstream of the cylinder, the speed is measured across a 2 m height interval. Hori¬ 
zontal speed measurements are gathered and modeled with the following relationship: 

U2iy) = 29 + (3/19) 

The width of the cylinder (perpendicular to the flow) is 2 m. The Mach number is very 
low, and the air density remains constant at p = 1,23 kg m“^; pressure is uniform all 
along the measurement field. 

• What is the drag force applying on the cylinder? 

• How would this value change if the flow in the cylinder wake was turbulent, 
and the function U 2 {y) above only modeled time-averaged values of the horizontal 
velocity? (briefly justify your answer, e.g. in 30 words or less) 
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Figure 3.11 - A cylinder profile set up in a wind tunnel, with the air flowing from left to right. 

Figure CC-0 o.c. 


3.7 Drag on a flat plate 

We wish to measure the drag applying on a thin plate positioned parallel to an air stream. 
In order to achieve this, measurements of the horizontal velocity u are made around the 
plate (fig. 3.12). 


Ujly) - U U2(yj 



Figure 3.12 - Side view of a plate positioned parallel to the flow. 


Figure CC-0 o.c. 


At the leading edge of the plate, the horizontal velocity of the air is uniform: Ui = U = 
10ms“^. 

At the trailing edge of the plate, we observe that a thin layer of air has been slowed down 
by the effect of shear. This layer, called boundary layer, has a thickness of 5 = 1 cm. The 
horizontal velocity profile can be modeled with the relation: 

U2{y) = u{^^y (3/20) 

The width of the plate (perpendicular to the flow) is 30 cm and it has negligible thickness. 
The flow is incompressible (p = 1,23 kg m“^) and the pressure is uniform. 
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• What is the drag force applying on the plate? 

• What is the power required to compensate the drag? 



























































































Under which form is the kinetic energy lost by the flow carried away? Can this 
new form of energy be measured? (briefly justify your answer, e.g. in 30 words or 
less) 


3.8 Drag measurements in a wind tunnel 


A group of students proceeds with speed measurements in a wind tunnel. The objective 
is to measure the drag applying on a wing profile positioned across the tunnel test section 
(fig. 3.13). 



Figure 3.13 - Wing profile positioned across a wind tunnel. The horizontal velocity distributions 
upstream and downstream of the profile are also shown. 

Figure CC-0 o.c. 


Upstream of the profile, the air flow velocity is uniform (ui = U = 50 m s ^). 

Downstream of the profile, horizontal velocity measurements are made every 5 cm across 
the flow; the following results are obtained: 


vertical position (cm) 

horizontal speed U 2 (m s ^) 

0 

50 

5 

50 

10 

49 

15 

48 

20 

45 

25 

41 

30 

39 

35 

40 

40 

43 

45 

47 

50 

48 

55 

50 

60 

50 


The width of the profile (perpendicular to the flow) is 50 cm. The airflow is incompressible 
{p = 1,23 kg m“^) and the pressure is uniform across the measurement surface. 

1. What is the drag applying on the profile? 

2. How would the above calculation change if vertical speed measurements were also 
taken into account? 
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3.9 Moment on gas deflector 

We revisit the exhaust gas deflector of exercise 3.3 p.75. Figure 3.14 below shows the 
deflector viewed from the side. The midpoint of the inlet is 2 m above and 5 m behind the 
wheel labeled “A”, while the midpoint of the outlet is 3,5 m above and 1,5 m behind it. 



3,5 

2 


Figure 3.14 - Side view of the mobile exhaust gas deflector which was shown in fig. 3.8 p.75 

Figure CC-0 o.c. 


What is the moment generated by the gas flow about the axis of the wheel labeled “A”? 


3.10 Helicopter tail moment 

In a helicopter, the role of the tail is to counter exactly the moment exerted by the main 
rotor about the main rotor axis. This is usually done using a tail rotor which is rotating 
around a horizontal axis. 

A helicopter, shown in flgure 3.15, is designed to use a tail without a rotor, so as to reduce 
risks of accidents when landing and taking off. The tail is a long hollow cylindrical tube 
with two inlets and one outlet. To simplify calculations, we consider that pressure is 
atmospheric at every inlet and outlet. 

• inlet A has a cross-section area of 0,2 m^. It contributes hot exhaust gases of density 
0,8 kg m“^ and velocity 12 m s“\ aligned with the (x) axis of the tail; 

• inlet B contributes 25kgs“^ of atmospheric air incoming at an angle a = 130° 
relative to the axis of the tail, with a velocity of 3 m s“k 

The mix of exhaust gases and atmospheric air is rejected at the tip of the tail (outlet C) 
with a flxed velocity of 45 m s“k The angle 9 at which gases are rejected is controlled by 
the flight computer. 

1. What is the rejection angle 9 required so that the tail generates a moment of -1-6 kN m 
around the main rotor (y) axis? 

2. Propose and quantify a modification to the tail geometry or operating conditions 
that would allow the tail to produce no thrust (that is to say, zero force in the 
x-axis), while still generating the same moment. 























Figure 3.15 - Top-view of a helicopter using a rotor-less tail. 

Figure derived from a figure CC-by-sa by Commons User.FOX 52 


Remark: this system is commercialized by MD Helicopters as the notar. The use of exhaust 
gases was abandoned, however, a clever use of air circulation around the tail pipe axis 
contributes to the generated moment; this effect is explored in chapter 8 (§8.3.4 p. 176). 


3.11 Pressure losses in a pipe flow ,, , , , , 

^ Non-examinable. Based on White [13] 

Water is circulated inside a cylindrical pipe with diameter 1 m (fig. 3.16). 


^2(r) 



Figure 3.16 - Velocity profiles at the inlet and outlet of a circular pipe. 


Figure CC-0 o.c. 


At the entrance of the pipe, the speed is uniform: Ui = Uav. = Sms h 

At the outlet of the pipe, the velocity profile is not uniform. It can be modeled as a 
function of the radius with the relationship: 


U2{r) - 


U, 


center 



(3/21) 
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The momentum flow at the exit may be described solely in terms of the average velocity 
f/av. with the help of a momentum flux correction factor such that p f u^dA = jfpAU^^ . 

When U2 = f/center (l - A Can be shown that ^ m= j, 

we have yS « 1,02. 

[difficult question] 

What is the net force applied on the water so that it may travel through the pipe? 


3.12 Thrust reverser 

non-examinable 

A turbofan installed on a civilian aircraft is equipped with a thrust inverter system. When 
the inverters are deployed, the “cold” outer flow of the engine is deflected outside of the 
engine nacelle. 


We accept that the following characteristics of the engine, when it is running at full 
power, and measured from a reference point moving with the engine, are independent of 
the aircraft speed and engine configuration: 


Inlet diameter 
Inlet speed 
Inlet temperature 
Cold flow outlet velocity 
Cold flow temperature 
Hot flow outlet diameter 
Hot flow outlet velocity 
Hot flow outlet temperature 
Air pressure at all exits 


D\ = 2,5 m 
Vi = 60ms“^ 

Ti = 288 K 
Vj, = 85ms“^ 

To, = 300 K 
D 4 = 0,5 m 
V4 = 200 ms“^ 

T 4 = 400 K 

p3 = pi = patm. = 1 bar 



Figure 3.17 - Conceptual layout sketch of a turbofan. The airflow is from left to right. 

Figure CC-by o.c. 


We first study the engine in normal (forward thrust) configuration. 

1. When mounted on a stationary test bench, what is the thrust provided by the 
engine? 

2. In the approximation that the operating parameters remain unchanged, what is 
the thrust generated when the engine moves with a speed of 40 m s“^ ? 

We now deploy the thrust reversers. They deflect the cold flow with a 75° angle relative 
to the aircraft longitudinal axis. 

3. What is the net thrust developed when the aircraft velocity is 40 m s“^? 

4. What is the net thrust developed when the aircraft is stationary? 
























Answers 

3.1 1) m = 1,0603 kg 2) Fnet = (-0,5681;-1,2184) N 

3) The force is quadrupled. 

3.2 1) iViet sys. — —120 N; 2) iViet sys. — —7,5 N, V^xit, absolute — (15; 5); 

fexit, absolute = 15,8niS 

3.3 1) Fnet X = -9,532 kN & Fnet y = +1,479 kN : fnet = 9,646 kN; 

2) Tnet2 = 8,525 kN. 

3.4 1) Fnet = -883,6 N; 2) Mnetx = |Tnet|f? = 883,6 N m; 

3) Wrotor = 0 W; 4) CO = 143,2 rpm (Fnet = 0 N); 

5) Wrotor = 0 W again; 

0) li^otor, max — 1,963 kW^ @ Vjbiade, optimal — ^kwaterjet- 

3.5 1) m = 2500kgs“^; Fnet X — +10,07 kN, Fnet z — +12,63 kN (force on blade is oppo¬ 
site); 2) W = FnefTplow = -FnetxlVi| = 69,94kW 

3.6 Fnetx = pLjg2 (“2 “ ^“ 2 ) dy = -95,78 N. 

3.7 Fnetx = pL ff (u^y^ - UU(y)) dy = -7,175 • lO-^ N : Wdrag = f7|J"netxl = 0,718 W. 

3.8 Fnetx » pL^y [{uj - UU 2 ) (5l/j = -64,8 N. 

3.9 Re-use rh = 67,76 kg s“\ Vi = 166,7 m s“\ Vz = 33,94 m s“^ from ex.3.3. From trigonom¬ 
etry, find ^ 22.^2 - 1,717m. Then, plug in eq. 3/12 p.67: Mnet = -18,64kNm in 
z-direction. 

3.10 1) Work eq. 3/11 down to scalar equation (in y-direction), solve for 9: 9 = 123,1°. 

2) There are multiple solutions which allow both moment and force equations to 
be solved at the same time, rc can be shortened, the flow in C can be split into 
forward and rearward components, or tilted downwards etc. Reductions in or 
Vc are also possible, but quantifying them requires solving both equations at once. 

3.11 Reenter = 1,224517; Fnet = +393 N (positive!) 

3.12 /hcold — 297 kgs , rilhot — 59,4kgs ; Feold flow, normal, bench & runway — +74,25kN, 

■Fhot flow, normal & reverse, bench & runway — +8,316 kN; 

■Fcold flow, reverse, bench & runway — —24,35 kN and F^ot flow, normal, bench & runway — +8,316 kN. 
Adding the net pressure force due to the (lossless) flow acceleration upstream of the 
inlet, we obtain, on the bench: Fenglne bench, normal — —93,26 kN, Fenglne bench, reverse — 
+5,344 kN; and on the runway: Fgnglne runway, normal — —92,5 kN and Fgnglne bench, reverse — 
+6,094 kN. 
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These lecture notes are based on textbooks by White [13], Qengel & al.[16], and Munson & al.[18]. 


4.1 Motivation 


In this chapter we assign ourselves the daunting task of describing the 
movement of fluids in an extensive manner. We wish to express formally, 
and calculate, the dynamics of fluids —the velocity field as a function of 
time— in any arbitrary situation. For this, we develop a methodology named 
derivateive analysis. 

Let us not shy away from the truth: the methods developed here are dis¬ 
proportionately complex in comparison to the solutions they allow us to 
derive by hand. Despite this, this chapter is extraordinarily important, for 
two reasons: 



Video: pre-lecture briefing for 
this chapter 

by o.c. (CC-by) 
https://youtu.be/wOGlqykD3_E 


• derivative analysis allows us to formally describe and relate the key 
parameters that regulate fluid flow, and so is the key to developing an 
understanding of any fluid phenomenon, even when solutions cannot 
be derived; 

• it is the backbone for computational fluid dynamics (cfd) in which 
approximate flow solutions are obtained using numerical procedures. 


4.2 Eulerian description of fluid flow 


4.2.1 Problem description 

From now on, we wish to describe the velocity and pressure fields of a fluid 
with the highest possible resolution. For this, we aim to predict and describe 
the trajectory of fluid particles (§0.2.2) as they travel. 
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Newton’s second law allows us to quantify how the velocity vector of a 
particle varies with time. If we know all of the forces to which one particle 
is subjected, we can obtain a streak of velocity vectors l^artide = {u, v, w) = 
f{x,y,z,t) as the particle moves through our area of interest. This is a 
description of the velocity of one particle; we obtain a function of time 
which depends on where and when the particle started its travel: I^artide = 
f{xo,yo,zoJoJ)- This is the process used in solid mechanics when we wish 
to describe the movement of an object, for example, a single satellite in orbit. 

This kind of description, however, is poorly suited to the description of fluid 
flow, for three reasons: 

• Firstly, in order to describe a given fluid flow (e.g. air flow around the 
rear-view mirror standing out of a car), we would need a large number 
of initial points {xo,yo,zoJo), and we would then obtain as many trajec¬ 
tories V^artide- It then becomes very difficult to describe to study and 
describe a problem that is local in space (e.g. the wake immediately 
behind the car mirror), because this requires finding out where the 
particles of interest originated, and accounting for the trajectories of 
each of them. 

• Secondly, the concept of a “fluid particle” is not well-suited to the draw¬ 
ing of trajectories. Indeed, not only can particles strain indefinitely, 
but they can also diffuse into the surrounding particles, “blurring” and 
blending themselves one into another. 

• Finally, the velocity and other properties of a given particle very 
strongly depend on the properties of the surrounding particles. We 
have to resolve simultaneously the movement equations of all of the 
particles. A space-based description of properties —one in which we 
describe properties at a chosen fixed point of coordinates Xpoint.i/point. 
Zpoint> fpoint~ is much more useful than a particle-based description 
which depends on departure points xo,yo,zo,to. It is easier to deter¬ 
mine the acceleration of a particle together with that of its current 
neighbors, than together with that of its initial (former) neighbors. 

What we are looking for, therefore, is a description of the velocity fields 
that is expressed in terms of a fixed observation point I^oint = {u,v,w) = 
f{xpoinuypoinuZpoinut), through which particles of many different origins 
may be passing. This is termed a Eulerian flow description, as opposed to the 
particle-based Lagrangian description. Grouping all of the point velocities in 
our flow study zone, we will obtain a velocity field I^oint that is a function of 
time. 

4.2.2 Total time derivative 

Let us imagine a canal in which water is flowing at constant and uniform 
speed u = I/canal (Ag- 4.1). The temperature T^ater of the water is constant 
(in time), but not uniform (in space). We measure this temperature with a 
stationary probe, reading Tprobe = Twater on the instrument. Even though the 
temperature Twater is constant, when reading the value measured at the probe, 
temperature will be changing with time: 
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In the case where the water temperature is not simply heterogeneous, but 
also unsteady (varying in time with a rate dT/dt), then that local change in 
time will also affect the rate which is measured at the probe: 


dTw 


dt 

Clprobe 


dt 


0 

dT^ 


dt 


Iprobe 


^^water ^. 

dx ^water 


(4/1) 


We keep in mind that the rate dT jdt can itself be a function of time and 
space; in equation 4/1, it is its value at the position of the probe and at the 
time of measurement which is taken into account. 

This line of thought can be generalized for three dimensions and for any 
property A of the fluid (including vector properties). The property A of 
one individual particle can vary as it is moving, so that it has a distribution 
A = f{x, y, z, t) within the fluid. The time rate change of A measured at 
a point fixed in space is named the total time derivative or simply total 
derivative^of A and written DA/Dt: 

(4/2) 

(4/3) 
(4/4) 

(4/5) 




u u u u u 

- =-h U -h V -h W- 

Dt dt dx dy dz 

DA dA dA dA dA 

- = —-h U— -h V— -h W-— 

Dt dt dx dy dz 

DA dA dA dA dA 

- = —-h U- -h V— -h W-— 

Dt dt dx dy dz 

I + \ 

dt dx dy dz 

dAy dAy dAy dAy 


dt 

dA, 

dt 


+ U 
+ U 


dx 

dx 


+ V 
+ V 


dA^ 

dy 


+ w 
+ w 


dz 
dA, 
dz / 


Figure 4.1 - A one-dimensional water flow, for example in a canal. The water has 
a non-uniform temperature, which, even if it is constant in time, translates in a 
temperature rate change in time at the probe. 

Figure CC-0 o.c. 


^Unfortunately this term has many denominations across the litetature, including ad- 
vective, convective, hydrodynamic, Lagrangian, particle, substantial, substantive, or Stokes 
derivative. In this document, the term total derivative is used. 


87 























In eq. equation 4/3, it is possible simplify the notation of the last three terms. 
We know the coordinates of the velocity vector V (by definition) and of the 
operator gradient V (as defined in eq. 1/8 p. 32) are: 


V = iu+jv + kw 
^ ^ d ^ d ^ d 

dx dy dz 


We thus define the advective operator, {V • V): 


^ d d d 

V • V = u- —I- V- —I- w— 
dx dy dz 


(4/6) 


We can now rewrite eqs. 4/2 and 4/3 in a more concise way: 


D 

Dt 


d ^ 

-.(y.v) 


(4/7) 


DA 


^4(F.V)A 


(4/8) 


As we explored in equation 4/1, the total derivative of a property can be 
non-null even if the flow is perfectly steady; similarly, the total derivative of 
a property can be zero (e.g. when a stationary probe provides a constant mea¬ 
surement) even if the properties of the particles have a non-zero time change 
(e.g. a fluid with falling temperature, but with a non-uniform temperature). 



Video: half-century-old, but 

timeless didactic exploration 
of the concept of Lagrangian 
and Eulerian derivatives, 
with accompanying notes by 
Lumley[3] 

by the National Committee for Fluid 
Mechanics Films (ncfmf, 1969[12]) (styl) 
https://youtu.be/mdN800kx2ko 


The total time derivative is the tool that we were looking for. Now, as we 
solve for the fluid properties in a stationary reference frame, instead of in the 
reference frame of a moving particle, we will replace all the time derivatives 
with total derivatives This allows us to compute the change in time 
of a property locally, without the need to track the movement of particles 
along our field of study. The first and most important such property we are 
instested in is velocity. Instead of solving for the acceleration of numerous 
particles (dV^article/ dt), we will focus on calculating the acceleration Geld 
DV/Dt. 


4.3 Mass conservation 


The first physical principle we wish to express in a derivative manner is the 
conservation of mass (eq. 0/22 p.l9). Let us consider a fluid particle of volume 
dW, at a given instant in time (fig. 4.2). We can reproduce our analysis from 
chapter 3 by quantifying the time change of mass within an arbitrary volume. 

In the present case, the control volume is stationary and the particle (our 
system) is flowing through it. We start with eq. 3/5 p.64: 
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Figure 4.2 - Conservation of mass within a fluid particle. In the x-direction, a mass 
flow rhi = JJP\U\ dz dy is flowing in, and a mass flow = ff P 2 U 2 dz dy is flowing 
out. These two flows may not be equal, since mass may also flow in the y- and 
z-directions. 

Figure CC-0 o.c. 


The first of these two integrals can be rewritten using the Leibniz integral 
rule: 



ffLi 

fflJ" 


(4/10) 


where Vs is the speed of the control volume wall; 

and where the termj[j^g pVg dA is simply zero because we chose a fixed control 
volume, such as a fixed computation grid I 


Now we turn to the second term of equation 4/9, p(i4el ■ which 

represents the mass flow thnet flowing through the control volume. 

In the direction x, the mass flow rhnet x flowing through our control volume 
can be expressed as: 


rhnetx= (I -pi|Mi|dzdt/+ (I p 2 \u 2 \dzdy 

JJcs JJcs 


=//J 

£ILi 


— (pu) dxdzdy 
ox 


(pu) dV 


(4/11) 


The same applies for directions y and z, so that we can write: 


IL 


p(Vrel ■n)dA= rhnet = rhnet X + rhnet y + fhnet : 


=IfL " s 


(pw) 


dV 

(4/12) 


Tn a CFD calculation in which the grid is deforming, this term pVs dA will have to 
be re-introduced in the continuity equation. 
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Now, with these two equations 4/10 and 4/12, we can come back to equa¬ 
tion 4/9, which becomes: 


df^particle ^ d 

dt dt 


flL 


pd^ + 


•= fffj" 'fIL 


J£^p(Vrd-n)dA 

V • (pV) d'V 


(4/13) 


Since we are only concerned with a very small volume d'V, we drop the 
integrals, obtaining: 


dp ^ ^ 

o = ^ + v.(py) 

(4/14) 

for all flows, with all fluids. 



This equation 4/14 is named continuity equation and is of crucial importance 
in fluid mechanics. All fluid flows, in all conditions and at all times, obey this 
law. 

Equation 4/14 is a scalar equation: it gives us no particular information about 
the orientation of V or about its change in time. In itself, this is insufficient 
to solve the majority of problems in fluid mechanics, and in practice it is used 
as a kinematic constraint to solutions used to evaluate their physicality or 
the quality of the approximations made to obtain them. 

The most assiduous readers will have no difficulty reading the following 
equation, 

^ + V • (pV) = ^ + V-Vp + p{V-V) = 0 (4/15) 

which allows us to re-write equation 4/14 like so: 

+ V -V = 0 (4/16) 

p Dt 


Therefore, we can see that for an incompressible flow, the divergent of 

—^ ^ 

velocity V • V (sometimes called volumetric dilatation rate) is zero: 




V-V = 0 

for any incompressible flow. 

(4/17) 


In three Cartesian coordinates, this is re-expressed as: 

du dv dw 
dx dy dz 


(4/18) 
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4.4 Change of linear momentum 


Now that we have seen what mass conservation tells us about the dynamics 
of fluids, we turn to the momentum equation, hoping to express the velocity 
field to the forces that are applied on the fluid. 


4.4.1 The Cauchy equation 


We start by writing Newton’s second law (eq. 0/23 p.l9) as it applies to a fluid 
particle of mass mpartidej as shown in fig. 4.3. Fundamentally, the forces on a 
fluid particle are of only three kinds, namely weight, pressure, and shear: ^ 


dV _ ^ 

^particle “ -^weight -fnet, pressure -fnet, sheai 


(4/19) 


We now apply this equation to a stationary cube of infinitesimal volume d'V 
which is traversed by fluid. Measuring the velocity from the reference frame 
of the cube, and thus using a Eulerian description, we obtain: 

^particle DE Ttl -> i i 0 

“ dT^ ^ pressure + -^Fnet, shear 

P ~ P9 “t pressure + ~^PyPnet, shear (4/20) 

Therefore, we will have obtained an expression for the acceleration field 
DE /Dt if we can find an expression for the net pressure and shear force. 
Let us quantify those efforts as they apply on each of the faces of the fluid 



Figure 4.3 - In our study of fluid mechanics, we consider only forces due to gravity, 
shear, or pressure. 

Figure CC-0 o.c. 


hn some very special and rare applications, electromagnetic forces may also apply. 
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particle. For pressure, we obtain: 


dy dz 


dx 


ax 


-dp 

dx 
-dp 

dy 
-dp 

-fnet, pressure,z ~ dH^ ~d^ 


Pnet, pressurepc 
Pnet, pressure,:/ 


d'V 

d^ 


(4/21) 

(4/22) 

(4/23) 


This pattern, we had seen in chapter 1 (with eq.1/8 p.32), can more ele¬ 
gantly be expressed with the gradient operator (see also Appendix A2 p.217), 
allowing us to write: 


1 


net, pressure 


= -Vp 


(4/24) 


As for shear, we can put our work from chapter 2 to good use now. There, 
we had seen that the effect of shear on a volume had eighteen components, 
the careful examination of which allowed us to quantify the net force due to 
shear (which has only three components). In the x-direction, for example, we 
had seen (with eq. 2/8 p.49, as illustrated again in fig. 4.4) that the net shear 
force in the x-direction Fsheaix was the result of shear in the x-direction on 
each of the six faces: 


fshear x — dx dy{T;^x 3 '^zx o) 

+ dx dz(ryx 2 ~ ^t/x 5) 
+ dz dy(Txx 1 - Txx 4) 


This, in turn, was more elegantly expressed using the divergent operator (see 
also Appendix A2 p.217) as equation 2/13: 



Figure 4.4 - Shear efforts on a cubic fluid particle (already represented in fig. 2.1 
p.48). 

Figure CC-0 o.c. 
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Finally, we had assembled the components in all three directions into a single 
expression for the net shear force with eq. 2/14 p.50: 

fshear = d^ V • Tij (4/27) 


Now, we can put together our findings, eq. 4/24 for the force due to pressure 
and eq. 4/27 for the force due to shear, back into equation 4/20, we obtain the 
Cauchy equation: 


II 

1 

+ 

<]i 

(4/28) 

for all flows, with all fluids. 



The Cauchy equation is an expression of Newton’s second law applied to a 
fluid particle. It expresses the time change of the velocity vector measured at 
a fixed point (the acceleration field DV/Dt) as a function of gravity, pressure 
and shear effects. This is quite a breakthrough. Within the chaos of an 
arbitrary flow, in which fluid particles are shoved, pressurized, squeezed, and 
distorted, we know precisely what we need to look for in order to quantify the 
time-change of velocity: gravity, the gradient of pressure, and the divergent 
of shear. 

Nevertheless, while it is an excellent start, this equation isn’t detailed enough 
for us. In our search for the velocity field V, the changes in time and space 
of the shear tensor tij and pressure p are unknowns. Ideally, those two 
terms should be expressed solely as a function of the flow’s other properties. 
Obtaining such an expression is what Claude-Louis Navier and Gabriel Stokes 
set themselves to in the 19* century: we follow their footsteps in the next 
paragraphs. 


4.4.2 The Navier-Stokes equation for incompressible flow 

The Navier-Stokes equation is the Cauchy equation (eq. 4/28) applied to 
Newtonian fluids. In Newtonian fluids, which we encountered in chapter 2 
(§2.3.3 p.52), shear efforts are simply proportional to the rate of strain; thus, 
the shear component of eq. 4/28 can be re-expressed usefully. 

We had seen with eq. 2/16 p.51 that the norm | | of shear component in 

direction j along a surface perpendicular to i can be expressed as: 

^ dui 


Thus, in a Newtonian fluid, on each of the six faces of the infinitesimal 
cube from fig. 4.4, the shear in x-direction is proportional to the partial 
derivative of u (the component of velocity in the x-direction) with respect to 
the direction perpendicular to each face: 



du 

T^xxl 



du 

?t/x2 

^ dy 


du 

'^zx3 



^ du 


T^xxA = ~P 

dx 

i 

4 

^ du 


V5 = ^ 

i 

5 

^ du 


OZ 

i 

6 
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These are the six constituents of the x-direction shear tensor Their net 
effect is a vector expressed as the divergent V • Tix, which, in incompressible 
flow\ is expressed as: 


V • nv = 


dfxx 


+ 


^ dTz. 


= A'- 


dx dy dz 
dx 

«(i) 


+ 


+ 


dx 

d^u 


dy dz 

i + y — - - 1 + y- 


+ 


dy 

d^u d^u 


dz 


+ 


{dxf {dyf {dz) 


(4/29) 


Using the Laplacian operator (see also Appendix A2 p.217) to represent the 
spatial variation of the spatial variation of an object: 


= 


V • V 


(4/30) 

V-VA 


(4/31) 



/ V • VA;, \ 



= 

V-VA, 

(4/32) 

V^A, 1 


1 V • VA J 



we can re-write eq. 4/29 more elegantly and generalize to three dimensions: 

^ ■ fix = yV^u i = yV^u (4/33) 

V -fiy = yV^v j = yV^v (4/34) 

V • = yV^w k = pV^w (4/35) 


The three last equations are grouped together simply as 

V-fij = yV^V (4/36) 


With this new expression, we can come back to the Cauchy equation (eq. 4/28), 
in which we can replace the shear term with eq. 4/36. This produces the 
glorious Navier-Stokes equation for incompressible How. 


p— = pg-Vp + pVV (4/37) 

for all incompressible flows of a Newtonian fluid. 


This formidable equation describes the property fields of all incompressible 
flows of Newtonian fluids. It expresses the acceleration field (left-hand side) 
as the sum of three contributions (right-hand side): those of gravity, gradient 
of pressure, and divergent of shear. The solutions we look for in equation 4/37 
are the velocity (vector) field V = {u, v, w) = fi {x, y, z, t) and the pressure field 
p = f 2 {x,y,z,t), given a set of constraints to represent the problem at hand. 
These constraints, called boundary conditions, may be expressed in terms of 

^When the flow becomes compressible, normal stresses Txx will see an additional term 
related to the change in density within the particle. An additional term, H-jpV (^V • v), would 
appear on the right side of eq. 4/37. This is well outside of the scope of this course. 
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velocities (e.g. the presence of a fixed solid body is expressed as a region for 
which y = 0) or pressure (e.g. a discharge into an open atmosphere may be 
expressed as a region of known constant pressure). 

The combination of equations 4/37 and 4/14 (continuity & Navier-Stokes) may 
not be enough to predict the solution of a given problem, most especially if 
large energy transfers take place within the flow. In that case, the addition of 
an energy equation and an expression of the second law of thermodynamics, 
both in a form suitable for derivative analysis, may be needed to provide as 
many equations as there are unknowns. 

Though it is without doubt charming, equation 4/37 should be remembered 
for what it is really: a three-dimensional system of coupled equations. In 
Cartesian coordinates this complexity is more apparent: 


du du du du 


P9x 


dv dv dv dv 


P9y 


dw dw dw dw 


P9z 


dp 

d^u 


(dxy 

dp 

d^v 

dy ^ 

_{dxf 

dp 



{dxY 


d^u 

d'^v 

{dyf 


d^u 

(4/38) 

d^v 

(4/39) 

(4/40) 


Today indeed, 150 years after it was first written, no general expression 
has been found for velocity or pressure fields that would solve this vector 
equation in the general case. Nevertheless, in this course we will use it 
directly: 


• to understand and quantify the importance of key fluid flow parameters, 
in chapter 6; 

• to find analytical solutions to flows in a few selected cases, in the other 
remaining chapters. 


After this course, the reader might also engage into computational fluid 
dynamics (cfd) a discipline entirely architectured around this equation, and 
to which it purposes to find solutions as fields of discrete values. 

As a finishing remark, we note that when the flow is strictly two-dimensional, 
the Navier-Stokes equation is considerably simplified, shrinking down to the 
system: 


P 

P 


du 

dt 

dv 

dt 


+ u- 


+ u- 


du du 

dp 

d^u d^u 

- -h V — 

dx dy 

= P9x-y^+P 

{dxY ' (dyy 

dv dv 

dp 

d^v d^v 

-h V - 

dx dy 

+ 

1 

II 

1 - 

to 

+ 

to 

■ 


(4/41) 

(4/42) 


Unfortunately, this simplification prevents us from accounting for turbulence, 
which is a strongly three-dimensional phenomenon. 


4.4.3 The Bernoulli equation 

We had made clear in chapter 3 (§3.7) that the Bernoulli equation was ex¬ 
traordinarily limited in scope, so much that starting from the more general 
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integral expressions for energy or momentum conservation was much safer 
in practice. As a reminder of this fact, and as an illustration of the bridges that 
can be built between integral and derivative analysis, it can be instructive to 
derive the Bernoulli equation directly from the Navier-Stokes equation. 

Let us start by following a particle along its path in an arbitrary flow, as 
displayed in fig. 4.5. The particle path is known (condition 5 in §3.7), but its 
speed V is not. 

We are now going to project every component of the Navier-Stokes equation 
(eq. 4/37) onto an infinitesimal portion of trajectory ds. Once all terms have 
been projected, the Navier-Stokes equation becomes simply a scalar equation: 

+ p{V • V)k = pg- Vp + pV^k 
p— • ds + p{V-V)V- ds = pg - ds-Vp - ds + ■ ds 


—> 

Because the velocity vector V of the particle, by definition, is always aligned 
with the path, its projection is always equal to its norm: V • ds = V ds. Also, 
the downward gravity g and the upward altitude z have opposite signs, so 
that g - ds = -p dz; we thus obtain: 

dV dV dp -9- ^ 

p— ds -I- p—^V ds = -ppdz-^ ds -t pV^V • ds 

ot ds ds 


When we restrict ourselves to steady flow (condition 1 in §3.7), the first 
left-hand term vanishes. Neglecting losses to friction (condition 4) alleviates 
us from the last right-hand term, and we obtain: 


dV dp 

p Yds = -pgdz-^ds 
ds ds 

pV dV = -pg dz - dp 


This equation can then be integrated from point 1 to point 2 along the 
pathline: 

P^dz-J^ dp 




Figure 4.5 - Different pathlines in an arbitrary flow. We follow one particle as it 
travels from point 1 to point 2. An infinitesimal path segment is named ds. 

Figure CC-0 o.c. 
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When no work or heat transfer occurs (condition 3) and the flow remains 
incompressible (condition 2), the density p remains constant, so that we 
indeed have returned to eq. 3/18: 

A + ghz + = 0 (4/43) 

Thus, we can see that if we follow a particle along its path, in a steady, 
incompressible, frictionless flow with no heat or work transfer, its change in 
kinetic energy is due only to the result of gravity and pressure, in accordance 
with the Navier-Stokes equation. 


4.5 Change of angular momentum 


There is no differential equation for conservation of angular momentum in 
fluid flow. Indeed, as the size of a considered particle tends to zero, it loses 
its ability to possess (or exchange) angular momentum about its center of 
gravity. Infinitesimal fluid particles are not able to rotate about themselves, 
only about an external point. So, in this chapter, no analysis is carried out 
using the concept of angular momentum. 


4.6 Energy conservation 
and increase in entropy 

This topic is well covered in Anderson [8] 

The last two key principles used in fluid flow analysis can be written together 
in a differential equation similar to the Navier-Stokes equation. 

Once again, we start from the analysis of transfers on an infinitesimal control 
volume. We are going to relate three energy terms in the following form, 
naming them A, B and C for clarity: 


the rate of change 

the net flux 

the rate of work done 

of energy inside 

= of heat into 

+ on the element due to 

the fluid element 

the element 

body and surface forces 

A 

= B 

+ c 


(4/44) 


Let us first evaluate term C. The rate of work done on the particle is the dot 
product of its velocity V and the net force F^et applying to it: 


C = V ■ I-its::: + --- 

\ dw drv 

C = v- [pg- Vp + V • Tij) dW 


WCigllL 

dW 


dW / 




and for the case of a Newtonian fluid in an incompressible flow, this expres¬ 
sion can be re-written as: 


C = V- [pg- Vp + pV^v) dW 


(4/45) 


We now turn to term B, the net flux of heat into the element. We attribute 
this flux to two contributions, the first named Qradiation from the emission 
or absorption of radiation, and the second, named ^conduction to thermal 
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conduction through the faces of the element. We make no attempt to quantify 
Eradiation, simply Stating that 


Eradiation ~ P^radiation d*!^ (4/46) 

in which ^'radiation is the local power per unit mass (in Wkg“') transferee! to the 
element, to be determined from the boundary conditions and flow temperature 
distribution. 


In the x-direction, thermal conduction through the faces of the element 
causes a net flow of heat Econduction,x expressed as a function of the power 
per area q (in W m~^) through each of the two faces perpendicular to x: 


Econductionpc 


dqx , 

qx ~ 1tlx 


dqx 

dx 

dqx 

dx 


dx 
dx dy dz 


d'V 


dy dz 


Summing contributions from all three directions, we obtain: 


(4/47) 

(4/48) 


Econduction ~ Econductionpc + Econduction,}/ + Econduction,z 

dqx 


_ dqy dqz 

dx dy dz 


d'V 


(4/49) 


So we finally obtain an expression for B: 


B — Eradiation "b Econduction 

_ dq, 

/^^radiation 


X ^ dqy 

dx 


dqz 
dy dz 


+ 


d'V 


(4/50) 


Finally, term A, the rate of change of energy inside the fluid element, can be 
expressed as a function of the specific kinetic energy and specific internal 
energy i: 


A = p 


D 

Dt 


('■ 


+ 



d^ 


(4/51) 


We are therefore able to relate the properties of a fluid particle to the principle 
of energy conservation as follows: 



"b 2 ^ j ~ (p^radiation 


dx dy dz j 


+ V • {^pg - Vp + 

(4/52) 


This is a scalar equation — it only has one dimension, and involves the length 
of the velocity vector, V = j ^ 

In an interesting hack, we are able to incorporate an expression for the second 
principle of thermodynamics in this equation simply by expressing the fluxes 
qi as a function of the temperature gradients. Indeed, expressing a heat flux 
as the result of a difference in temperature according to the Fourier law, 

,dT ^ ^ 

qi = -k— (4/53) 

oi 
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we constrain the direction of the heat fluxes and thus ensure that all dissipa¬ 
tive terms resulting in temperature increases cannot be fed back into other 
energy terms, thus increasing the overall entropy. 

Equation 4/53 inserted into eq. 4/52 yields: 


P 2^ ) ~ (/^^radiation + k 


q2j q2j q2j 

+ + 


{dxf {dyY {dzf 


^+V-(pg-Vp + yV^v) 


Making use of a Laplacian operator, we come to: 


P~^ 2^^) ~ (p^radiation + + V ■ {^pg - Vp + 

(4/54) 


This equation has several shortcomings, most importantly because the term 
^radiation Is not expressed in terms of fluid properties, and because p is typ¬ 
ically not independent of the temperature T. Nevertheless, it in principle 
brings closure to our system of continuity and momentum equations, and 
these two influences may be either neglected, or modeled numerically. 


4.7 CFD: the Navier-Stokes equations in practice 

This topic is well covered in Versteeg & Malalasekera [11] 

4.7.1 Principle 

In our analysis of fluid flow from a derivative perspective, our five physical 
principles from §0.7 have been condensed into three equations (often loosely 
referred together to as the Navier-Stokes equations). Out of these, the first 
two, for conservation of mass (4/17) and linear momentum (4/37) in incom¬ 
pressible flows, are often enough to characterize most free flows, and should 
in principle be enough to find the primary unknown, which is the velocity 
field V: 


0 = V • T 
DF ^ 

P-^ = PP - Vp + pV F 

We know of many individual analytical solutions to this mathematical prob¬ 
lem: they apply to simple cases, and we shall describe several such flows in 
the upcoming chapters. However, we do not have one general solution: one 
that would encompass all of them. For example, in solid mechanics we have 
long understood that all pure free fall movements can be described with the 
solution X = xo + Uot and y = yo + vot + \gT, regardless of the particularities 
of each fall. In fluid mechanics, even though our analysis was carried out in 
the same manner, we have yet to find such one general solution — or even to 
prove that one exists at all. 

It is therefore tempting to attack the above pair of equations from the nu¬ 
merical side, with a computer algorithm. If one discretizes space and time in 
small increments 5x, Sy, 5z and 5t, we could re-express the x-component of 
eq. 4/37 as: 
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p 


Su\t 

~5r 


+ u- 


Su\x 

5x 


+ V- 


5u\y 

Sy 


+ w 


Su\z 

Sz 


- P9x 


Sp\x 

5x 


+ P 


S 


5u\x I 
5x \x 


8x 


+ 


Su\ 

~Sy 




+ 


5 ^1 

OZ I2: 

5z 


(4/55) 


If we start with a known (perhaps guessed) initial field for velocity and pres¬ 
sure, this equation 4/55 allows us to isolate and solve for 5u\t, and therefore 
predict what the u velocity field would look like after a time increment 5t. 
The same can be done in the y- and z-directions. Repeating the process, we 
then proceed to the next time step and so on, marching in time, obtaining at 
every new time step the value of u, v and w at each position within our com¬ 
putation grid. This is the fundamental working principle for computational 
fluid dynamics (cfd) today. 


4.7.2 Two problems with CFD 

There are two flaws with the process described above. 

The first flaw is easy to identify: while we may be able to start the calculation 

with a known pressure field, we do not have a mean to predict p (and thus 

SpISx in eq. 4/55) at the next position in time (t + 5t). 

This is not a surprise. Going back to the Cauchy equation (eq. 4/28), we 

notice that we took care of the divergent of shear V • Tij, but not of the 
—^ ^ 

gradient of pressure Vp. We have neither an expression p = f{V), nor an 

equation that would give us a value for dpidt. Cfd software packages deal 

with this problem in the most uncomfortable way: by guessing a pressure 

field p(f+^f), calculating the new velocity field V(^t+st)^ then correcting and 

re-iterating with the help of the continuity equation, until continuity errors 

are acceptably small. This scheme and its implementation are a great source 

of anguish amongst numerical fluid dynamicists. 

The second flaw appears once we consider the effect of grid coarseness. Every 
decrease in the size of the grid cell and in the length of the time step increases 
the total number of equations to be solved by the algorithm. Halving each 
of 5x, Sy, Sz and St multiplies the total number of equations by 16, so that 
soon enough the experimenter will wish to know what maximum (coarsest) 
grid size is appropriate or tolerable. Furthermore, in many practical cases, 
we may not even be interested in an exhaustive description of the velocity 
field, and just wish to obtain a general, coarse description of the fluid flow. 

Using a coarse grid and large time step, however, prevents us from resolving 
small variations in velocity: movements which are so small they fit in between 
grid points, or so short they occur in between time steps. Let us decompose 
the velocity field into two components: one is the average flow {u,v,w), 
and the other the instantaneous fluctuation flow {u',v',w'), too fine to be 
captured by our grid: 


Ui = Ui -h u'l (4/56) 
Ui = 0 (4/57) 
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With the use of definition 4/56, we re-formulate equation 4/38 as follows: 


d{u + u') d{u + u') d{u + u') d{u + u') 

+ {u + U ) - -- h (v + V ) - - - h (w + W ) 


dt 


~ P9x 


dx 
d(p + p') 
dx 


dy 


dz 


+ y 


d^{u + u') ^ d^{u + u') ^ d^{u + u') 


{dxY 


{dyf 


{dzy 


Taking the average of this equation —thus expressing the dynamics of the 
flow as we calculate them with a finite, coarse grid— yields, after some 
intimidating but easily conquerable algebra: 


P 


du _du _du _du 

-Z -1 U- -h V- -h w — 

dt dx dy dz 


+ P 


du' 

f 

dx 


+ v'- 


du' 


dy 


+ w'- 


_du' 

dz 


dp 

= PQx-^+P 
dx 


d^u d^u d^u 
(dxY (dyy {dzy 


(4/58) 


Equation 4/58 is the x-component of the Reynolds-averaged Navier-Stokes 
equation (rans). It shows that when one observes the flow in terms of the 
sum of an average and an instantaneous component, the dynamics cannot be 
expressed solely according to the average component. Comparing eqs. 4/58 
and 4/38 we find that an additional term has appeared on the left side. This 
term, called the Reynolds stress, is often re-written as p du'pi'./dj] it is not zero 
since we observe in practice that the instantaneous fluctuations of velocity 
are strongly correlated. 

The difference between eqs. 4/58 and 4/38 can perhaps be expressed differ¬ 
ently: the time-average of a real flow cannot be calculated by solving for 
the time-average velocities. Or, more bluntly: the average of the solution 
cannot be obtained with only the average of the flow. This is a tremendous 
burden in computational fluid dynamics, where limits on the available com¬ 
putational power prevent us in practice from solving for these fluctuations. 
In the overwhelming majority of computations, the Reynolds stress has to be 
approximated in bulk with schemes named turbulence models. 

Diving into the intricacies of cfd is beyond the scope of our study. Never¬ 
theless, the remarks in this last section should hopefully hint at the fact that 
an understanding of the mathematical nature of the differential conservation 
equations is of great practical importance in fluid dynamics. It is for that 
reason that we shall dedicate the exercises of this chapter solely to playing 
with the mathematics of our two main equations. 
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Exercise sheet 4 (Differential analysis) 

last edited March 17, 2017 


These lecture notes are based on textbooks by White [13], Qengel & al.[16], and Munson & al.[18]. 


Except otherwise indicated, we assume that fluids are Newtonian, and that: 

Pwater = 1 000 kg m“^ patm. = 1 bar; Patm. = 1,225 kgm“^ Tatm. = ll,3°C;patm. = 1,5 • 10‘^Nsm“2; 
g = 9,81 ms“^. Air is modeled as aperfect gas (i?air = 287JK“^ kg“'; pair = 1,4; Cpair = 1005Jkg“^ K“^). 


Continuity equation: 


1 Dp ^ ^ 

—c + v-v = o 

p Dt 

(4/16) 

Navier-Stokes equation for incompressible flow: 


DV ^ 

P Df ^ ^ 

(4/37) 


4.1 Revision questions 

For the continuity equation (eq. 4/16), and then for the incompressible Navier-Stokes 
equation (eq. 4/37), 

1. Write out the equation in its fully-developed form in three Cartesian coordinates; 

2. State in which flow conditions the equation applies. 

Also, in order to revise the notion of total (or substantial) derivative: 

3. Describe a situation in which the total derivative of a property is non-zero although 
the fluid property is independent of time. 

4. Describe a situation in which the total derivative of a property is zero although the 
time rate of change of this property is non-zero. 


4.2 Acceleration field 

Qengel & al. [16] E4-3 

A flow is described with the velocity field V = (0,5 -l- 0,8x)i -l- (1,5 - 0,8y)j (in si units). 
What is the acceleration measured by a probe positioned at (2; 2; 2) at 1 = 3 s ? 
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4.3 Volumetric dilatation rate 


A flow is described by the following field (in si units): 

u = + y‘^ + z 

V = xy + yz + 
w = - z^ + 4 


der. Munson & al. [ 18 ] 6.4 


What is the volumetric dilatation rate field? What is the value of this rate at {2;2;2}? 


4.4 Incompressibility 

Qengel & al. [ 16 ] 9-28 

Does the vector field V = (1,6 + l,8x)i + (1,5 - l,8y)j satisfy the continuity equation for 
two-dimensional incompressible flow? 


4.5 Missing components 

^ ^ Munson & al. [ 18 ] E6.2 + Qengel & al. [ 16 ] 9-4 

Two flows are described by the following fields: 

2 , 2,2 
Ui = X + y z 

vi = xy -h yz -h z 

Wi = ? 

U2 = ax^ -I- by^ cz^ 

V2 = ? 

W 2 = axz -h byz^ 

What must wi and V 2 be so that these flows be incompressible? 


4.6 Another acceleration field 

White [ 13 ] E4.1 

Given the velocity field V = {3t)i + {xz)j -h {ty^)k (si units), what is the acceleration field, 
and what is the value measured at {2;4;6} and t = 5s? 


4.7 Vortex 

A vortex is modeled with the following two-dimensional flow: 


-I- 


Qengel & al. [ 16 ] 9.27 


Verify that this field satisfies the continuity equation for incompressible flow. 


104 









4.8 Pressure fields 

Qengel & al. [ 16 ] E9-13, White [ 13 ] 4.32 & 4.34 

We consider the four (separate and independent) incompressible flows below: 

Vi = {ax + b)i + {—ay + cx)j 

V 2 = {2y)i + {8x)f 

Vs = {ax + bt)i + {cx^ + ey)j 

V4 = Uoil + jy-Uo^f 

The influence of gravity is neglected on the first three fields. 

Does a function exist to describe the pressure field of each of these flows, and if so, what 
is it? 
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Answers 


4.1 1) For continuity, use eqs. 4/2 and 1/8 in equation 4/16. For Navier-Stokes, see 

eqs. 4/38, 4/39 and 4/40 p. 95; 2) Read §4.3 p. 88 for continuity, and 

§4.4.2 p.93 for Navier-Stokes; 3) and 4) see §4.2.2 p. 86. 

4.2 ^ = (0,4 -I- 0,64x)/ + (-1,2 + 0,64y)j. At the probe it takes the value l,68i + 0,08^ 
(length 1,682 ms“^). 

4.3 V • V = -x^ + X - z; thus at the probe it takes the value (^V • = ~4 s“h 

—^ 

4.4 Apply equation 4/18 p.90 to V: the answer is yes. 

4.5 1) Applying equation 4/18: Wi = -3x2 - + f{x,y,ty, 

2) idem, V 2 = -3axy - bzy^ + 

4.6 ^ = {3)i+(3z+y^x)tj+(y^+2xyzt)k. At the probe it takes the value 3j-l-250f-l-496k. 

4.7 Apply equation 4/18 to V to verify incompressibility. 


4.8 Note: the constant (initial) value po is sometimes implicitly written in the unknown 
functions /. 

1) p = -p ^abx + la^x^ + bey + +po + /(tp 

2) p = -p (8x2 ^ gy 2 j 3 ) ^ ^ we cannot 

describe the pressure with a mathematical function; 

4)p=-p ^ + 2L 2l) ~ dxX - Pyy +p0+/(f)- 
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These lecture notes are based on textbooks by White [13], Qengel & al.[16], and Munson & al.[18]. 


5.1 Motivation 


In this chapter we focus on fluid flow within ducts and pipes. This topic 
allows us to explore several important phenomena with only very mod¬ 
est mathematical complexity. In particular, we are trying to answer two 
questions: 

1. What does it take to describe fluid flow in ducts? 

2. How can we quantify pressure losses in pipes and the power necessary 
to overcome them? 



Video: pre-lecture briefing for 
this chapter 

by o.c. (CC-by) 
https://youtu.be/7Ba0vkAx6ds 


5.2 Perfect flow in ducts 


We begin with the simplest possible ducted flow case: a purely hypothetical 
fully-inviscid, incompressible, steady fluid flow in a one-dimensional pipe. 
Since there are no viscous effects, the velocity profile across the duct remains 
uniform (flat) all along the flow (fig. 5.1). 

If the cross-sectional area A is changed, then the principle of mass conserva¬ 
tion (eqs. 0/22, 3/5) is enough to allow us to compute the change in velocity 
u = V-. 


pV\Ai — PV2A2 


(5/1) 


in uniform pipe flow. 
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Figure 5.1 - In viscid fluid flow in a one-dimensional duct. In this purely hypothetical 
case, the velocity distribution is uniform across a cross-section of the duct. The 
average velocity and pressure change with cross-section area, but the total pressure 
po = p + remains constant. 


Figure CC-0 o.c. 


This information, in turn allows us to compute the pressure change between 
two sections of different areas by using the principle of energy conservation. 
We notice that the flow is so simple that the five conditions associated with 
the use of the Bernoulli equation (§3.7) are fulfilled: the flow is steady, 
incompressible, one-dimensional, has known trajectory, and does not feature 
friction or energy transfer. A simple application of eq. 3/18 p. 69 yields: 



(5/2) 


in steady, incompressible, inviscid pipe flow without heat or work transfer. 

Thus, in this kind of flow, pressure increases everywhere the velocity de¬ 
creases, and vice-versa. 

Another way of writing this equation is by stating that at constant altitude, 
the total or dynamic pressure ptotal = po = p + remains constant: 


(5/3) 


po = CSt. 


at constant altitude, in laminar inviscid straight pipe flow. 

We shall, of course, turn immediately to more realistic flows. With vis¬ 
cous effects, one key assumption of the Bernoulli equation breaks down, 
and additional pressure losses Apfriction will occur, which we would like to 
quantify. 





































5.3 Viscous laminar flow in ducts 


5.3.1 The entry zone 

Let us observe the velocity profile at the entrance of a symmetrical duct in a 
steady, laminar, viscous flow (fig. 5.2). 


Irrotational 


Velocity 


Developing 


Fully Developed 



Figure 5.2 - Fluid flow velocity distributions at the entrance of a duct. All throughout 
the entrance region, the core region of the flow is accelerated, and the outer region 
decelerated, even though the flow itself is steady. 

Diagram CC-by-sa by Wikimedia Commons user.Devender Kumar5908 


The fluid enters the duct with a uniform (flat) velocity profile. Because of the 
no-slip condition at the wall (§2.3.1 p.50), the particles in contact with the 
duct walls are immediately stopped. Through viscosity, shear stress effects 
propagate progressively inwards. A layer appears in which viscosity effects 
are predominant, which we name boundary layer, this layer grows until it 
reaches the center of the duct. Past this point, the flow is entirely dictated by 
viscous effects and the velocity profile does not change with distance. The 
flow is then said to be fully developed. 


5.3.2 Laminar viscous flow in a one-dimensional duct 


Once the flow is fully developed (but still steady and laminar), how can 
we describe the velocity and pressure distribution in our one-dimensional 
horizontal duct? We can start with a qualitative description. Because of the 
no-slip condition at the walls, the velocity distribution within any cross-flow 
section cannot be uniform (fig. 5.3). Shear occurs, which translates into a 
pressure decrease along the flow. The faster the flow, and the higher the 
gradient of velocity. Thus, shear within the flow, and the resulting pressure 
loss, both increase when the cross-sectional area is decreased. 

How can we now describe quantitatively the velocity profile and the pressure 
loss? We need a powerful, extensive tool to describe the flow — and we turn 
to the Navier-Stokes equation which we derived in the previous chapter as 
eq. 4/37 p.94: 



Video: a young fluid dynami- 
cist explores the mechanics of a 
granite ball fountain. The aes¬ 
thetics are questionable but the 
physics are enjoyable (see ex¬ 
ercise 5.4): a thin layer of wa¬ 
ter flowing in a laminar flow 
is enough to support a sphere 
weighing several tons. 

by World Travel (styl) 
https://youtu.be/23rl8zOX30k 
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DV 


pg - Vp -h p\V 


(5/4) 


Since we are applying this tool to the simple case of fully-developed, two- 
dimensional incompressible fluid flow between two parallel plates (fig. 5.4), 

















































pressure 


total pressure 



distance 


Figure 5.3 - Viscous laminar fluid flow in a one-dimensional pipe. This time, the 
no-slip condition at the wall creates a viscosity gradient across the duct cross-section. 
This in turn translates into pressure loss. Sudden duct geometry changes such as 
represented here would also disturb the flow further, but the effect was neglected 
here. 

Figure CC-0 o.c. 


we need only two Cartesian coordinates, reproducing eqs. 4/41 & 4/42: 


du 

du 

du 

- -h U - 

+ V — 

dt 

dx 

dy\ 

dv 

dv 

dv 

- “h U - 

+ V — 

dt 

dx 

dy_ 


dp 

d^u d^u 

II 

-h 

{dxY ' (dyY 

dp 

d^v d^v 

= pgy-Yy+r 

1- 

to 

-h 

to 

■ 


(5/5) 

(5/6) 


In this particular flow, we have restricted ourselves to a fully-steady (-^ = 0), 
horizontal (g = gy), one-directional flow {v = 0). When the flow is fully 
developed, = 0 and = 0, and the system above shrinks down to: 


dp 

0 = -/+/i 

ox 

dp 


d^u 

{dyf 


(5/7) 

(5/8) 



Figure 5.4 - Two-dimensional laminar flow between two plates, also called Couette 
flow. We already studied this flow case in fig. 2.4 p.55; this time, we wish to derive 
an expression for the velocity distribution. 


Figure CC-0 o.c. 



































































We only have to integrate equation 5/7 twice with respect to y to come to 
the velocity profile across two plates separated by a height 2H: 



(5/9) 


Now, the longitudinal pressure gradient can be evaluated by working out 
the volume flow rate W for any given width Z with one further integration 
of equation 5/9: 



(5/10) 


In this section, the overall process is more important than the result: by 
starting with the Navier-Stokes equations, and adding known constraints 
that describe the flow of interest, we can predict analytically all of the 
characteristics of a laminar flow. 

5.4 Viscous laminar flow in cylindrical pipes 


We now turn to studying flow in cylindrical pipes, which are widely used; 
first considering laminar flow, and then expanding to turbulent flow. 

5.4.1 Principle 

The process is identical to above, only applied to cylindrical instead of 
Cartesian coordinates. We focus on the fully-developed laminar flow of a 
fluid in a horizontal cylindrical pipe (fig. 5.5). 

For this flow, we wish to work out the velocity profile and calculate the 
pressure loss related to the flow. 



Figure 5.5 - A cylindrical coordinate system to study laminar flow in a horizontal 
cylindrical duct. 


Figure CC-0 o.c. 
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5.4.2 Velocity profile 


We once again start from the Navier-Stokes vector equation, choosing this 
time to develop it using cylindrical coordinates: 


dv, 


r uur V0 dVy dVr 
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(5/13) 


This mathematical arsenal does not frighten us, for the simplicity of the flow 
we are studying allows us to bring in numerous simplifications. Firstly, we 
have Vy = 0 and ve = 0 everywhere. Thus, by continuity, dvzidz = 0. 
Furthermore, since our flow is symmetrical, Vz is independent from 9. With 
these two conditions, the above system shrinks down to: 


0 = 
0 = 
0 = 


dp 

P^"-yr 

1 dp 
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P9e- 
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i_5 
r dr 



(5/14) 

(5/15) 

(5/16) 


Equations (5/14) and (5/15) hold the key to describing pressure: 


P = -P9y + f{z) 


(5/17) 


this expresses the (unsurprising) fact that the pressure distribution within a 
vertical cross-section of the pipe is simply the result of a hydrostatic gradient. 

Now, with equation 5/16, we work towards obtaining an expression for Vz by 
integrating twice our expression for dvz/dr: 


d / dvz\ 
Yr V'^) 



Vz 


r dp 
p dz 

2 p \dzj 
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+ ki 

+ ki In r + ^2 


(5/18) 
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We have to use boundary conditions so as to unburden ourselves from 
integration constants ki and ^ 2 - 
































By setting Vz (r=o) as finite, we deduce that ki = 0 (because ln(0) ^ -oo). 

By setting (r=R) = 0 (no-slip condition), we obtain ^2 = 

This simplifies eq. (5/18) and brings us to our objective, an extensive ex¬ 
pression for the velocity profile across a pipe of radius R when the flow is 
laminar; 


(i) 

This equation is parabolic (fig. 5.6). It tells us that in a pipe of given length L 
and radius R, a given velocity profile will be achieved which is a function 
only of the ratio Ap/p. 



Figure 5.6 - The velocity profile across a cylindrical pipe featuring laminar viscous 
flow. 

Figure CC-0 o.c. 


In the same way we proceeded before, we can express the pressure gradient 
in the pipe as a function of the volume flow rate, which itself is obtained 
through integration: 

O' 

dp _ Apfi-iction 
dx ~ L 

This equation is interesting in several respects. For a given pipe length L and 
pressure drop ApfI.iction^ the volume flow 'V increases with the power 4 of the 
diameter D. In other words, the volume flow is multiplied by 16 every time 
the diameter is doubled. 

We also notice that the pipe wall roughness does not appear in equation 5/20. 
In a laminar flow, increasing the pipe roughness has no effect on the velocity 
of the fluid layers at the center of the pipe. 

When the pipe is inclined relative to the horizontal with an angle a, eqs. 
(5/19) and (5/20) remain valid except for the pressure drop term Ap which is 
simply replaced by (Apfriction + pgL sin a), since gravity then contributes to 
the longitudinal change in pressure. 

5.4.3 Quantification of losses 

Hydraulics is the oldest branch of fluid dynamics, and much of fhe nofation 
used to describe pressure losses predates modern applications. The most 
widely-used parameters for quantifying losses due to friction in a duct are 
the following: 


/ \ j (dp 

Jo 128 p \dx 

128 p 


n D4 




(5/20) 
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The elevation loss noted |Ae| is defined as 


i 


Video: very basic visualization 
of laminar and turbulent flow 
regimes in a transparent pipe 

by Engineering Fundamentals (styl) 
https://youtu.be/56AyTIhNQBo 


I A I I^Pfrictionl 

|Ae| = - 

pg 


(5/21) 


It represents the hydrostatic height loss (with a positive number) caused 
by the fluid flow in the duct, and is measured in meters. It is often writ¬ 
ten Ah (without bars, but as a positive number), but in this document 
we reserve h to denote enthalpy (as in eq. 3/17 p.68). 


The Darcy friction factor noted f is defined as 

I ^Pfriction I 
--oV^ 

D 2 ” av. 

where Vav. is the average flow velocity in the pipe. 


(5/22) 


This dimensionless term is widely used in the industry because it 
permits the quantification of losses independently of scale effects (L/D) 
and inflow conditions (pV^^). 


The loss coefficient noted Ki is defined as 


Kl = 


I /^pfriction 
-pV^ 

2r^a.Y. 


(5/23) 


It is simply equivalent to a Darcy coefficient in which no value of D/L 
can be measured, and it is used mostly in order to describe losses in 
single components found in pipe systems, such as bends, filter screens, 
or valves. 


For a laminar flow in a cylindrical duct, for example, we can easily compute /, 
inserting equation 5/20 into the definition 5/22: 


fiaminar cylinder flow 


32V^v.ph ^ . P 

- — o4- 


in which we inserted the average velocity Vav. = = 


Ap d2 
32/iL • 


(5/24) 


5.5 Viscous turbulent flow in cylindrical pipes 


5.5.1 Predicting the occurrence of turbulence 


It has long been observed that pipe flow can have different regimes. In 
some conditions, the flow is unable to remain laminar (one-directional, fully- 
steady); it becomes turbulent. Although the flow is steady when it is averaged 
over short time period (e.g. a few seconds), it is subject to constant, small- 
scale, chaotic and spontaneous velocity field changes in all directions. 

In 1883, Osborne Reynolds published the results of a meticulous investigation 
into the conditions in which the flow is able, or not, to remain laminar 
(figs. 5.7 and 5.8). He showed that they could be predicted using a single 
non-dimensional parameter, later named Reynolds number, which, as we 
have seen already with eq. 0/27 p.21, is expressed as: 


[Re] = 


pVL 

P 
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(5/25) 











Figure 5.7 - Illustration published by Reynolds in 1883 showing the installation 
he set up to investigate the onset of turbulence. Water flows from a transparent 
rectangular tank down into a transparent drain pipe, to the right of the picture. 
Colored die is injected at the center of the pipe inlet, allowing for the visualization 
of the flow regime. 

Image by Osborne Reynolds (1883, public domain) 



Laminar flow 





Turbulent flow 




— ■ 

Turbulent flow (observed with an electric spark) 


Figure 5.8 - Illustration published by Reynolds in 1883 showing two different flow 
regimes observed in the installation from fig. 5.7. 

Image by Osborne Reynolds (1883, public domain) 
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In the case of pipe flow, the representative length L is conventionally set to 
the pipe diameter D and the velocity to the cross-section average velocity: 

[Rein = (5/26) 

where Vav, is the average velocity in the pipe (m s“'), 
and D is the pipe diameter (m). 

The occurrence of turbulence is very well documented (although it is still not 
well-predicted by purely-theoretical models). The following values are now 
widely accepted: 

• Pipe flow is laminar for [Re]^ < 2 300 ; 

• Pipe flow is turbulent for [ReJ^ > 4 000. 

The significance of the Reynolds number extends far beyond pipe flow; we 
shall explore this in chapter 6 (pp.l34 & 139). 

5.5.2 Characteristics of turbulent flow 

This topic is well covered in Tennekes & Lumley [4] 

Turbulence is a complex topic which is still not fully described analytically 
today. Although it may display steadiness when time-averaged, a turbulent 
flow is highly three-dimensional, unsteady, and chaotic in the sense that 
the description of its velocity field is carried out with statistical, instead of 
analytic, methods. 

In the scope of our study of fluid dynamics, the most important characteristics 
associated with turbulence are the following: 

• A strong increase in mass and energy transfer within the flow. Slow and 
rapid fluid particles have much more interaction (especially momentum 
transfer) than within laminar flow; 

• A strong increase in losses due to friction (typically by a factor 2). The 
increase in momentum exchange within the flow creates strong dissipa¬ 
tion through viscous effects, and thus transfer (as heat) of macroscopic 
forms of energy (kinetic and pressure energy) into microscopic forms 
(internal energy, translating as temperature); 

• Internal flow movements appear to be chaotic (though not merely 
random, as would be white noise), and we do not have mathematical 
tools to describe them analytically. 

Consequently, solving a turbulent flow requires taking account of flow in all 
three dimensions even for one-directional flow! 

5.5.3 Velocity profile in turbulent pipe flow 

In order to deal with the vastly-increased complexity of turbulent flow in 
pipes, we split each velocity component vi in two parts, a time-averaged 
component vl and an instantaneous fluctuation v’^: 

Vr = + v'^ 

V0 = V0+v'q 

3^2 = + t;' 
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In our case, and vq are both zero, but the fluctuations v'^ and v'q are not, and 
will cause Vz to differ from the laminar flow case. The extent of turbulence is 
often measured with the concept of turbulence intensity I: 


I = 



Vi 


(5/27) 


Regrettably, we have not found a general analytical solution to turbulent pipe 
flow — note that if we did, it would likely exhibit complexity in proportion 
to that of such flows. A widely-accepted average velocity profile constructed 
from experimental observations is: 


U(^) — Vz — Vz max 



(5/28) 


While it closely and neatly matches experimental observations, this model 
is nowhere as potent as an analytical one and must be seen only as an 
approximation. For example, it does not allow us to predict internal energy 
dissipation rates (because it describes only time-averaged velocity), or even 
wall shear stress (because it yields {duldr)r=R = oo, which is not physical). 

The following points summarize the most important characteristics of turbu¬ 
lent velocity profiles: 

• They continuously fluctuate in time and we have no means to predict 
them extensively; 

• They are much “flatter” than laminar profiles (fig. 5.9); 

• They depend on the wall roughness; 

• They result in shear and dissipation rates that are markedly higher 
than laminar profiles. 



ABC 


Figure 5.9 - Velocity profiles for laminar (A), and turbulent (B and C) flows in 
a cylindrical pipe. B represents the time-averaged velocity distribution, while C 
represents an arbitrary instantaneous distribution. Turbulent flow in a pipe also 
features velocities in the radial and angular directions, which are not shown here. 

Figure CC-0 o.c. 


5.5.4 Pressure losses in turbulent pipe flow 

Losses caused by turbulent flow depend on the wall roughness e and on the 
diameter-based Reynolds number [Re]^). 

For lack of an analytical solution, we are not able to predict the value of the 
friction factor f anymore. Several empirical models can be built to obtain /, 
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the most important of which is known as the Colebrook equation expressed 
as: 


1 , / 1 e 2,51 \ 


(5/29) 


The structure of this equation makes it inconvenient to solve for f. To 
circumvent this difficulty, equation 5/29 can be solved graphically on the 
Moody diagram, fig. 5.10. This classic document allows us to obtain numerical 
values for f (and thus predict the pressure losses Apfriction) if we know the 
diameter-based Reynolds number [Re]^) and the relative roughness e/D. 



Reynolds Number, Re = 


Figure 5.10 - A Moody diagram, which presents values for f measured experi¬ 
mentally, as a function of the diameter-based Reynolds number [Re]^), for different 
values of the relative roughness e/D. This figure is reproduced with a larger scale as 
figure 5.12 p.l20. 


Diagram CC-by-sa S Beck and R Collins, University of Sheffield 














































































































































































Exercise sheet 5 (Pipe flow) 

last edited April 4, 2017 


These lecture notes are based on textbooks by White [13], Qengel & al.[16], and Munson & al.[18]. 


Except otherwise indicated, we assume that fluids are Newtonian, and that: 

Pwater = 1 000 kg m“^ patm. = 1 bar; Patm. = 1,225 kgm“^ Tatm. = ll,3°C;patm. = 1,5 • 10‘^Nsm“2; 
g = 9,81 ms“^. Air is modeled as aperfect gas (i?air = 287JK“^ kg“'; pair = 1,4; Cpair = 1005Jkg“^ K“^). 


In cylindrical pipe flow, we accept the flow is always laminar for [ReJ^ < 2 300, and 
always turbulent for [Re]^) > 4 000. The Darcy friction factor / is defined as: 




|Ap| 


--oV 


ai. 


The loss coefficient Ki is defined as: 


(5/22) 


Kl = 



(5/23) 


Viscosities of various fluids are given in fig. 5.11. Pressure losses in cylindrical pipes 
can be calculated with the help of the Moody diagram presented in fig. 5.12 p.l20. 


o 



XXX 
jUi/s . ^ Ajisoosta siniosqy 


Figure 5.11 - Viscosity of various fluids at a pressure of 1 bar (in practice viscosity is almost 
independent of pressure). 

Figure © White 2008 [13] 
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Moody Diagram 
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Figure 5.12 - A Moody diagram, which presents values for f measured experimentally, as a 
function of the diameter-based Reynolds number [Re]/), for different relative roughness values. 

Diagram CC-by-sa S Beck and R Collins, University of ShefReld 
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5.1 Revision questions 

non-examinable 

The Moody diagram (fig. 5.12 p.l20) is simple to use, yet it takes practice to understand 
it fully... here are three questions to guide your exploration. They can perhaps be 
answered as you work through the other examples. 

1. Why is there no zero on the diagram? 

2. Why are the curves sloped downwards — should friction losses not instead increase 
with increasing Reynolds number? 

3. Why can the pressure losses Ap be calculated given the volume flow 'V, but not 
the other way around? 


5.2 Air flow in a small pipe 

Munson & al. [IS] E8.5 

A machine designed to assemble micro-components uses an air jet. This air is driven 
through a 10 cm-long cylindrical pipe with a 4 mm diameter, roughness 0,0025 mm, at an 
average speed of 50 m s“^. 

The inlet air pressure and temperature are 1,2 bar and 20 °C; the air viscosity is quantified 
in fig. 5.11 p. 119. 

1. What is the pressure loss caused by the flow through the pipe? 

2. Although that is not possible in practice, what would the loss in the case where 
laminar flow could be maintained throughout the pipe? 


5.3 Couette flow 

CC-O o.c. 

We consider the laminar flow of a fluid between two parallel plates (named Couette flow), 
as shown in fig. 5.13. 



Figure 5.13 - Two-dimensional laminar flow between two plates, also called Couette How. 

Figure CC-O o.c. 


1. Starting from the Navier-Stokes equations for incompressible, two-dimensional 
flow. 
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show that the velocity profile in a horizontal, laminar, fully-developed flow between 
two horizontal plates separated by a gap of height 2H is: 


u = — (^) 

2y\dxl^^ ’ 


(5/9) 


2. Why would this equation fail to describe turbulent flow? 


5.4 Kugel fountain 


Derived from Munson & al. [is] 6.91 


A Kugel fountain is erected to entertain students of fluid mechanics (figs. 5.14 and 5.15). 
A granite block (2 700kgm“^) is sculpted into a very smooth sphere of diameter 1,8 m. 
The sphere is laid on a cylindrical concrete stand of internal diameter 1,2 m, whose edges 
are also smoothed out. 


The stand is filled with water, and a pump creates a flow which lifts the sphere to make 
a decorative fountain. Water flows between the stand and the sphere along a length 
of 10 cm, with a thickness of 0,13 mm. Entrance effects are negligible, and because the 
Reynolds number is low, the flow is entirely laminar. 


1. What is the pressure exerted by the sphere in the central part of the fountain? 

2. The velocity distribution for flow in between two flat plates of width Z separated by 
an interval 2H, known as Couette flow, is modeled as a function of the stream-wise 
pressure gradient dp/dx with the equation: 


3. 

4. 




(5/9) 


Starting from this equation 5/9, show that the pressure gradient can be quantified 
as a function of the volume flow 'V: 


dp 

dx 


3 p 


2ZH3 




(5/10) 


What is the pressure drop across the interval between the sphere and the concrete 
foundation? 


What is the pumping power required for the fountain to work? 




Figure 5.14 - Schematic layout for a Kugel fountain. A pump (bottom left) increases the water 
pressure to a value high enough that it can support the weight of the sphere. A thin strip of water 
flows between the concrete foundation and the sphere. 
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Figure CC-0 o.c. 





























5. How would this power change if the diameter of the concrete support (but not of 
the sphere) was increased? (briefly justify your answer, e.g. in 30 words or less) 



Figure 5.15 - A Kugel fountain, featuring questionable aesthetics but providing the opportunity 
for amusing interaction. A granite sphere is supported by the flow of water between the sphere 
and a solid strip of concrete. Photo CC-by-sa by Commons User.Atamari 


5.5 Water piping 

A long pipe is installed to carry water from one large reservoir to another (fig. 5.16). The 
total length of the pipe is 10 km, its diameter is 0,5 m, and its roughness is e = 0,5 mm. It 
must climb over a hill, so that the altitude changes along with distance. 



Figure 5.16 - Layout of the water pipe. For clarity, the vertical scale is greatly exaggerated. The 
diameter of the pipe is also exaggerated. 

Figure CC-0 o.c. 


The pump must be powerful enough to push 1 m^ s ^ of water at 20 °C. 

Figure 5.11 p.ll9 quantifies the viscosity of various fluids, and fig. 5.12 p. 120 quantifies 
losses in cylindrical pipes. 

1. Will the flow in the water pipe be turbulent? 

2. What is the pressure drop generated by the water flow? 

3. What is the pumping power required to meet the design requirements? 

4. What would be the power required for the same volume flow if the pipe diameter 
was doubled? 
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5.6 Pipeline 

^ CC-O o.c. 

The Trans-Alaska Pipeline System is a cylindrical smooth steel duct with 1,22 m diameter, 

average roughness e = 0,15 mm, and length 1200 km. Approximately 700 thousand 
barrels of oil (110 000 m^) transit through the pipe each day. 

The density of crude oil is approximately 900 kg m“^ and its viscosity is quantified in 
fig. 5.11 p. 119. The average temperature of the oil during the transit is 60 °C. In industrial 
pumps, oil starts to cavitate (change state, a very undesirable behavior) when its pressure 
falls below 0,7 bar. 

The pipeline is designed to withstand ground deformations due to seismic movements in 
several key zones, and crosses a mountain range with a total altitude variation of 1 400 m. 

1. How much time does an oil particle need to travel across the line? 

2. Propose a pumping station arrangement, and calculate the power required for each 
pump. 

3. How would the pumping power change if the speed was increased? 


5.7 Pump with pipe expansion 


A pump is used to carry a volume flow of 200 L s ^ from one large water reservoir to 
another (fig. 5.17). The altitude of the water surface in both reservoirs is the same. 



Figure 5.17 - Layout of the water pipe. For clarity, the diameter of the pipe and the vertical scale 
are exaggerated. 

Figure CC-O o.c. 


The pipe connecting the reservoirs is made of concrete (e = 0,25 mm); it has a diameter 
of 50 cm on the first half, and 100 cm on the second half. In the middle, the conical 
expansion element induces a loss coefficient of 0,8. At the outlet (at point D), the pressure 
is approximately equal to the corresponding hydrostatic pressure in the outlet tank. 

The inlet is 14 m below the surface. The total pipe length is 400 m; the altitude change 
between inlet and outlet is 12 m. 

1. Represent qualitatively (that is to say, showing the main trends, but without 
displaying accurate values) the water pressure as a function of pipe distance, when 
the pump is turned off. 

2. On the same graph, represent qualitatively the water pressure when the pump is 
switched on. 
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3. What is the water pressure at points A, B, C and D? 












5.8 A more complex ducted flow 

non-examinable. From institute archives 

A more complicated laminar flow case can be studied with the following case. 

In the middle of a vertical container of width 2h filled with oil (fig. 5.18), a plate of width b 
(perpendicular to the plane) and length L with negligible thickness is sinking at constant 
velocity Vp. 

Under the assumption that the flow is laminar everywhere, steady and fully developed, 
give an analytical expression for the friction force applied on the plate. 


cn 







Figure 5.18 - The hypothetical case of a plate sinking vertically in an oil reservoir. 


5.9 Wind tunnel 

non-examinable 

Describe the main characteristics of a wind tunnel that could be installed and operated 
in the room you are standing in. 

In order to do this: 

• Start by proposing key characteristics for the test section; 

• From these dimensions, draw approximately an air circuit to feed the test section 
(while attempting to minimize the size of the fan, whose cost increases exponen¬ 
tially with diameter). 

• Quantify the static and stagnation pressures along the air circuit, by estimating the 
losses generated by wall shear and in the bends (you may use data from fig. 5.19); 

• Quantify the minimum power required to generate your chosen test section flow 
characteristics. 
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T1 = O.IIO T| = 0.138 ^ = 0.200 


Filter screen: rj = 0,05 


Figure 5.19 - Loss coefficients Kl (here noted rj) generated by the use of various components 
within wind tunnel ducts. 


Figure © Barlow, Rae & Pope 1999 [9] 


5.10 Politically incorrect fluid mechanics 


non-examinahle 


In spite of the advice of their instructor, a group of students attempts to apply fluid 
mechanics to incommendable activities. Their objective is to construct a drinking straw 
piping system that can mix a drink of vodka and tonic water in the correct proportions 
(fig. 5.20). They use a “Strawz” kit of connected drinking straws, two bottles, and a glass 
full of ice and liquid water to cool the mix. 

For simplicity, the following information is assumed about the setup: 

• Vodka is modeled as 40 % pure alcohol (ethanol) with 60 % water by volume; 

• Ethanol density is 0,8 kg m“^; 

• Tonic water is modeled as pure water; 

• The viscosity of water and alcohol mixes is described in table 5.1 (use the nearest 
relevant value); 

• The pipe bends induce a loss coefficient factor Kibend = 0^5 each; 

• The pipe T-junction induces a loss coefficient factor Ki = 0,3 in the line direction 
and 1 in the branching flow; 



vodka 


tonic water 


ice water 


Figure 5.20 - Conceptual sketch of a student experiment. 
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Figure CC-0 o.c. 



























Percentage of alcohol by weight 

Viscosity in centipoise 

0 

1,005 

10 

1,538 

20 

2,183 

30 

2,71 

40 

2,91 

50 

2,87 

60 

2,67 

70 

2,37 

80 

2,008 

90 

1,61 

100 

1,2 


Table 5.1 - Viscosity of a mix of ethanol and water at 20 °C. 

data from Bingham, Trans. Chem. Soc, 1902, 81, 179. 


• The pipe has inner diameter D = 3 mm and roughness rj = 0,0025 mm. 

The students wish to obtain the correct mix: one quarter vodka, three quarters tonic 
water. For given levels of liquid in the bottles, is there a straw pipe network configuration 
that will yield the correct mix, and if so, what is it? 
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Answers 


5.2 1) Calculating inlet density with the perfect gas model, [ReJ^ = 14 263 (turbulent), 
a Moody diagram read gives / « 0,029, so Apfnction = “1 292,6 Pa = -0,0129 bar. 

2) If the flow were (magically) kept laminar, with equation 5/24, Apfnction = 200 Pa. 

5.3 The structure is given in the derivation of equation 5/9 p. Ill, and more details 
about the math are given in the derivation of the (very similar) equation 5/19 p. 113. 

5.4 Ap = =71,51 kPa. Then 'V = 0,494 L s“^ which enables us to obtain Vl^ump = 

35,3 W. 

5.5 lApaitl = pp(26-8+5-7) = l,57bar and |Apfrictionl = 51,87 bar: VVpump = 5,345 MW. 

5.6 The total pumping power is W = 10 MW (!). Be careful not to cavitate the oil in 
the ascending sections, and not to burst the pipe in the descending sections! 

5.7 Pressure losses to friction are -3 735 Pa in the first half, +71 Pa in the throat, 
-117 Pa in the second half. Add hydrostatic pressure to obtain final result. 

5.8 Calculate the velocity distribution in the same way as for equation 5/9 p.lll. 
Evaluate the pressure gradient ^ using mass conservation (total cross-section 
mass flow is zero). With the full velocity distribution, derivate v with respect to x 
to obtain E,- = pSgh + ^Vp. 

5.9 This is a fun (and quick) exercise, but the results strongly depend on your proposed 
design! Get help during office hours. 

5.10 The author cannot remember which exercise you are referring to. 
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These lecture notes are based on textbooks by White [13], Qengel & al.[16], and Munson & al.[18]. 


6.1 Motivation 


In this chapter, we use the tools we derived in chapters 3 and 4 to analyze 
scale effects in fluid mechanics. This study should allow us to answer two 
questions: 

• How can a flow be adequately reproduced in a smaller or larger scale? 

• How do forces and powers change when the flow is scaled? 


6.2 Scaling forces 


6.2.1 Force coefficients 


Let us begin as a designer working for an automobile manufacturer. We wish 
to investigate the air flow around a full-size car, using a small model in a wind 
tunnel. The question we seek to answer is: if we succeed in “reproducing the 
exact real flow” in the wind tunnel, how large (or small) will the aerodynamic 
forces be on the model? 

For this, we wish to have a sense of how flow-induced forces scale when 
fluid flows are scaled. A look back on chapter 3, and in particular eq. 3/8 p.65, 
reminds us that for a steady flow through a given control volume, the net 
force induced on the fluid is expressed by: 


net 


"fL‘ 


pV (Vrel ■n)dA = 





Video: pre-lecture briefing for 
this chapter 

by o.c. (CC-by) 
https://youtu.be/ZMlmOgXPruI 


( 6 / 1 ) 
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This equation tells us that the norm of the force vector fnet is directly related 
to the norm of the vector T p\Vx\JW. In the selection of a scale by which to 
measure ^’net^ it is therefore sensible to include a term proportional to the 
the density p, a term proportional to the area A, and a term proportional to 
the square of velocity V. This “scale of fluid-induced force” is conventionally 
measured using the force coefficient Cp: 


Cp = 


F 


where F is the considered fluid-induced force (N); 
p is a reference fluid density (kg m“^); 

5 is a reference surface area (nr); 
and y is a reference velocity (m s“^). 


( 6 / 2 ) 


In effect, the force coefficient relates the magnitude of the force exerted by 
the fluid on an object (F) to the rate of flow of momentum towards the object 
{pSV^ = mV). 

It is worth making a few remarks about this equation. First, it is important 
to realize that eq. 6/2 is a definition: while the choice of terms is guided 
by physical principles, it is not a physical law in itself, and there are no 
reasons to expect Cp (which has no dimension, and thus no unit) to reach 
any particular value in any given case. The choice of terms is also worth 
commenting: 


• F can be any fluid-induced force; generally we are interested in quan¬ 
tifying either the drag Fq (the force component parallel to the free- 
stream velocity) or the lift Fp (force component perpendicular to the 
free-stream velocity); 

• the reference area S is chosen arbitrarily. It matters only that this 
area grow and shrink in proportion to the studied flow case. In the 
automobile industry, it is customary to choose the vehicle frontal area 
as a reference, while in the aeronautical industry the top-view wing 
area is customarily used. The square of a convenient reference length L 
can also be chosen instead; 

• the choice of reference velocity V and density p is also arbitrary, since 
both these properties may vary in time and space within the studied 
fluid flow case. It matters merely that the chosen reference values are 
representative of the case; typically the free-stream (faraway) condi¬ 
tions I4o and p^o are used; 

• the term 1/2 in the denominator is a purely arbitrary and conventional 
value. 

Force coefficients are meaningful criteria to compare and relate what is going 
on in the wind tunnel and on the full-size object: in each case, we scale the 
measured force according to the relevant local flow conditions. 


6.2.2 Power coefficient, and other coefficients 

During our investigation of the flow using a small-scale model, we may be 
interested in measuring not only forces, but also other quantities such as 
power — this would be an important parameter, for example, when studying 
a pump, an aircraft engine, or a wind turbine. 
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Going back once again to chapter 3, and in particular eq. 3/15 p.68, we recall 
that power gained or lost by a fluid flowing steadily through a control volume 
could be expressed as: 


dE. 


sys 


dt 


ff pe{V,,i-n)dA = J]{{p\V^\A)e} 
JJCS net 


in which the specific energy e grows proportionally to 


(6/3) 


Thus, the power gained or lost by the fluid is directly related to the magnitude 
of the scalar S p\V±\AV^. A meaningful ’’scale” for the power of a machine 
can therefore be the amount of energy in the fluid that is made available to it 
every second, a quantity that grows proportionally to pAV^. 

This “scale of fluid flow-related power” is conventionally measured using the 
power coefficient Cp: 


Cp = 


W 

1 /^ 


’where W is the po’wer added to or subtracted from the fluid (W). 


(6/4) 


Other quantities related to fluid flow, such as pressure loss or shear friction, 
are measured using similarly-defined coefficients. We have already used the 
pressure loss coefficient Ki = \Ap\/^pV^^ in chapter 5 (eq. 5/23 p.ll4), and we 
shall soon use the shear coefficient Cf in the coming chapter (eq. 7/5 p.l51). 


6.2.3 Physical similarity 

In the paragraphs above, we have already hinted that under certain condi¬ 
tions, the force and power coefficients in experiments of differing scales may 
have the same value. This happens when we have attained physical similarity. 
Several levels of similarity may be achieved between two experiments: 

• geometrical similarity (usually an obvious requirement), by which we 
ensure similarity of shape between the solid boundaries, which are 
then all related by a scale factor r; 

• kinematic similarity, by which we ensure that the motion of all solid 
and fluid particles is reproduced to scale — thus that for a given time 
scale ratio rt the relative positions are scaled by r, the relative velocities 
by r/rf and the relative accelerations by (r/rf)^; 

• dynamic similarity, by which we ensure that not only the movement, 
but also the forces are all scaled by the same factor. 

In more complex cases, we may want to observe other types of similarity, 
such as thermal or chemical similarity (for which ratios of temperature or 
concentrations become important). 

The general rule that we follow is that the measurements in two experiments 
can be meaningfully compared if the experiments are physically similar. A 
direct consequence for our study of fluid mechanics is that in two different 
experiments, the force and power coefficients will be the same if the flows are 
physically similar. 

In the design of experiments in fluid mechanics, considerable effort is spent to 
attain at least partial physical similarity. Once this is done, measurements can 
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volume compared to 100% 



Video: practical application of 
scale effects: translating mea¬ 
surements made a 60 % size for¬ 
mula one car wind tunnel model 
into their “real size” race car val¬ 
ues, as in exercise 6.4 

by the Sauber FI team (styl) 
https://youtu.be/FQxmOQmnaGw 



Video: when their movement is 
filmed and then viewed sped-up, 
fog banks can appear to flow 
like water. What time-lapse rate 
is required to achieve physical 
similarity? 

by Simon Christen (stvl) 
https://vimeo.com/69445362 


easily be extrapolated from one case to the other because all forces, velocities, 
powers, pressure differences etc. can be connected using coefficients. 

For example, a flow case around a car A may be studied using a model B. 
If dynamic similarity is maintained (and that is by no means an easy task!), 
then the flow dynamics will be identical. The drag force Fb b measured on 
the model can then be compared to the force Fd a on the real car, using 
coefficients: since we have: 


Fda = = Cfd^-PaSaVI = 

-lpA^Yl\r 

IpbSbV?) 


Fd] 


IpbSbV^ 2 


PaSa^a 


It should be noted that even when perfect dynamic similarity is achieved, 
fluid-induced forces do not scale with other types of forces. For example, if 
the length of model B is 1/lOth that of car A, then Sb = 0,1^Sa and at identical 
fluid density and velocity, we would have Fd b = 0,1^Fd a- Nevertheless, the 
weight Fw b of model B would not be one hundredth of the weight of car A. 
If the same materials were used to produce the model, weight would stay 
proportional to volume, not surface, and we would have Fwb = OT^Fwa- 
This has important consequences when weight has to be balanced by flow- 
induced forces such as aerodynamic lift on an airplane. It is also the reason 
why large birds such as condors or swans do not look like, and cannot fly as 
slowly as mosquitoes and bugs!^ 


6.3 Scaling flows 


We have just seen that when two flows are dynamically similar, their compar¬ 
ison is made very easy, as measurements made in one case can be transposed 
to the other with the use of coefficients. This naturally raises the question: 
how can we ensure that two flows are flows are dynamically similar? The 
answer comes from a careful re-writing of the Navier-Stokes equation: two 
flows are dynamically similar when all the relevant flow coefficients have 
the same value. 

6.3.1 The non-dimensional Navier-Stokes equation 

In this section we wish to obtain an expression for the Navier-Stokes equation 
for incompressible flow that allows for an easy comparison of its constituents. 
We start with the original equation, which we derived in chapter 4 as eq. 4/37 
p.94: 

+ p{V-V)^ = pg-Vp + (6/5) 
ot 

What we would now like to do is separate the geometry from the scalars in 

this equation. This is akin to re-expressing each vector A as the multiple of 

—^ 

its length A and a non-dimensional vector A*, which has the same direction 
as A but only unit length. 
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'These ideas are beautifully and smartly explored in Tennekes [7, 14]. 










In order to achieve this, we introduce a series of non-dimensional physical 
terms, starting with non-dimensional time t*, defined as time t multiplied by 
the frequency / at which the flow repeats itself: 

t* = ft (6/6) 

A flow with a very high frequency is highly unsteady, and the changes in 
time of the velocity field will be relatively important. On the other hand, 
flows with very low frequencies are quasi-steady. In all cases, as we observe 
the flow, non-dimensional time t* progresses from 0 to 1, after which the 
solution is repeated. 

We then introduce non-dimensional velocity V*, a unit vector field equal to 
the velocity vector field divided by its own (scalar field) length V: 

V 

V*^- (6/7) 

Pressure p is non-dimensionalized differently, since in fluid mechanics, it is 

the pressure changes across a field, not their absolute value, that influence 

—^ ^ 

the velocity field. For example, in eq. 6/5 Vp can be replaced by V(p - poo) 
(in which poo can be any constant faraway pressure). Now, the pressure 
field p - poo is non-dimensionalized by dividing it by a reference pressure 
difference po - poo, obtaining: 


_ P poo 

p =- 

po poo 


( 6 / 8 ) 


If po is taken to be the maximum pressure in the studied field, then p* is a 
scalar field whose values can only vary between 0 and 1. 

Non-dimensional gravity g* is simply a unit vector equal to the gravity 
vector g divided by its own length g: 


g 


(6/9) 


And finally, we define a new operator, the non-dimensional del V*, 

V* = L V (6/10) 

which ensures that vector fields obtained with a non-dimensional gradient, 
and the scalar fields obtained with a non-dimensional divergent, are “scaled” 
by a reference length L. 

These new terms allow us to replace the constituents of equation 6/5 each by 
a non-dimensional “unit” term multiplied by a scalar term representing its 
length or value: 

t* 

^“7 

p = pp* 

p ~ poo — p {po ~ poo) 

9 = g g 

^ 1 ^ 

V = -V* 

L 


Now, inserting these in equation 6/5, and re-arranging, we obtain: 
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{vv*)+p{vv* • ^vj {vv*) = pgr- [p* (Po -Poo) +Poo]+P^ (vf*) 

^*=P3^- l^* -pco)]+pV^T^V* 

(V*-T)v* = [pg]r- 


dV* 

wn ^ + 


o 

1 

8 

Tp* + 

\py] 

L J 


_L2 _ 


V*2y* 


( 6 / 11 ) 


In equation 6/11, the terms in brackets each appear in front of non-dimen¬ 
sional (unit) vectors. These bracketed terms all have the same dimension, i.e. 
kgm“^ s“^. Multiplying each by (of dimension m^ s^ kg“^), we obtain a 
purely non-dimensional equation: 


'[L 

V 


dV* 

'W 


-I- 


[ 1 ] 


'gL 


o 

1 

8 

_1 


p 

y\ 

9 - 

pV^ 

V p -1- 

pVL_ 
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Equation 6/12 does not bring any information on top of the original incom¬ 
pressible Navier-Stokes equation (eq. 6/5). Instead, it merely separates it into 
two distinct parts. The first is a scalar field (purely numbers, and noted in 
brackets), which indicates the magnitude of the acceleration field. The other 
is a unit vector field (a field of oscillating vectors, all of length one, and noted 
with stars), which represents the geometry (direction) of the acceleration 
field. In this form, we can more easily observe and quantify the weight of the 
different terms relative to one another. This is the role of the flow parameters. 


6.3.2 The flow parameters of Navier-Stokes 


From here on, we convene to call the terms written in brackets in eq. 6/12 
flow parameters, and label them with the following definitions: 

[St] = ^ (6/13) 

[Eu] s (6/14) 

p 

[Fr] = (6/15) 

yjgL 

[Re] = P- - (6/16) 

P 


This allows us to re-write eq. 6/12 as: 


[St] — + [1] y* • v*y* 
ot* 


1 

F? 


r - [Eu] + 


1 




And finally, we arrive at the simple, elegant and formidable non-dimensional 
incompressible Navier-Stokes equation expressed with flow parameters: 
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[St] 


dV* ^ ^ 

— + [1] r • 

dt* 


1 


g 


[Eu] Tp* + 


1 




(6/17) 


Equation 6/17 is an incredibly potent tool in the study of fluid mechanics, for 

two reasons: 

1. It allows us to quantify the relative weight of the different terms in a 
given flow. 

In this way, it serves as a compass, helping us determine which terms 
can safely be neglected in our attempts to find a particular flow solution, 
merely by quantifying the four parameters [St], [Eu], [Fr], and [Re], 

2. It allows us to obtain dynamic similarity between two flows at dif¬ 
ferent scales (fig. 6.1). 

In order that a fluid flow be representative of another flow (for example 
in order to simulate airflow around an aircraft with a model in a wind 
tunnel), it must generate the same vector field V*. In order to do this, 
we must generate an incoming flow with the same four parameters [St], 
[Eu], [Fr], and [Re]. 

Let us therefore take the time to explore the signification of these four 

parameters. 

f L 

The Strouhal number [St] = ^ (eq. 6/13) quantifies the influence of un¬ 
steadiness effects over the acceleration field. It does this by evaluating 
the importance of the representative frequency f at which the flow 
pattern repeats itself. Very high frequencies denote highly unsteady 
flows. When the frequency is very low, [St] is very small; the flow is 
then quasi-steady and it can be solved at a given moment in time as if 
it was entirely steady. 

The Euler number [Eu] = (eq. 6/14) quantifies the influence of the 
pressure gradient over the acceleration field. It does this by comparing 
the largest relative pressure po - poo to the flow of momentum in 
the flow field. The greater [Eu] is, and the more the changes in the 
velocity field V are likely to be caused by pressure gradients rather 
than viscosity, convection or unsteadiness. 

The Froude number [Fr] = (eq. 6/15) quantifies the relative impor- 

yg L 

tance of gravity effects. In practice, gravity effects only play an impor¬ 
tant role in free surface flows, such as waves on the surface of a water 
reservoir. In most other cases, gravity contributes only to a hydrostatic 
effect, which has little influence over the velocity field. 

The Reynolds number [Re] = (eq. 6/16, also eqs. 0/27 p.21 & 5/25 
p.ll4) quantifies the influence of viscosity over the acceleration field. 
It does this by comparing the magnitude of inertial effects (pVL) to 
viscous effects (p). When [Re] is very large, viscosity plays a negligible 
role and the velocity field is mostly dictated by the inertia of the fluid. 
We return to the significance of the Reynolds number in §6.5 below. 
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The Mach number [Ma] ^Vjc (eq. 0/10 p.l5) compares the flow speed V 
with that of the molecules within the fluid particles (the speed of 
sound c). [Ma] does not appear in equation 6/17 because from the start 
in chapter 4, we chose to restrict ourselves to non-compressible flows. 
If we hadn’t, an additional term would exist on the right-hand side, 
-l-|pV which is null when the divergent of velocity V • T is zero 

(which is always so at low to moderate flow speeds). 

We shall explore the significance of the Mach number in chapter 9. 

These five flow parameters should be thought of scalar fields within the 
studied flow domain: there is one distinct Reynolds number, one Mach 
number etc. for each point in space and time. Nevertheless, when describing 
fluid flows, the custom is to choose for each parameter a single representative 
value for the whole flow. For example, when describing pipe flow, it is 
customary to quantify a representative Reynolds number [Re]^) based on 
the average flow velocity Tav. and the pipe diameter D (as we have seen 
with eq. 5/26 p. 116), while the representative Reynolds number for flow 
over an aircraft wing is often based on the free-stream velocity Too and the 
wing chord length. Similarly, the flight Mach number [MaJoo displayed on 
an aircraft cockpit instrument is computed using the relative free-stream 
air speed Too and the free-stream speed of sound Coo, rather than particular 
values measured closer to the aircraft. 

6.3.3 Scaling flows 

Although the above mathematical treatment may be intimidating, it leads 
us to a rather simple and useful conclusion, which can be summarized as 
follows. 

When studying one particular fluid flow around or within an object, it may be 
useful to construct a model that allows for less expensive or more comfortable 
measurements (fig. 6.1). In that case, we need to make sure that the fluid flow 
on the model is representative of the real case. This can be achieved only if 
the relevant flow parameters have the same value in both cases. If this can 
be achieved —although in practice this can rarely be achieved exactly!— then 
the force and power coefficients also keep the same value. The scale effects 
are then fully-controlled and the model measurements are easily extrapolated 
to the real case. 



Figure 6.1 - In order for the flow around a wind tunnel model to be representative of 
the flow around the real-size aircraft (here, a 48 m-wide Lockheed C-141 Star lifter), 
dynamic similarity must be obtained. The value of all flow parameters must be kept 
identical. This is not always feasible in practice. 

Wind tunnel photo by NASA (public domain) 
Full-size aircraft photo CC-by by Peter Long 
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6.4 Flow parameters as force ratios 


This topic is well covered in Massey [5] 

Instead of the above mathematical approach, the concept of flow parameter 
can be approached by comparing forces in fluid flows. 

Fundamentally, understanding the movement of fluids requires applying 
Newton’s second law of motion: the sum of forces which act upon a fluid 
particle is equal to its mass times its acceleration. We have done this in 
an aggregated manner with integral analysis (in chapter 3, eq. 3/8 p.65), 
and then in a precise and all-encompassing way with differential analysis 
(in chapter 4, eq. 4/37 p.94). With the latter method, we obtain complex 
mathematics suitable for numerical implementation, but it remains difficult 
to obtain rapidly a quantitative measure for what is happening in any given 
flow. 

In order to obtain this, an engineer or scientist can use force ratios. This 
involves comparing the magnitude of a type of force (pressure, viscous, 
gravity) either with another type of force, or with the mass-times-acceleration 
which a fluid particle is subjected to as it travels. We are not interested in 
the absolute value of the resulting ratios, but rather, in having a measure of 
the parameters that influence them, and being able to compare them across 
experiments. 


6.4.1 Acceleration vs. viscous forces: 
the Reynolds number 


The net sum of forces acting on a particle is equal to its mass times its 
acceleration. If a representative length for the particle is L, the particle 
mass grows proportionally to the product of its density p and its volume L^. 
Meanwhile, its acceleration relates how much its velocity V will change over 
a time interval At: it may be expressed as a ratio AVI At. In turn, the time 
interval At may be expressed as the representative length L divided by the 
velocity V, so that the acceleration may be represented as proportional to the 
ratio V AT/L. Thus we obtain: 


|net force I = |mass x acceleration! ~ pL^ 


VAV 

L 


Ifnetl - pL^VAV 


We now observe the viscous force acting on a particle: it is proportional to 
the the shear effort and a representative acting surface L^. The shear can 
be modeled as proportional to the viscosity p and the rate of strain, which 
will grow proportionally to AT/L. We thus obtain a crude measure for the 
magnitude of the shear force: 


viscous force j 


\Fy 



pAVL 


The magnitude of the viscous force can now be compared to the net force: 

jnet force! pL^VAV _ pVL _ 

!viscous force! pAVL 




(6/18) 







and we recognize the ratio as the Reynolds number (6/16). We thus see that 
the Reynolds number can be interpreted as the inverse of the influence of 
viscosity. The larger [Re] is, and the smaller the influence of the viscous 
forces will be on the trajectory of fluid particles. 

6.4.2 Acceleration vs. gravity force: the Fronde number 

The weight of a fluid particle is equal to its mass, which grows with pV", 
multiplied by gravity g-. 


I weight force I = |Fwl ~ P^9 


The magnitude of this force can now be compared to the net force: 


|netforce| pL^V^ 



(6/19) 


I weight force! pL'^g Lg 


and here we recognize the square of the Froude number (6/15). We thus see 
that the Froude number can be interpreted as the inverse of the influence of 
weight on the flow. The larger [Fr] is, and the smaller the influence of gravity 
will be on the trajectory of fluid particles. 

6.4.3 Acceleration vs. elastic forces: the Mach number 

In some flows called compressible flows the fluid can perform work on itself, 
and and the fluid particles then store and retrieve energy in the form of 
changes in their own volume. In such cases, fluid particles are subject to 
an elastic force in addition to the other forces. We can model the pressure 
resulting from this force as proportional to the bulk modulus of elasticity K 
of the fluid (formally defined asK = p dpldp); the elastic force can therefore 
be modeled as proportional to KL^: 


ry 

I elasticity force! = l-Pelastkl ~ KL 


The magnitude of this force can now be compared to the net force: 


!net force! pL^V^ pV^ 


!elasticity force! KL^ K 


This ratio is known as the Cauchy number; it is not immediately useful 
because the value of FC in a given fluid varies considerably not only according 
to temperature, but also according to the type of compression undergone by 
the fluid: for example, it grows strongly during brutal compressions. 

During isentropic compressions and expansions (isentropic meaning that the 
process is fully reversible, i.e. without losses to friction, and adiabatic, i.e. 
without heat transfer), we will show in chapter 9 (with eq. 9/7 p.l95) that 
the bulk modulus of elasticity is proportional to the square of the speed of 
sound c: 



( 6 / 20 ) 
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The Cauchy number calibrated for isentropic evolutions is then 


I net force I 

I elasticity force |i-eversible 


pV'^ 


^ = [Ma] 

C2 


2 


( 6 / 21 ) 


and here we recognize the square of the Mach number (0/10). We thus see 
that the Mach number can be interpreted as the influence of elasticity on the 
flow. The larger [Ma] is, and the smaller the influence of elastic forces will be 
on the trajectory of fluid particles. 


6.4.4 Other force ratios 

The same method can be applied to reach the definitions for the Strouhal and 
Euler numbers given in §6.3. Other numbers can also be used which relate 
forces that we have ignored in our study of fluid mechanics. For example, the 
relative importance of surface tension forces or of electromagnetic forces are 
quantified using similarly-constructed flow parameters. 

In some applications featuring rotative motion, such as flows in centrifugal 
pumps or planetary-scale atmospheric weather, it may be convenient to 
apply Newton’s second law in a rotating reference frame. This results in the 
appearance of new reference-frame forces, such as the Coriolis or centrifugal 
forces; their influence can then be studied using additional flow parameters. 

In none of those cases can flow parameters give enough information to predict 
solutions. They do, however, provide quantitative data to indicate which 
forces are relevant in which places: this not only helps us understand the 
mechanisms at work, but also distinguish the negligible from the influential, 
a key characteristic of efficient scientific and engineering work. 


6.5 The Reynolds number in practice 


Among the five non-dimensional parameters described above, the Reynolds 
number [Re] is by far the most relevant in the study of most fluid flows, and 
it deserves a few additional remarks. As we have seen, the Reynolds number 
is a measure of how little effect the viscosity has on the time-change of the 
velocity vector field: 

• With low [Re], the viscosity p plays an overwhelmingly large role, and 
the velocity of fluid particles is largely determined by that of their own 
neighbors; 

• With high [Re], the momentum pV of the fluid particles plays a more 
important role than the viscosity p, and the inertia of fluid particles 
affects their trajectory much more than the velocity of their neighbors. 



Video: half-century-old, but 

timeless didactic exploration 
how the dynamics of fluids 
change with Reynolds number, 
with accompanying notes by 
Taylor[2] 

by the National Committee for Fluid 
Mechanics Films (ncfmf, 1967[12]) (styl) 
https://youtu.be/51-6QCJTAjU 


In turn, this gives the Reynolds number a new role in characterization of 
flows: it can be thought of as the likeliness of the flow being turbulent 
over the length L. Indeed, from a kinematic point of view, viscous effects are 
highly stabilizing: they tend to harmonize the velocity field and smooth out 
disturbances. On the contrary, when these effects are overwhelmed by inertial 
effects, velocity non-uniformities have much less opportunity to dissipate, 
and at the slightest disturbance the flow will become and remain turbulent 
(fig. 6.2). This is the reason why the quantification of a representative 
Reynolds number is often the first step taken in the study of a fluid flow. 
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Figure 6.2 - A viscous opaque fluid is dropped into a clearer receiving static fluid 
with identical viscosity. The image shows four different experiments photographed 
after the same amount of time has elapsed. The viscosity is decreased from left to 
right, yielding Reynolds numbers of 0,05, 10, 200 and 3 000 respectively. 

As described in eq. 6/18, low Reynolds number indicate that viscous effects dominate 
the change in time of the flow field. As the Reynolds number increases, the nature 
of the velocity field changes until it becomes clearly turbulent. 

By the National Committee for Fluid Mechanics Films (ncfmf, 1967[12]), with accompanying notes by Taylor[2] 
Image ©1961-1969 Educational Development Center, Inc., reproduced under Fair Use doctrine 
A screen capture from film Low Reynolds Number Flow at http://web.mitedu/hml/ncfmfhtml 





Exercise sheet 6 (Scale effects) 

last edited April 2, 2018 


Except otherwise indicated, we assume that fluids are Newtonian, and that: 

Pwater = 1 000 kg m“^ patm. = 1 bar; Patm. = 1,225 kgm“^ Tatm. = ll,3°C;Patm. = 1,5 • lO'^Nsm”^; 
g = 9,81 ms“^. Air is modeled as a perfect gas (i?air = 287 JK“' kg“^; Yaa = 1,4; Cpair = 1 005Jkg“' K“'). 


Non-dimensional incompressible Navier-Stokes equation: 


^1/^* 1 1 
[St] ^ + [1] 7* • TV* = —r ^ - [Eu] Tp* + — T^V* 
dt* [Fr]^ [Re] 

in which [St] = [Eu] = [Fr] = and [Re] = 

P ^ aJg L P 


The force coefficient Cf and power coefficient Cp are defined as: 
F _ W 


Cf 


Cp 


IpSV^ IpSV^ 

The speed of sound c in air is modeled as: 

c = Vr^ 


(6/17) 


(6/4) 


( 0 / 11 ) 


Figure 6.3 quantifies the viscosity of various fluids as a function of temperature. 


o 



XXX 
jUi/s . ^ Ajisoosta 9fn[osqv 


Figure 6.3 - Viscosity of various fluids at a pressure of 1 bar (in practice viscosity is almost 
independent of pressure). 

Figure repeated from Gg. 5.11 p.ll9;© White 2008 [13j 
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6.1 


Scaling with the Reynolds number 


CC-O o.c. 


1. Give the definition of the Reynolds number, indicating the si units for each term. 

2. What is the consequence on the velocity field of having a low Reynolds number? 
A high Reynolds number? 

3. Give an example of a high Reynolds number flow, and of a low Reynolds num¬ 
ber flow. 

The standard golf ball has a diameter of 42,67 mm and a mass of 45,93 g. A typical 

maximum velocity for such balls is 200 kmh“h 

A student wishes to investigate the flow over a golf ball using a model in a wind tunnel. 

S/he prints an enlarged model with a diameter of 50 cm. 

4. If the atmospheric conditions are identical, what flow velocity needs to be generated 
during the experiment in order to reproduce the flow patterns around the real ball? 

5. Would the Mach number for the real ball then be reproduced? 

6. Would this velocity need to be adjusted if the experiment was run on a very hot 
summer day? 

7. If the model was made of the same material as the real ball, how heavy would 
it be? 


6.2 Fluid mechanics of a giant airliner 

non-examinable, CC-O o.c. 

In this exercise, we consider that that viscosity is entirely independent from pressure 
and density (this is a reasonable first approximation). 

An Airbus A380 airliner has an 80 m-wingspan and is designed to cruise at [Ma] = 0,85 
where the air temperature and pressure are -40 °C and 0,25 bar. 

A group of students wishes to study the flow field around the aircraft, with a wind tunnel 
whose test section has a diameter of 1 m (in which, obviously, the model has to fit!). 

1. If the temperature inside the tunnel is maintained at Ttunnel = 20 °C, propose a 
combination of wind tunnel velocity and pressure which would enable the team to 
adequately model the effects of viscosity. 

2. Which flow conditions would be required in a water tunnel with the same dimen¬ 
sions? 

Shocked by their answers to the above questions, the group of students decides instead 
to study the effect of compressibility (while accepting that viscous effects may not be 
adequately modeled). 

3. If the maximum velocity attainable in the wind tunnel is 80ms“\ which air 
temperature is required for compressibility effects to be modeled? 
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6.3 Scale effects on a dragonfly 


A dragonfly (sketched in fig. 6.4) has a 10 cm wingspan, a mass of 80 mg, and cruises 
at 4ms“^, beating its four wings 20 times per second. 


10 cm 


Figure 6.4 - Plan view of a dragonfly 



Figure by o.c., A. Plank, Drury Dru, Westwood, J. O. (CC-0) 


You are tasked with the investigation of the flight performance of the dragonfly, and 
have access to a wind tunnel with a test section of diameter 1 m. 

1. Which model size and flow velocity would you use for this experiment? 

2. How many wing beats per second would then be required on the model? 

3. What would be the lift developed by the model during the experiment? 

4. How much mechanical power would the model require, compared to the real 
dragonfly? 


6.4 Formula One testing 

You are leading the Aerodynamics team in a successful Formula One racing team. Your 
team races a car that is 5,1 m long, 1,8 m wide, has a 610 kW power plant, and a mass 
of 750kg. The car can reach a top speed of 310kmh“\ at which speed the drag force 
is 7,1 kN. 

In order to test different aerodynamic configurations for the car, your team invests in a 
50-million euro wind tunnel in which you will run tests on a model. You currently have 
the choice between two possible model sizes: a 60 % model and a (smaller) 50 % model. 

1. What would be the frontal area of each model, in proportion to the frontal area of 
the real car? 

2. What would be the volume of each model, compared to the volume of the real car? 

3. How much less weight would the 50 % model have than the 60 % model? 

In the end, the regulations for the Formula One racing change, and your team is forced 
to opt for a 50 % car model. 


143 









4. If the ambient atmospheric conditions cannot be changed, which flow speed in 
the wind tunnel is required, so that the air flow around the real car is reproduced 
around the model? 

In practice the regulations change once again, and the maximum speed authorized in 
wind tunnel tests is now 50 m s“h 

5. Your team considers modifying the air temperature to compensate for the limit in 
the air speed. If the temperature in the tunnel can be controlled between -10 °C 
and 40 °C, but the pressure remains atmospheric, what is the maximum race-track 
speed that can be reproduced in the wind tunnel? 

6. In that case, by which factor should the model drag force measurements be multi¬ 
plied in order to correspond to the real car? 

7. [non-examinable question] The wind tunnel has a 2 m by 2 m square test cross- 
section. What is the power required to bring atmospheric air (Tatm. = 15 °C, V"atm. = 
0ms“^) to the desired speed and temperature in the test section? 

NB: this exercise is inspired by an informative and entertaining video by the SauberFl team about 
their wind tunnel testing, which the reader is encouraged to watch at 
https://www.youtube.com/watch?v=KC0E0wU6inU. 
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Answers 


6.1 1 ) See equation 6/16 and subsequent comments; 

2 ) See §6.5, in practice this can be summarized in two or three sentences; 

3) High-Re: air flow around an airliner in cruise; Low-Re: air flow around a dust 

particle falling to the ground. 4) Match [Re]: V 2 = 4,74ms“^ 

5) No, but it is very low anyway (no compressibility effects to be reproduced) 

6 ) Yes! (discuss in class) 7) m 2 = 73,9 kg 

6.2 1 ) With a half-aircraft model of half-length y = 80 cm, at identical speed, we 

would need an air density p 2 = 24,9 kg m“^. At ambient temperature, this requires 
a pressure p 2 = 20,9 bar! 2) In water, the density cannot be rea¬ 

sonably controlled, and we need a velocity V 3 = 418 m s“M 

3) In air, at V 4 = 80ms“\ Mach number can be reproduced at T 4 = -251 °C (al¬ 
though the Reynolds number is off). Wind tunnels used to investigate compressible 
flow around aircraft have very powerful coolers. 

6.3 1 ) With e.g. a model of span 60 cm, match the Reynolds number: V 2 = 0,67 m s“^ 

2 ) Match the Strouhal number: /2 = 0,56 Hz. Mach, Froude and Euler numbers will 
have no effect here; 3) and 4) are left as a surprise for the 

student. 

6.4 3) The third model would have 42 % less mass than the second; 

4) Maintaining [Re] requires Vs = 172,2 ms“^ (fast!); 

5) Now the fastest track speed that can be studied is 27 m s“^ (assuming identical 

viscosity everywhere); 6 ) Multiply force measurements by 1,74 

to scale up to reality; 7) One needs powerful coolers! 
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These lecture notes are based on textbooks by White [13], Qengel & al.[16], and Munson & al.[18]. 


7.1 Motivation 


In this chapter, we focus on fluid flow close to solid walls. In these regions, 
viscous effects dominate the dynamics of fluids. This study should allow us 
to answer two questions: 

• How can we quantify shear-induced friction on solid walls? 

• How can we describe and predict flow separation? 


7.2 Describing the boundary layer 


7.2.1 Concept 

At the very beginning of the 20* century, Ludwig Prandtl observed that for 
most ordinary fluid flows, viscous effects played almost no role outside of a 
very small layer of fluid along solid surfaces. In this area, shear between the 
zero-velocity solid wall and the outer flow dominates the flow structure. He 
named this zone the boundary layer. 

We indeed observe that around any solid object within a fluid flow, there 
exists a thin zone which is significantly slowed down because of the object’s 
presence. This deceleration can be visualized by measuring the velocity 
profile (fig. 7.1). 

The boundary layer is a concept, a thin invisible layer whose upper limit 
(termed S, as we will see in §7.2.3) is defined as the distance where the fluid 
velocity in direction parallel to the wall is 99 % of the outer (undisturbed) 
flow velocity. 



Video: pre-lecture briefing for 
this chapter 

by o.c. (CC-by) 
https://youtu.be/kUlm5ojOnUQ 


147 








Figure 7.1 - A typical velocity profile in a boundary layer. Only the horizontal 
component of velocity {Vx = u) is represented. 

Figure CC-by Olivier Cleynen 


Upon investigation, we observe that the boundary layer thickness depends 
strongly on the main flow characteristics. In particular, it decreases when 
speed increases or when viscosity is decreased (fig. 7.2). 

As we travel downstream along a boundary layer, we observe experimentally 
that the flow regime is always laminar at first. Then, at some distance 
downstream which we name transition point, the boundary layer becomes 
turbulent. The flow-wise position of the transition point depends on the flow 
properties and is somewhat predictable. Further downstream, the boundary 
layer becomes fully turbulent. It has larger thickness than in the laminar 
regime, and grows at a faster rate. Like all turbulent flows, it then features 
strong energy dissipation and its analytical description becomes much more 
difficult. 

The flow within the boundary layer, and the main external flow (outside of it) 
affect one another, but may be very different in nature. A laminar boundary 
layer may exist within a turbulent main flow; while turbulent boundary 
layers are commonplace in laminar flows. 


U 



U 



Figure 7.2 - The thickness of the boundary layer depends strongly on the faraway 
incoming flow velocity Uco- 
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Figure CC-by-sa Olivier Cleynen & Commons User.F 1 an k e r 



































7.2.2 Why do we study the boundary layer? 


Expending our energy on solving such a minuscule area of the flow may 
seem counter-productive, yet three great stakes are at play here: 

• First, a good description allows us to avoid having to solve the 
Navier-Stokes equations in the whole flow. 

Indeed, outside of the boundary layer and of the wake areas, viscous 
effects can be safely neglected. Fluid flow can then be described with 
= -Vp (the Euler equation, which we will introduce as eq. 8/3 
in the next chapter) with acceptably low inaccuracy. As we shall 
soon see, this allows us to find many interesting solutions, all within 
reach of human comprehension and easily obtained with computers. 
Unfortunately, they cannot account for shear on walls. 

Thus, solving the flow within the boundary layer, whether analytically 
or experimentally, allows us to solve the rest of the flow in a simplified 
manner (fig. 7.3). 



Figure 7.3 - Fluid flow around a wing profile. When analyzing the flow, whether 
analytically or within a computational fluid dynamics (cfd) simulation, the flow 
domain is frequently split into three distinct areas. In the boundary layer (B), fluid 
flow is dominated by viscosity. Outside of the boundary layer (A), viscous effects 
are very small, and the flow can be approximately solved using the Euler equation. 
Lastly, in the turbulent wake (C), characterization of the flow is very difficult and 
typically necessitates experimental investigations. 

Figure CC-by-sa Olivier Cleynen 


• Secondly, the boundary layer is the key to quantifying friction. A 
good resolution of the boundary layer allows us to precisely quantify 
the shear forces generated by a fluid on an object. 

• Finally, a good understanding of the mechanisms at hand within the 
boundary layer allows us to predict flow separation, which is the di¬ 
vergence of streamlines relative to the object. Control of the boundary 
layer is key to ensuring that a flow will follow a desired trajectory! 

7.2.3 Characterization of the boundary layer 

Three different parameters are typically used to quantify how thick a bound¬ 
ary layer is at any given position. 


The first is the thickness 5, 


<5 = y\u=0,99U (7/1) 

which is, as we have seen above, equal to the distance away from the wall 
where the speed u is 99 % of U. 

The second is the displacement thickness 5*, which corresponds to the 
vertical distance by which streamlines outside of the boundary layer have 
been displaced. This vertical “shifting” of the flow occurs because the inner 
fluid is slowed down near the wall, creating some blockage for the outer flow, 
which then proceeds to avoid it partially by deviating outwards (fig. 7.4). 

Integral analysis performed on a control volume enclosing the boundary 
layer (as for example in exercise 3.7 from chapter 3) allows us to quantify the 
displacement thickness as: 

In practice, this integral can be calculated on a finite interval (instead of using 
the oo limit), as long as the upper limit exceeds the boundary layer thickness. 

The third and last parameter is the momentum thickness 9, which is equal 
to the thickness of a corresponding layer of fluid at velocity U which would 
possess the same amount of momentum as the boundary layer. The momen¬ 
tum thickness can be thought of as the thickness of the fluid that would need 
to be entirely stopped (for example by pumping it outside of the main flow) 
in order to generate the same drag as the local boundary layer. 

A review of the concepts of chapter 3 (e.g. in exercises 3.6 to 3.8) allows us to 
to quantify the momentum thickness as: 

with the same remark regarding the upper limit. The momentum thickness, in 
particular when compared to the displacement thickness, is a key parameter 
used in prediction models for boundary layer separation. 

Once these three thicknesses have been quantified, we are generally looking 
for a quantification of the shear term r^aii- Since we are working with the 
hypothesis that the fluid is Newtonian, we merely have to know U(y) to 



Figure 7.4 - Thickness S and displacement thickness <5* of a boundary layer. It is 
important to understand that the boundary layer is not a closed zone: streamlines 
(drawn blue) penetrate it and the vertical velocity v, although very small compared 
to M, is not zero. 

Figure CC-by-sa Olivier Cleynen & Commons User:F lanke r 
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quantify shear according to equation 2/16 (p.51): 


hvall yx ~ ('7/4) 

this expression is a function of x, and typically r^aii decreases with increasing 
distance inside a laminar boundary layer. 

The wall shear exerted by the boundary layer is typically non-dimensionalized 
with the shear coefficient Cf, 


'^wall 

where U is the outer-layer (free-stream) velocity. 


(7/5) 


The shear coefficient, just like the shear, remains a function of the flow-wise 
distance x. 

Once these parameters are known, the wall shear force is simply obtained by 
integrating r over the entire surface S: 


shestryx 


^fs" 


hvall yx dx dz 


(7/6) 


7.3 The laminar boundary layer 


7.3.1 Governing equations 


What is happening inside a laminar, steady boundary layer? We begin by 
writing out the Navier-Stokes for incompressible isothermal flow in two 
Cartesian coordinates (eqs. 4/41 & 4/42 p.95): 


P 

P 


du 

dt 

dv 

dt 


+ u- 


+ u- 


du du 

dp 

d^u d^u 

-h V - 

dx dy 


(dx)2 ' (dyY 

dv dv 

dp 

d^v d^v 

-h V - 

dx dy 

= pgy-Yy+p 

1 

CM 

+ 

CM 

52 


(7/7) 

(7/8) 


Building from these two equations, we are going to add three simplifica¬ 
tions, which are hypotheses based on experimental observation of fluid flow 
in boundary layers: 


1. Gravity plays a negligible role; 


2. The component of velocity perpendicular to the wall (in our conven¬ 
tion, v) is very small {v <^u). 

Thus, its stream-wise spatial variations can also be neglected: ~ 0 


and 




{dxp 


0. The same goes for the derivatives in the y-direction: 

0 . 


3. The component of velocity parallel to the wall (in our convention, u) 

varies much more strongly in the y-direction than in the x-direction: 

d'^u d'^u 

















With all of these simplifications, equation 7/8 shrinks down to 


dp 

dy 


« 0 


(7/9) 


which tells us that pressure is a function of x only (|| = 

We now turn to equation 7/7, first to obtain an expression for pressure by 
applying it outside of the boundary layer where u = U: 


dp 

dx 



(7/10) 


and secondly to obtain an expression for the velocity profile: 


du du 


\ dp p d^u 
pdx^ p {dyf 
^ dU ^ p d^u 
cbr p {dyY 


(7/11) 


Thus, the velocity field V = {u-,v) = f{x,y) in a steady laminar boundary 
layer is driven by the two following equations, answering to Newton’s second 
law and continuity: 


du du dU p d^u 

du dv 
dx dy 


(7/12) 

(7/13) 


The main unknown in this system is the longitudinal speed profile across the 
layer, U(^x,y)- Unfortunately, over a century after it has been written, we still 
have not found an analytical solution to it. 


7.3.2 Blasius’ solution 

Heinrich Blasius undertook a PhD thesis under the guidance of Prandtl, in 
which he focused on the characterization of laminar boundary layers. As 
part of his work, he showed that the geometry of the velocity profile (i.e. 
the velocity distribution) within such a layer is always the same, and that 
regardless of the flow velocity or the position, u can be simply expressed as a 
function of non-dimensionalized distance away from the wall termed y: 


y^y 



(7/14) 


Blasius was able to show that u is a function such that ^ , with f'" + 

\ff'' = 0 ~ unfortunately, no known analytical solution to this equation is 
known. However, it has now long been possible to obtain numerical values 
for f at selected positions y. Those are plotted in fig. 7.5. 

Based on this work, it can be shown that for a laminar boundary layer flowing 
along a smooth wall, the four parameters about which we are interested are 
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Figure 7.5 - The velocity profile obtained by Blasius (an exact solution to the Navier- 
Stokes equations simplified with laminar boundary-layer hypothesis). 

Figure CC-0 o.c. 


solely function of the distance-based Reynolds number [Re]^: 


5 _ 

4,91 

(7/14a) 
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5* 

1,72 

(7/14b) 
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0,664 

(7/14d) 
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7.3.3 Pohlhausen’s model 

Contrary to Blasius who was searching for an exact solution, Ernst Pohlhausen 
made an attempt at modeling the velocity profile in the boundary layer with 
a simple equation. He imagined an equation of the form: 

^ = g^Y) =aY + bY^ + cY^ + dY^ (7/15) 

in which Y = y/5 represents the distance away from the wall, and he then 
set to find values for the four coefficients a, b, c and d that would match the 
conditions we set in eqs. (7/12) and (7/13). 

In order to obtain these coefficients, Blasius simply calibrated his model using 
known boundary conditions at the edges of the boundary layer: 

• for y = 0 (meaning Y = 0), we have both u = 0 and z; = 0. Equation 7/12 
becomes: 


P „ 


(7/16) 
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for y = 5 (meaning 7=1) we have u = U, du/dy = 0 and d^ul{dy)^ = 0 . 
These values allow us to evaluate three properties of function g ; 


u = U = Ug(i) 

(7/17) 

du U , 

(7/18) 

U .. 

{dyY , 52 ^( 1 ) 

(7/19) 


This system of four equations allowed Pohlhausen to evaluate the four terms 
a, b, c and d. He introduced the variable A, a non-dimensional measure of 
the outer velocity flow-wise gradient: 

cpp dU 

A = (5-— (7/20) 

fi ax 

and he then showed that his velocity model in the smooth-surface laminar 
boundary layer, equation 7/15, had finally become: 

- = !-{! +Y){l-Y)^+A-{1-Y^) (7/21) 

U 6 

Pohlhausen’s model is remarkably close to the exact solution described by 
Blasius (fig. 7.6). We thus have an approximate model for U(x,t/) within the 
boundary layer, as a function of 17(x). 



Figure 7.6 - Velocity profiles in the Pohlhausen and Blasius boundary layer profiles. 

Figure © from Richecoeeur 2012 [17] 


It can be seen now that neither Blasius nor Pohlhausen’s models are fully 
satisfactory. The first is drawn from a handful of realistic, known boundary 
conditions, but we fail to find an analytical solution to match them. The 
second is drawn from the hypothesis that the analytical solution is simple, 
but it is not very accurate. And we have not even considered yet the case for 
turbulent flow! 

Here, it is the process which matters to us. Both of those methods are 
frequently-attempted, sensible approaches to solving problems in modern 
fluid dynamics research, where, just as here, one often has to contend with 
approximate solutions. 
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7.4 Transition 


After it has traveled a certain length along the wall, the boundary layer 
becomes very unstable and it transits rapidly from a laminar to a turbulent 
regime (fig. 7.7). We have already briefly described the characteristics of a 
turbulent flow in chapter 5; the same apply to turbulence within the boundary 
layer. Again, we stress that the boundary layer may be turbulent in a globally 
laminar flow (it may conversely be laminar in a globally turbulent flow): that 
is commonly the case around aircraft in flight, for example. 


laminar turbulent 



Figure 7.7 - Transition of a boundary layer from laminar to turbulent regime. 

Figure CC-by Olivier Cleynen 


We observe that the distance Xtransition at which the boundary layer changes 
regime is reduced when the velocity is increased, or when the viscosity is 
decreased. In practice this distance depends on the distance-based Reynolds 
number [Re];c = The most commonly accepted prediction for the transi¬ 
tion position is: 


[ReJx transition « 5 • 10^ (7/22) 

Transition can be generated earlier if the surface roughness is increased, or 
if obstacles (e.g. turbulators, vortex generators, trip wires) are positioned 
within the boundary layer. Conversely, a very smooth surface and a very 
steady, uniform incoming flow will result in delayed transition. 


7.5 The turbulent boundary layer 


The extensive description of turbulent flows remains an unsolved problem. 
As we have seen in chapter 5 when studying duct flows, by contrast with 
laminar counterparts turbulent flows result in 

• increased mass, energy and momentum exchange; 

• increased losses to friction; 


• apparently chaotic internal movements. 


Instead of resolving the entire time-dependent flow in the boundary layer, we 
satisfy ourselves with describing the average component of the longitudinal 
speed, u. A widely-accepted velocity model is: 


u 

U 




(7/23) 
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for turbulent boundary layer flow over a smooth surface. 


This profile has a much flatter geometry near the wall than its laminar 
counterpart (fig. 7.8). 

In the same way that we have worked with the laminar boundary layer 
profiles, we can derive models for our characteristics of interest from this 
velocity profile: 


(5 ^ 

0,16 

1 

(7/23a) 
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(7/23b) 
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[Re]| 
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(7/23d) 


[Re]l 



Figure 7.8 - Comparison of laminar and turbulent boundary layer profiles, both 
scaled to the same size. It is important to remember that in practice turbulent 
boundary layers are much thicker than laminar ones, and that u is only a time- 
average velocity. 


Figure CC-0 o.c. 










7.6 Separation 


7.6.1 Principle 

Under certain conditions, fluid flow separates from the wall. The boundary 
layer then disintegrates and we observe the appearance of a turbulent wake 
near the wall. Separation is often an undesirable phenomenon in fluid 
mechanics: it may be thought of as the point where we fail to impart a 
desired trajectory to the fluid. 

When the main flow speed U along the wall is varied, we observe that the 
geometry of the boundary layer changes. The greater the longitudinal speed 
gradient (dU / dx > 0), and the flatter the profile becomes. Conversely, when 
the longitudinal speed gradient is negative, the boundary layer velocity 
profile straightens up. When it becomes perfectly vertical at the wall, it 
is such that streamlines separate from the wall: this is called separation 
(fig. 7.9). 

The occurrence of separation can be predicted if we have a robust model 
for the velocity profile inside the boundary layer. For this, we go back to 
fundamentals, stating that at the separation point, the shear effort on 
the surface must be zero: 


_ _ ,du 

hvall at separation ~ b = /i I 7 — , 

\^yly=o 


(7/24) 


At the wall surface (u = 0 and v = 0), equation 7/12 p.l52 becomes: 





dp dU 

-J- = -pU - 

cbc cbc 


(7/25) 


Thus, as we progressively increase the term dp/ dx, the term reaches 
higher (positive) values on the wall surface. Nevertheless, we know that it 
must take a negative value at the exterior boundary of the boundary layer. 



Figure 7.9 - Separation of the boundary layer. The main flow is from left to right, 
and flowing towards a region of increasing pressure. For clarity, the y-scale is greatly 
exaggerated. 


Figure CC-by Olivier Cleynen 





















Video: a football with a rough 
edge on one side will see turbu¬ 
lent boundary layer on that side 
and laminar boundary layer on 
the other. When the laminar 
layer separates before the turbu¬ 
lent one, the main flow around 
the ball becomes asymmetric 
and deviates the ball sideways, 
a phenomenon that can be ex¬ 
ploited to perform impressive 
tricks. 

by freekikerz (styl) 
https://youtu.be/rzRuZNvkMbc 


Therefore, it must change sign somewhere in the boundary. This point where 
changes sign is called inflexion point. 

The existence of the inflection point within the boundary layer tells us that 
the term tends towards a smaller and smaller value at the wall. Given 
enough distance x, it will reach zero value, and the boundary layer will 
separate (fig. 7.10). Therefore, the longitudinal pressure gradient, which in 
practice determines the longitudinal velocity gradient, is the key factor in 
the analytical prediction of separation. 

We shall remember two crucial points regarding the separation of boundary 
layers: 

1. Separation occurs in the presence of a positive pressure gradi¬ 
ent, which is sometimes named adverse pressure gradient 
Separation points along a wall (e.g. a car bodywork, an aircraft wing, 
rooftops, mountains) are always situated in regions where pressure 
increases (positive dp/dx or negative dU/ dx). If pressure remains 
constant, or if it decreases, then the boundary layer cannot separate. 

2. Laminar boundary layers are much more sensitive to separa¬ 
tion than turbulent boundary layers (fig. 7.11). 

A widely-used technique to reduce or delay the occurrence of separa¬ 
tion is to make boundary layers turbulent, using low-height artificial 
obstacles positioned in the flow. By doing so, we increase shear-based 
friction (which increases with turbulence) as a trade-off for better 
resistance to stall. 

Predicting in practice the position of a separation point is difficult, because 
an intimate knowledge of the boundary layer profile and of the (external- 
flow-generated) pressure field are required — and as the flow separates, these 
are no longer independent. Resorting to experimental measurements, in this 
case, is often a wise idea! 



Figure 7.10 - The inflection point within a boundary layer about to separate. 

Figure CC-by Olivier Cleynen 


158 















decreasing 

Reynolds 

number 




V 


Figure 7.11 - The effect of decreasing Reynolds number on flow attachment over an 
airfoil at constant angle of attack, with the transition point highlighted. Laminar 
boundary layers are much more prone to separation than turbulent boundary layers. 

Figure CC-by-sa Olivier Cleynen, based on Barlow & Pope 1999 [9] 


7.6.2 Separation according to Pohlhausen 


In the case of a laminar flow along a smooth surface, Pohlhausen’s model 
allows for interesting results. Indeed, at the separation point, Twall = 0, thus: 


^wall — 0 — p 



t /=0 



0 — — 0 . 
A = -12 


(7/26) 




-12 


(7/27) 


This equation allows us to predict the thickness of the boundary layer at 
the point where separation occurs, and therefore its position, if the velocity 
gradient dU / dx is known. 

We can see here that once we have obtained an extensive description of the 
flow within the boundary layer, we can generate a large amount of useful 
information regarding the general flow around the object. 

Unfortunately, Pohlhausen’s model’s agreement with experimental measure¬ 
ments is not very good, and it can only usefully predict the position of the 
separation point in a few simple cases. Nevertheless, the overall thought 
process, i.e. 

intuition of a simple model 

i 

non-dimensionalization 

i 

correlation with known boundary conditions 
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i 

exploitation of the model to predict flow behavior 


is what should be remembered here, for it is typical of the research process 
in fluid dynamics. 
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Exercise sheet 7 (Flow near walls) 

last edited April 4, 2017 


These lecture notes are based on textbooks by White [13], Qengel & al.[16], and Munson & al.[18]. 


Except otherwise indicated, we assume that fluids are Newtonian, and that: 

Pwater = 1 000 kg m“^ patm. = 1 bar; Patm. = 1,225 kgm“^ Tatm. = ll,3°C;patm. = 1,5 • 10‘^Nsm“2; 
g = 9,81 ms“^. Air is modeled as aperfect gas (i?air = 287JK“^ kg“'; pair = 1,4; Cpair = 1005Jkg“^ K“^). 


In boundary layer flow, we accept that transition occurs at [Re];,. > 5 • 105. 

The wall shear coefficient cp, a function of distance x, 
is defined based on the free-stream flow velocity U: 


Iwall 


(7/5) 


Exact solutions to the laminar boundary layer along a smooth surface yield: 


8 

X 

0 

X 


4,91 

0,664 


8* _ 
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1,72 

VK 

0,664 

VK 


(7/14b) 

(7/14d) 


Solutions to the turbulent boundary layer along a smooth surface yield the following 

time-averaged characteristics: 
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(7/23b) 

(7/23d) 


[Re]; [Re]; 

Figure 7.12 quantifies the viscosity of various fluids as a function of temperature. 
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Figure 7.12 - Viscosity of various fluids at a pressure of 1 bar (in practice viscosity is almost 
independent of pressure). 

Figure repeated from fig. 5.11 p.ll9; © White 2008 [13] 


7.1 Water and air flow 

White [13] E7.2 

A flat plate of length 0,3 m is placed parallel to a uniform flow with speed 0,3 m s“^. How 
thick can the boundary layer become: 

1 . if the fluid is air at 1 bar and 20 °C? 

2 . if the fluid is water at 20 °C? 

7.2 Boundary layer sketches ccooc 

A thin and long horizontal plate is moved horizontally through a stationary fluid. 

1 . Sketch the velocity profile of the fluid: 

• at the leading edge; 

• at a point where the boundary layer is laminar; 

• and at a point further downstream where the boundary layer is turbulent. 

2 . Draw a few streamlines, indicate the boundary layer thickness 5, and the displace¬ 
ment thickness 8*. 

3. Explain shortly (e.g. in 30 words or less) how the transition to turbulent regime 
can be triggered. 

4. Explain shortly (e.g. in 30 words or less) how the transition to turbulent regime 
could instead be delayed. 
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7.3 Shear force due to boundary layer 


M^hite [13] E7.3 


A thin and smooth plate of dimensions 0,5 X 3 m is placed with a zero angle of attack in a 
flow incoming at 1,25 m as shown in fig. 7.13. 



Figure 7.13 - A thin plate positioned parallel to an incoming uniform flow. Two configurations 
are studied in this exercise. 

Figure CC-0 o.c. 


1 . what is the shear force exerted on the plate for each of the two configurations 
shown in fig. 7.13, when the fluid is air of viscosity pair = h5 • 10 “^ Pa s? 

2 . What are the shear forces when the fluid is water of viscosity //water = 1 • 10”^ Pa s? 

3. [difficult question] How would these shear efforts evolve if the plate was tilted with 
an angle of 20 ° relative to the flow? 


7.4 Wright Flyer I 

C_^C^~0 o.c. 

The Wright Flyer I, the first airplane capable of sustained controlled flight (1903), was a 
biplane with a 12 m wingspan (fig. 7.14). It had two wings of chord length 1,98 m stacked 
one on top of the other. The wing profile was extremely thin and it could only fly at very 
low angles of attack. Its flight speed was approximately 40kmh“h 



Figure 7.14 - The Wright Flyer I, first modern airplane in history. Built with meticulous care and 
impeccable engineering methodology by two bicycle makers, it made history in December 1903. 

Photo by Orville Wright, 1908 (public domain) 


1 . If the flow over the wings can be treated as if they were flat plates, what is the 
power necessary to compensate the shear exerted by the airflow on the wings 
during flight? 

2. Which other forms of drag would also be found on the aircraft? (give a brief 
answer, e.g. in 30 words or less) 
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7.5 Power lost to shear on an airliner fuselage 

C_^C^~0 o.c. 

An Airbus A340-600 is cruising at [Ma] = 0,82 at an altitude of 10 000 m (where the air 
has viscosity 1,457 • 10“^Nsm“^, temperature 220 K, density 0,4kgm“^). 

The cylindrical part of the fuselage has diameter 5,6 m and length 65 m. 

1 . What is approximately the maximum boundary layer thickness around the fuse¬ 
lage? 

2 . What is approximately the average shear applying on the fuselage skin? 

3. Estimate the power dissipated to friction on the cylindrical part of the fuselage. 

4. In practice, in which circumstances could flow separation occur on the fuselage 
skin? (give a brief answer, e.g. in 30 words or less) 


7.6 Separation according to Pohlhausen 

^ non-examinable, based on Richecoeur 2012 [17] 

Air at 1 bar and 20 °C flows along a smooth surface, and decelerates slowly, with a 
constant rate of -0,25 m s“^ m“h 

According to the Pohlhausen model (eq. 7/21 p.l54), at which distance downstream will 
separation occur? Is the boundary layer still laminar then? 

How could one generate such a deceleration in practice? 


7.7 Laminar wing profile 

^ ^ non-examinable. Based on a diagram from Bertin et al 2010 [15] 

The characteristics of a so-called “laminar” wing profile are compared in figs. 7.15 to 7.17 
with those of an ordinary profile. 


On the graph representing the pressure coefficient Cp 
spending to each profile. 


P poo 

'W 


, identify the curve corre- 


What advantages and disadvantages does the laminar wing profile have, and how can 
they be explained? In which applications will it be most useful? 



Figure 7.15 - Comparison of the thickness distribution of two uncambered wing profiles: an 
ordinary medium-speed naca 0009 profile, and a “laminar” naca 66-009 profile. 

Figure © Bertin & Cummings 2010 [15] 
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Figure 7.16 - Static pressure distribution (represented as a the local non-dimensional pressure 
coefhcient Cp = f"^) as a function of distance x (non-dimensionalized with the chord c) over 

2pV 

the surface of the two airfoils shown in fig. 7.15. 

Figure © Bertin & Cummings 2010 [15] 


□ NACA 0009 
ONACA 66-009 



Figure 7.17 - Values of the section drag coefficient C 4 
coefficient C/ = for both airfoils presented in fig. 


_ d 
\cpV‘^ 

7.15. 


as a function of the section lift 


Figure © Bertin & Cummings 2010 [15], based on data by Abott & Von Doenhoff 1949 [1 ] 
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7.8 Separation mechanism 

non-examinable, CC-0 o.c. 

Sketch the velocity profile of a laminar or turbulent boundary layer shortly upstream of, 
and at a separation point. 

The two equations below describe flow in laminar boundary layer: 


du du dU u d^u 

u — x) — — U -"t - T 

dx dy dx p {dyY 

du dv 

dx dy 


(7/12) 

(7/13) 


Identify these two equations, list the conditions in which they apply, and explain shortly 
(e.g. in 30 words or less) why a boundary layer cannot separate when a favorable pressure 
gradient is applied along the wall. 


7.9 Thwaites’ separation model 

non-examinable, from 'White [ 13] E7.5 

In 1949, Bryan Thwaites explored the limits of the Pohlhausen separation model. He 
proposed a different model to describe the laminar boundary layer velocity profile, which 
is articulated upon an expression for the momentum thickness 9 and is more accurate 
than the Reynolds-number based descriptions that we have studied. Thwaites proposed 
the model: 

= +0,45^^J^ U^dx (7/28) 

in which the 0-subscripts denote initial conditions (usually 6o = 0) 

Thwaites also defined a non-dimensional momentum thickness Ag: 

Ag = 0^-^ (7/29) 

p dx 

and then went on to show that for the model above, separation (cy = 0) occurs when 
Xe = -0,09. 

A classical model for flow deceleration is the Howarth longitudinal velocity profile 
U(x) = (l - i)> ill which L is a reference length of choice. In this velocity distribution 
with linear deceleration, what is the distance at which the Thwaites model predicts the 
boundary layer separation? 
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Answers 


7.1 1) At trailing edge [ReJ^f = 5 348 thus the layer is laminar everywhere. 5 will grow 

from 0 to 2,01 cm (eq. 7/14a p.l53); 2) For water: ^trailing edge = 4,91 mm. 

7.2 1) See fig. 7.7 p.l55. At the leading-edge the velocity is uniform. Note that the 

y-direction is greatly exaggerated, and that the outer velocity U is identical for 
both regimes; 2) See fig. 7.4 p.l50. Note that stream¬ 
lines penetrate the boundary layer; 3) and 4) See §7.4 p. 155. 


7.3 


^transition, air = 4,898 m and Xtransition, water = 0,4 m. In a laminar boundary layer, 
inserting equation 7/14d into equation 7/5 into equation 7/6 yields 

f.= 0,664LC7>.= Vw[Vr];;’'""“. 

In a turbulent boundary layer, we use equation 7/23d instead and get 

= 0,01575Lp7[/^/t^ [x?] These expressions allow the calculation of 

^transition 

the forces below, for the top surface of the plate: 

1) (air) First case: F = 3,445 • 10“^ N; 
second case F = 8,438 • 10“^ N (who would have thought eh?); 

2 ) (water) First case: F = 3,7849 N; sec¬ 


ond case F = 2,7156 N. 


7.4 Using the expressions developed in exercise 7.3, Wfuction ~ 255 W. 


7.5 1) Xtransition = 7,47 cm (the laminar part is negligible). With the equations developed 

in exercise 7.3, we get F = 24,979 kN and W = 6,09 MW. Quite a jump from the 
Wright Flyer I! 

2 ) When the longitudinal pressure gradient is zero, the boundary layer cannot 
separate. Thus separation from the fuselage skin can only happen if the fuselage is 
flown at an angle relative to the flight direction (e.g. during a low-speed maneuver). 


7.9 


Once the puzzle pieces are put together, this is an algebra exercise: 
0,1231. Bryan beats Ernst! 



separation 
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Fluid Mechanics 

Chapter 8 - Modeling large- and small-scale flows 

last edited October 30, 2017 
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8.3.5 Modeling lift with circulation 177 

8.4 Flow at very small scales 180 

8.5 Exercises 183 


These lecture notes are based on textbooks by White [13], Qengel & al.[16], and Munson & al.[18]. 


8.1 Motivation 


This exploratory chapter is not a critical component of fluid dynamics; in¬ 
stead, it is meant as a brief overview of two extreme cases: flows for which 
viscous effects are negligible, and flows for which they are dominant. This 
exploration should allow us to answer two questions: 

• How can we model large-scale flows analytically? 

• How can we model small-scale flows analytically? 


8.2 Flow at large scales 


8.2.1 Problem statement 


In this section, we are interested in flow at very large scales: those for which 
a representative length L is very large. In particular, when L is large enough, 
the influence of viscosity is reduced. Formally, this corresponds to the case 
where the Reynolds number [Re] = pVT/p is very large. 

To examine the mechanics of such a flow, we turn to our beloved non- 
dimensional Navier-Stokes equation for incompressible flow derived as 
eq. 6/17p.l35, 


[St] 


dV* 

'W 


+ [1] V* • r V* 


^ r - [Eu] vy + ^ V* V* 

[Fr]^ [Re] 

( 8 / 1 ) 



Video: pre-lecture briefing for 
this chapter 

by o.c. (CC-by) 
https://youtu.be/qEE-tyXgWLI 
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in which we see that when the Reynolds number [Re] is very large, the last 
term becomes negligible relative to the other four terms. Thus, our governing 
equation can be reduced as follows: 

dV* ^ ^ 1 

[St] ^ + [1] y* • «-r g - [Eu] Vy (8/2) 

at* [Fr]"^ 

Now, converting eq. 8/2 back to dimensional terms, our governing momentum 
equation for inviscid flow becomes: 

- - .X 

P-^ = pg-^p (8/3) 

We see that with the starting proposition that [Re] was large, we have re¬ 
moved altogether the viscous (last) term from the Navier-Stokes equation. 
Flows governed by this equation are called inviscid flows. Equation 8/3 is 
named the Euler equation; it stipulates that the acceleration field is driven 
only by gravity and by the pressure field. 

8.2.2 Investigation of inviscid flows 

From what we studied in chapter 6, we recognize immediately that flows 
governed by eq. 8/3 are troublesome: the absence of viscous effects facilitates 
the occurrence of turbulence and makes for much more chaotic behaviors. 
Although the removal of shear from the Navier-Stokes equation simplifies the 
governing equation, the solutions to this new equation become even harder 
to find and describe. 

What can be done in the other two branches of fluid mechanics? 

• Large-scale flows are difficult to investigate experimentally. As we have 
seen in chapter 6, scaling down a flow (e.g. so it may fit inside a labora¬ 
tory) while maintaining constant [Re] requires increasing velocity by a 
corresponding factor. 

• Large scale flows are also difficult to investigate numerically. At high 
[Re], the occurrence of turbulence makes for either an exponential 
increase in computing power (direct numerical simulation comput¬ 
ing power increases with [Re]^’^) or for increased reliance on hard-to- 
calibrate turbulence models (in Reynolds-averaged simulations). 

All three branches of fluid mechanics, therefore, struggle with turbulent 
flows. 

In spite of this, large-scale flows are undeniably important. In chapter 7 we 
were able to understand and describe fluid flow very close to walls. Now we 
wish to be able to to the same for large structures, for example, in order to 
describe the broad patterns of fluid flow (in particular, pressure distribution) 
in the wake of an aircraft or of a wind turbine, or within a hurricane. It is 
clear that we have no hope of accounting easily for turbulence, but we can 
at least describe the main features of such flows by restricting ourselves to 
laminar cases. In the following sections, we will model such laminar solutions 
directly, based on intuition and observation, and make sure that they match 
the condition described by eq. 8/3 above. 
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8.3 Plotting velocity with functions 


8.3.1 Kinematics that fit conservation laws 

For simple flow structures, the velocity field can be simply described based 
on observation and intuition. This is done for example in exercise 8.8, where 
we reconstruct the flow field within and around a tornado using very simple, 
almost primitive, kinematics. 

Without so much as a small increase in complexity, this approach becomes 
untenable. In the last exercises of chapter 4 (ex. 4.7 & 4.8 p.l05) we have seen 
that it is easy to propose a velocity field that does not satisfy (is not a solution 
of) either the continuity or the Navier-Stokes equation. For example, if one 
considers two of the tornado flows mentioned above together, the velocity 
field cannot be described easily anymore. 

One approach has been developed in the 17* century to overcome this 
problem. It consists in finding a family of flows, all steady, that always satisfy 
the conservation equations. Those flows can then be added to one another 
to produce new conservation-abiding flows. Such flows are called potential 
flows. 

Two conditions need to be fulfilled for this approach to work: 

1. The velocity field must always be describable with a function; that is, 
there must correspond a single value of u, of v and of w at each of 
the coordinates x,, y,, or z, in space. This means in practice that we 
cannot account for flows which “curl up” on themselves, occasionally 
recirculating back on their path. 

2. The velocity field must conserve mass. In incompressible flow, this is 

—^ ^ 

achieved if the continuity equation V • V = 0 (eq. 4/17 p.90) is respected. 

With potential flow, these two conditions are addressed as follows: 

1. We restrict ourselves to irrotational flows, those in which the curl of 
velocity (see Appendix A2 p.217) is always null: 


Vx V = 0 


(8/4) 


by definition, for an irrotational flow. 

It can be shown that flows are irrotational when there exists a scalar 


function (j) (pronounced “phi” and named potential function) of which 
the gradient is the velocity vector field: 



(8/5) 


In the case of two-dimensional flow, this translates as: 

d(j) _ 


— = u 


( 8 / 6 ) 


d(f) 


(8/7) 


dy 


When plotted out, lines of constant ^ (named equipotential lines) are 
always perpendicular to the streamlines of the flow. 
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2. The continuity equation is satisfied by referring to stream functions. It 
can be shown that the divergent of velocity is null when there exists a 
vector field function ij/ (pronounced “psi” and named stream function) 
of which the curl is the velocity vector field: 

Vxf=V (8/8) 


In the case of two-dimensional flow, ^ is a scalar field and eq. 8/8 
translates as: 




(8/9) 

( 8 / 10 ) 


When plotted out, lines of of constant \]/ value are streamlines - in other 
words, as they travel along, fluid particles follow paths of constant ^ 
value. 


In summary, we have shifted the problem from looking for u and v, to looking 
for ^ and (j). The existence of such functions ensures that flows can be added 
and subtracted from one another yet will always result in mass-conserving, 
mathematically-describable flows. If such two functions are known, then the 
velocity components can be obtained easily either in Cartesian coordinates. 


dij/ 

dx dy 
d(j) dxj/ 

or angular coordinates: 

d(j) 1 dxj/ 

^ dr r dO 

1 d^ dxj/ 


( 8 / 11 ) 

( 8 / 12 ) 


(8/13) 

(8/14) 


8.3.2 Strengths and weaknesses of potential flow 

The potential flow methodology allows us to find solutions to the Euler 
equation: flows in which the Reynolds number is high enough that viscosity 
has no significant role anymore. Potential flows are the simplest solutions 
that we are able to come up with. They are: 

• strictly steady; 

• inviscid; 

• incompressible; 

• devoid of energy transfers; 

• two-dimensional (in the scope of this course at least). 

This is quite convenient for the academician, who recognizes immediately 
four of the five criteria which we set forth in chapter 3 for using the Bernoulli 


equation. Along a streamline, the fifth condition is met, and Euler’s equation 
reduces to eq. 3/18 (p.69): 

— + -V]^ + gzi = — + -V 2 + gz 2 = cst. (8/15) 

p 2 P 2 

along a streamline in a steady incompressible inviscid flow. 

and so it follows that if the solution to an inviscid flow is known, the forces 
due to pressure can be calculated with relative ease. 

Nevertheless, from a science and engineering point of view, potential flows 
have only limited value, because they are entirely unable to account for 
turbulence, which we have seen is an integral feature of high-[Re] flows. 
We should therefore use them only with great caution. Potential flows help 
us model large-scale structures with very little computational cost, but this 
comes with strong limitations. 

8.3.3 Superposition: the lifting cylinder 

It is possible to describe a handful of basic potential flows called elementary- 
flows as fundamental ingredients that can be added to one another to create 
more complex and interesting flows. Without going into much detail, the 
most relevant elementary flows are: 

• Uniform longitudinal flow. 


(j) = V rcosd (8/16) 

^ = V r sinO (8/17) 


• Sources and sinks (fig. 8.1) which are associated with the appearance 
of a (positive or negative) volume flow rate W from a single point in 
the flow: 


— Inr 

2 n 

(8/18) 

In 

(8/19) 


• Irrotational vortices, rotational patterns which impart a rotational 
velocity V 0 = /(r, 9) on the flow, in addition to which there may exist 



Figure 8.1 - Concept of a source inside a two-dimensional potential flow. A sink 
would display exactly opposed velocities. 

Figure CC-by-sa Commons User.Nicoguaro 
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a radial component v^: 


(f> = ^e ( 8 / 20 ) 

In 

r 

f = -Inr (8/21) 

2 n 

in which F (termed circulation) is a constant proportional to the strength of 
the vortex. 


Doublets, which consist in a source and a sink of equal volume flow 
rate positioned extremely close one to another (fig. 8.2); 


^ = K 


cos 6 


\l/ = -K 


r 

sin0 


( 8 / 22 ) 

(8/23) 


in which is a constant proportional to the source/sink volume flow rate 'V. 


It was found in the 17* Century that combining a doublet with uniform 
flow resulted in flow patterns that matched that of the (inviscid, steady, 
irrotational) flow around a cylinder (fig. 8.3). That flow has the stream 
function: 


^ = Uco sin 6 



(8/24) 


This stream function allows us to describe the velocity everywhere: 

Idil/ ^ ( R^\ 

Ur = = UooCOse{l - — \ 8/25) 

r da \ } 

dll/ / B?\ 

U0 =-—=-UooSinO il +—\ (8/26) 

We can even calculate the lift and drag applying on the cylinder surface. 
Indeed, along the cylinder wall, r = R and 

tlr\r=R ~ b (8/27) 

=-2UooSin9 (8/28) 

Since the Bernoulli equation can be applied along any streamline in this 
(steady, constant-energy, inviscid, incompressible) flow, we can express the 



Figure 8.2 - Left: a simple uniform steady flow; Right: a doublet, the result of a 
source and a sink brought very close one to another 

Figure CC-0 o.c. 


174 

















Figure 8.3 - The addition of a doublet emd a uniform flow produces streamlines for 
an (idealized) flow around a cylinder. 


Figure CC-by-sa by Commons User.Kraaiennest 


pressure ps on the cylinder surface as a function of 0 : 



(8/29) 


Now a simple integration yields the net forces exerted by the fluid on the 
cylinder per unit width L, in each of the two directions x and y. 



(8/30) 


ps cos 6 RAO = Q 


(8/31) 


ps sinO RdO = 0 


The results are interesting, and at the time they were obtained by their author, 
Jean le Rond D’Alembert, were devastating: both lift and drag are zero. This 
inability to reproduce the well-known phenomena of drag is often called the 
d’Alembert paradox. 

The pressure distribution on the cylinder surface can be represented graphi¬ 
cally and compared to experimental measurements (flg. 8.4). Good agreement 
is obtained with measurements on the leading edge of the cylinder; but as 
the pressure gradient becomes unfavorable, in practice the boundary layer 
separates -a phenomenon that cannot be described with inviscid flow— and 
a low-pressure area forms on the downstream side of the cylinder. 



















Figure 8.4 - Pressure distribution (relative to the far-flow pressure) on the surface of 
a cylinder, with flow from left to right. On the left is the potential flow case, purely 
symmetrical. On the right (in blue) is a measurement made at a high Reynolds 
number. Boundary layer separation occurs on the second half of the cylinder, which 
prevents the recovery of leading-edge pressure values, and increases drag. 

Figure CC-by-sa Commons User.BoH & Olivier Cleynen 


8.3.4 Circulating cylinder 



Video: watch the Brazilian foot¬ 
ball team show their French 
counterparts how circulation 
(induced through friction by 
ball rotation) is associated to dy¬ 
namic lift on a circular body 

by TFl, 1997 (styl) 
https://youtu.be/oGeMZ3t8jn4 


An extremely interesting hack can be played on the potential cylinder flow 
above if an irrotational vortex of stream function ^ Inr (eq. 8/21) is 

added to it. The overall flow field becomes: 

\ Y 

r --In r (8/32) 

r ) 2n 

With this function, several key characteristics of the flow field can be obtained. 
The first is the velocity field at the cylinder surface: 

^r\r—R ~ b (8/33) 

r 

ue\r=R = -2Uoo sin 6 + —— (8/34) 

2nR 


^ = {7co sin 9 


and we immediately notice that the velocity distribution is no longer sym¬ 
metrical with respect to the horizontal axis (fig. 8.5): the fluid is deflected, 
and so there will be a net force on the cylinder. 

This time, the net pressure forces on the cylinder have changed: 

coseRde = 0 (8/35) 

sin 9 Rd9 = -p Uoo F (8/36) 



We thus find out that the drag is once again zero —as for any potential 
flow— but that lift occurs which is proportional to the free-stream velocity U 
and to the circulation F. 

In practice, such a flow can be generated by spinning a cylindrical part in 
a uniform flow. A lateral force is then obtained, which can be used as a 
propulsive or sustaining force. Several boats and even an aircraft have been 
used in practice to demonstrate this principle. Naturally, flow separation 
from the cylinder profile and the high shear efforts generated on the surface 











Figure 8.5 - The addition of an irrotational vortex on top of the cylinder flow 
described in fig. 8.3 distorts the flow field and it becomes asymmetrical: a lift force is 
developed, which depends directly on the circulation F. 

Figure © White 2008 [13] 


cause real flows to differ from fhe ideal case described here, and if turns out 
that rotating cylinders are a horribly uneconomical and unpractical way of 
generating lift. 


8.3.5 Modeling lift with circulation 


Fluid flow around cylinders may have little appeal for the modern student 
of fluid mechanics, buf the methodology above has been taken much further. 
With further mathematical manipulation called conformal mapping, potential 
flow can be used to described flow around geometrical shapes such as airfoils 
(fig. 8.6). Because the flow around such streamlined bodies usually does 
not feature boundary layer separation, the predicted flow fields everywhere 
except in the close vicinity of the solid surface are accurately predicted. 

There again, it is observed that regardless of the constructed geometry, no lift 
can be modeled unless circulation is also added within the flow. The amount 
of circulation needed so that results may correspond to experimental obser¬ 
vations is found by increasing it progressively until the the rear stagnation 
point reaches the rear trailing edge of the airfoil, a condition known as the 
Kutta condition. Regardless of the amount of circulation added, potential 
flow remains entirely reversible, both in a kinematic and a thermodynamic 
sense, thus, care must be taken in the problem setup to make sure the model 
is realistic (fig. 8.7). 



Video: acting on swimming 
pool water with a round plate 
sheds a half-circular vortex that 
is extremely stable and can 
be interacted with quite easily. 
Such stable laminar structures 
are excellent candidates for 
analysis using potential flow. 

by Physics Girl (Dianna Cowern) (styl) 
https://youtu.be/pnbJEg9rlo8 


It is then observed in general that any dynamic lift generation can be modeled 
as the superposition of a free-stream flow and a circulation effect (fig. 8.8). 
Wifh such a fool, pofential flow becomes an extremely useful tool, mathemat¬ 
ically and computationally inexpensive, in order to model and understand the 
cause and effect of dynamic lift in fluid mechanics. In particular, it has been 
paramount in the description of aerodynamic lift distribution over aircraft 
wing surfaces (fig. 8.9), with a concept called the Lifting-line theory. 
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Figure 8.6 - Potential flow around an airfoil without (top) and with (bottom) circu¬ 
lation. Much like potential flow around a cylinder, potential flow around an airfoil 
can only result in a net vertical force if an irrotational vortex (with circulation F) is 
added on top of the flow. Only one value for F will generate a realistic flow, with the 
rear stagnation point coinciding with the trailing edge, a occurrence named Kutta 
condition. 


Figure CC-0 o.c. 



Figure 8.7 - Potential flow allows all velocities to be inverted without any change in 
the flow geometry. Here the flow around an airfoil is reversed, displaying unphysical 
behavior. 


Figure CC-0 o.c. 



Figure 8.8 - A numerical model of flow around an airfoil. In the left figure, the 
velocity vectors are represented relative to a stationary background. In the right 
figure, the velocity of the free-stream flow has been subtracted from each vector, 
bringing the circulation phenomenon into evidence. 


Figures 1 & 2 CC-by-sa by en:Wikipedia User.Crowsnest 


























































Figure 8.9 - From top to bottom, the lift distribution over the wings of a glider is 
modeled with increasingly complex (and accurate) lift and circulation distributions 
along the span. The Lifting-line theory is a method associating each element of 
lift with a certain amount of circulation. The effect of each span-wise change of 
circulation is then mapped onto the flow field as a trailing vortex. 

Figures 1, 2 & 3 CC-by-sa Olivier Cleynen 
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8.4 Flow at very small scales 


At the complete opposite of the spectrum, we find flow at very small scales: 
flows around bacteria, dust particles, and inside very small ducts. In those 
flows the representative length L is extremely small, which makes for small 
values of the Reynolds number. Such flows are termed creeping or Stokes 
flows. What are their main characteristics? 

Looking back once again at the non-dimensional Navier-Stokes equation for 
incompressible flow derived as eq. 6/17 p.l35. 


Ot 7'>N -1 

[St] ^ + [1] • Tv* = 


dt* 


[Fr]‘ 


g 


[Eu]vy + —(8/37) 


we see that creeping flow will occur when the Reynolds number is much 
smaller than 1. The relative weight of the term ^ V*^V* then becomes 
overwhelming. 

In addition to cases where [Re] <?c 1, we focus our interest on flows for which: 


• gravitational effects have negligible influence over the velocity field; 

• the characteristic frequency is extremely low (quasi-steady flow). 

Wifh these characteristics, the terms associated with the Strouhal and Froude 
numbers become very small with respect to the other terms, and our non- 
dimensionalized Navier-Stokes equation (eq. 8/37) is approximately reduced 
to: 


0 « -[Eu] Tp* + ^ T^V* (8/38) 

We can now come back to dimensionalized equations, concluding that for a 
fluid flow dominated by viscosity, the pressure and velocity fields are linked 
together by the approximate relation: 


Vp = pV^V (8/39) 

In this type of flow, the pressure field is entirely dictated by the Laplacian 
of velocity, and the fluid density has no importance. Micro-organisms, for 
which the representative length L is very small, spend their lives in such 
flows (fig. 8.10). At the human scale, we can visualize the effects of these 
flows by moving an object slowly in highly-viscous fluids (e.g. a spoon in 
honey), or by swimming in a pool filled wifh plastic balls. The inertial effects 
are almost inexistent, drag is extremely important, and the object geometry 
has comparatively small influence. 

In 1851, George Gabriel Sfokes worked fhrough equation 8/39 for flow around 
a sphere, and obtained an analytical solution for the flow field. This allowed 
him to show that the drag Fd sphere applying on a sphere of diameter D in 
creeping flow (fig. 8.11) is: 


Fd sphere — 37rpUooD 


(8/40) 
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Figure 8.10 - Micro-organisms carry themselves through fluids at extremely low 
Reynolds numbers, since their scale L is very small. For them, viscosity effects 
dominate inertial effects. 

Photo by Yutaka Tsutsumi, M.D., Fujita Health University School of Medicine 


Inserting this equation 8/40 into the definition of the drag coefficient CfD = 
PollpSirontalU^ (from eq. 6/2 p.l30) then yields: 


Cf D = 


Pd sphere 




24/i 

p UcoD 


24 


(8/41) 


These equations are specific to flow around spheres, but the trends they 
describe apply well to most bodies evolving in highly-viscous flows, such 
as dust or liquid particles traveling through the atmosphere. Drag is only 
proportional to the speed (as opposed to low-viscosity flows in which it 
grows with velocity squared), and it does not depend on fluid density 



Figure 8.11 - Flow at very low Reynolds numbers around a sphere. In this regime, 
the drag force is proportional to the velocity. 

Figure CC-by-sa by Olivier Cleynen & Commons User.Kraaiennest 
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Exercise sheet 8 (Modeling large- and small-scale flows) 

last edited June 18, 2017 


These lecture notes are based on textbooks by White [13], Qengel & al.[16], and Munson & al.[18]. 


Except otherwise indicated, we assume that fluids are Newtonian, and that: 

Pwater = 1 000 kg m“^ patm. = 1 bar; Patm. = 1,225 kgm“^ Tatm. = ll,3°C;patm. = 1,5 • 10‘^Nsm“2; 
g = 9,81 ms“^. Air is modeled as aperfect gas (i?air = 287JK“^ kg“'; pair = 1,4; Cpair = 1005Jkg“^ K“^). 


In a highly-viscous (creeping) steady flow, the drag Fd exerted on a spherical body of 
diameter D at by flow at velocity Uoo is quantified as: 

Fd sphere = 3;Tpl7oo-D (8/40) 


8.1 Volcanic ash from the Eyiafiallaiokull 

Qengel & al. [16] E10.2 

In 2010, a volcano with a complicated name and unpredictable mood decided to ground 
the entire European airline industry for five days. 

We consider a microscopic ash particle released at very high altitude (-50 °C, 0,55 bar, 
1,474 • 10“^Nsm“^). We model it as a sphere with 50 pm diameter. The density of 
volcanic ash is 1 240 kgm“^. 

1. What is the terminal velocity of the particle? 

2. Will this terminal velocity increase or decrease as the particle progresses towards 
the ground? (give a brief, qualitative answer) 


8.2 Water drop 

Qengel & al. [16] 10-21 

A rainy day provides yet another opportunity for exploring fluid dynamics (fig. 8.12). A 
water drop with diameter 42,4 pm is falling through air at 25 °C and 1 bar. 

1. Which terminal velocity will it reach? 

2. Which velocity will it reach once its diameter will have doubled? 
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< 5 mm 




> 5 mm 


(D 


Figure 8.12 - A sketched diagram showing the geometry of water drops of various sizes in free 
fall. When their diameter is lower than 2 mm, water drops are approximately spherical (B). As 
they grow beyond this size, their shape changes and they eventually break-up (C-E). They never 
display the “classical” shape displayed in A, which is caused only by surface tension effects when 
they drip from solid surfaces. 

Figure CC-by-sa by Ryan Wilson 


8.3 Hangar roof 


based on White [13] P8.54 


Certain flows in which both compressibility and viscosity effects are negligible can be 
described using the potential flow assumption (the hypothesis that the flow is everywhere 
irrotational). If we compute the two-dimensional laminar steady fluid flow around a 
cylinder profile, we obtain the velocities in polar coordinates as: 


Uy — — Uqo cos U I 1 17 

r o9 




U0 = 




= -Ucx> sin d 1 + 


where the origin (r = 0) is at the center of the cylinder profile; 

6 is measured relative to the free-stream velocity vector; 
Uco is the incoming free-stream velocity; 
and R is the (fixed) cylinder radius. 


(8/25) 

(8/26) 


Based on this model, in this exercise, we study the flow over a hangar roof 

Wind with a nearly-uniform velocity Uc^ = 100kmh“^ is blowing across a 50m-long 
hangar with a semi-cylindrical geometry, as shown in fig. 8.13. The radius of the hangar 
is J? = 20 m. 


1. If the pressure inside the hangar is the same as the pressure of the faraway 
atmosphere, and if the wind closely follows the hangar roof geometry (without 
any flow separation), what is the total lift force on the hangar? 

(hint: we accept that f sin^ xdx = ^ cos^ x - cos x + k). 
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Figure 8.13 - A semi-cylindrical hangar roof. Wind with uniform velocity U flows perpendicular 
to the cylinder axis. 

Figure CC-0 o.c. 


2. At which position on the roof could we drill a hole to negate the aerodynamic lift 
force? 

3. Propose two reasons why the aerodynamic force measured in practice on the 
hangar roof may be lower than calculated with this model. 


8.4 Cabling of the Wright Flyer 

c> y derived from Munson & al. [18] 9.106 

The Wright Flyer I, the first powered and controlled aircraft in history, was subjected 
to multiple types of drag. We have already studied viscous friction on its thin wings in 
exercise 7.4. The data in figure 8.14 provides the opportunity to quantify drag due to 
pressure. 

A network of metal cables with diameter 1,27 mm criss-crossed the aircraft in order to 
provide structural rigidity. The cables were positioned perpendicularly to the air flow, 
which came at 40kmh“h The total cable length was approximately 60 m. 

What was the drag generated by the cables? 
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{A) (B) 



Oscillating Karman vortex street wake 
(C) 



Laminar boundary layer, 
wide turbulent wake 



Turbulent boundary layer, 
narrow turbulent wake 


iD) 


iE) 


Figure 8.14 - Experimental measurements of the drag coefficient applying to a cylinder and 
to a sphere as a function of the diameter-based Reynolds number [Re]^), shown together with 
schematic depictions of the flow around the cylinder. By convention, the drag coefficient 

Cb = Cp B = 1 6/2 p.l30) compares the drag force Fb with the frontal area S. 

Both figures © from Munson & al.[18] 
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8.5 Ping pong ball 


Munson & al. [is] E9.16 


A series of experiments is conducted in a wind tunnel on a large cast iron ball with a 
smooth surface; the results are shown in fig. 8.15. These measurement data are used to 
predict the behavior of a ping pong ball. Table tennis regulations constrain the mass of 
the ball to 2,7 g and its diameter to 40 mm. 


1. Is it possible for a ball thrown at a speed of 50 kmh ^ to have a perfectly horizontal 
trajectory? 

2. If so, what would be its deceleration? 

3. How would the drag and lift applying on the ball evolve if the air viscosity was 
progressively decreased to zero? 



CODI2U 


Figure 8.15 - Experimental measurements of the lift and drag coefficients applying on a rotating 
sphere in an steady uniform flow. 

Figure © from Munson & aL[18] 
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8.6 Lift on a symmetrical object 

non-examinable 

Briefly explain (e.g. with answers 30 words or less) how lift can be generated on a sphere 
or a cylinder, 

• with differential control boundary layer control; 

• with the effect of rotation. 

Draw a few streamlines in a two-dimensional sketch of the phenomenon. 


8.7 Air flow over a wing profile 

^ Non-examinable. From Munson & al. [18] 9.109 

The characteristics of a thin, flat-bottomed airfoil are examined by a group of students 
in a wind tunnel. The first investigations focus on the boundary layer, and the research 
group evaluate the boundary layer thickness and make sure that it is fully attached. 

Once this is done, the group proceeds with speed measurements all around the airfoil. 
Measurements of the longitudinal speed u just above the boundary layer on the top 
surface are tabulated below: 


x/c (%) 

y/c (%) 

u/U 

0 

0 

0 

2,5 

3,72 

0,971 

5 

5,3 

1,232 

7,5 

6,48 

1,273 

10 

7,43 

1,271 

20 

9,92 

1,276 

30 

11,14 

1,295 

40 

10,49 

1,307 

50 

10,45 

1,308 

60 

9,11 

1,195 

70 

6,46 

1,065 

80 

3,62 

0,945 

90 

1,26 

0,856 

100 

0 

0,807 


On the bottom surface, the speed is measured as being constant {u = U) to within 
experimental error. 

What is the lift coefficient on the airfoil? 
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8.8 Flow field of a tornado 


From Qengel & al. [16] E9-5, E9-14 & ElO-3 

The continuity and Navier-Stokes equations for incompressible flow written in cylindrical 
coordinates are as follows: 


dv. 


r V0 dVr dVy 

- “h Z)r - “F-— - “F 'O7 - 

dt dr r 89 r " dz 


dVr 

dp 




1 d 


r dr 

1 dr) 


Vr 1 d^^Vr 2 dVQ d^Vr 

{dOy 89 {dzy 
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1 d^v 


z d'^Vz 
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(8/42) 


(8/43) 


(8/44) 

(8/45) 


In this exercise, we are interested in solving the pressure field in a simplified description 
of a tornado. For this, we consider only a horizontal layer of the flow, and we consider 
that all properties are independent of the altitude z and of the time t. 

We start by modeling the tornado as a vortex imparting an angular velocity such that: 


U0 = 


T 

2nr 


in which F is the circulation (measured in s h and remains constant everywhere. 


(8/46) 


1. What form must the radial velocity Uf take in order to satisfy continuity? 

Among all the possibilities for Ur, we choose the simplest in our study, so that: 

Ur = 0 (8/47) 

all throughout the tornado flow field. 

2. What is the pressure field throughout the horizontal layer of the tornado? 

The flow field described above becomes unphysical in the very center of the tornado 
vortex. Indeed, our model for uq is typical of an irrotational vortex, which, like all 
irrotational flows, is constructed under the premise that the flow is inviscid. However, in 
the center of the vortex, we are confronted with high velocity gradients over very small 
distances: viscous effects can no longer be neglected and our model breaks down. 

It is observed that in most such vortices, a vortex core region forms that rotates just like 
a solid cylindrical body. This flow is rotational and its governing equations result in a 
realistic pressure distribution. 
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3. What is the pressure field within the rotational core of the tornado? 

Make a simple sketch showing the pressure distribution as a function of radius 
throughout the entire tornado flow field. 

4. [An opening question for eurious students] 

What is the basic mechanism of vortex stretching? How would it modify the flow 
field described here? 
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Answers 


8.1 1) At terminal velocity, the weight of the sphere equals the drag. This allows us to 
obtain U = gpspheref^ = 0,1146 m s“h unbearably slow when you are stuck in an 
airport! With U, check that the Reynolds number indeed corresponds to creeping 
flow: [Re]^) = 0,334. 

8.2 Same as previous exercise: Ui = 4,578 • 10“^ m s“^ and U 2 = 0,183 m s“^, with Reynolds 
numbers of 0,113 and 0,906 respectively (thus creeping flow hypothesis valid). 

8.3 1) express roof pressure as a function of 9 using eq. 8/29 p.l75 on eq. 8/26, then 
integrate the vertical component of force due to pressure: Fp roof = 1,575 MN. 

2) 9\p=o = 54,7° 

8.4 A simple reading gives Fd = 6,9 N, W = 76 W. 

8.5 Yes — a reading gives ro = 83 rev/s. 
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Fluid Mechanics 

Chapter 9 - Introduction to compressible flow 

last edited April 2, 2018 
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9.8 Exercises 207 


These lecture notes are based on textbooks by White [13], Qengel & al.[16], and Munson & al.[18]. 


9.1 Motivation 


In this chapter we investigate the fundamentals of compressible flows: those 
in which density varies significantly. This study should allow us to answer 
two questions: 

• What happens in ideal (reversible) compressible flows? 

• What is a shock wave, and how does it influence a flow? 



Video: pre-lecture briefing for 
this chapter 

by o.c. (CC-by) 
https ://youtu.be/Lk8BMRC_Cy w 


9.2 Compressibility and its consequences 


9.2.1 Problem description 

Ever since we have begun describing fluid flow in an extensive manner in 
chapter 4, we have restricted ourselves to incompressible flows — those in 
which the fluid density p is uniform and constant. Now, p becomes one 
additional unknown property field. 

In some cases, solving for an additional unknown is not particularly challeng¬ 
ing. For example, if we wish to calculate water temperature in a heat transfer 
case, or species concentration in a reacting flow, it is possible to first solve 
for the flow, and then solve for those particular property fields within the 
solution; thus the additional unknowns are solved separately with additional 
equations. 

The mathematical treatment of density, however, is much worse, because p 

appears in the momentum and mass conservation equations. Density directly 

—^ 

affects the velocity fields, and so we have solve for p, p and V simultaneously. 
From a physical point of view, compressibility brings two new effects: 
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• Density changes in fluids translate into large changes in fluid mechani¬ 
cal energy; these can be come very large in comparison to other forms 
of fluid energy, in particular kinetic energy. When p varies, fluids 
“perform work upon themselves”, accumulating and expending this 
energy, which translates into kinetic energy changes, and irreversible 
losses as heat transfer. 

• When density is allowed to vary, pressure travels at finite speed within 
the fluid — one could say that the fluid becomes appreciably “squishy”. 
The delay resulting from this finite pressure propagation speed be¬ 
comes very significant if the fluid velocity is high, and it may even 
prevent the propagation of pressure waves in the upstream direction 
altogether. 

In this chapter, we wish to quantify both of those effects. By now, after 
eight chapters already where p was just an input constant, the reader should 
readily agree we have no hope of coming up with a general solution to 
compressible flow. Nevertheless, it is possible to reproduce and quantify 
the most important aspects of compressibility effects with relatively simple, 
one-dimensional air flow cases — and this will suffice for us in this last 
exploratory chapter. 

9.2.2 The speed of pressure 

We start by considering the speed of pressure. In compressible fluids, pressure 
changes propagate with a speed that increases with the magnitude of the 
pressure change. When this magnitude tends to zero (and the pressure wave 
reduces to that of a weak sound wave, with fully-reversible fluid compression 
and expansion) then the speed tends to a given speed which we name speed 
of sound and noted c. At ambient conditions, this speed is approximately 
1 000 km/h in air and roughly 5 000 km/h in water. 

Let us, as a thought experiment, create a small pressure change inside a static 
fluid, and then travel along with the pressure wave. As we follow a wave 
which travels from right to left, we perceive fluid movement from left to right 
(fig. 9.1). We construct a small control volume encompassing and moving 
with the wave — at the speed of sound c. 

Control volume fluid flow analysis does not frighten us anymore. From mass 
conservation (eq. 3/5 p.64) we write: 

pcA ={p+ dp){c - dV)A (9/1) 


flow moving 
from left to right 


cross-section area A 


P 

P 


c + dV 

-► 


p + dp 
p + dp 


control volume around pressure wave 


Figure 9.1 - An infinitesimal pressure-difference wave is represented traveling from 
right to left inside a stationary fluid, in a one-dimensional tube. From the point of 
view of a control volume attached to (and traveling with) the wave, the flow is from 
left to right. 

Figure CC-0 o.c. 
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This equation 9/1 relates the density change across the wave to the velocity 
change between inlet and outlet. Re-arranging, and focusing only on the case 
where the pressure wave is extremely weak (as are sound waves in practice), 
we see the product dp dV vanish and obtain; 

pc = pc - p dV + c dp - dpdV 
pdV = c dp (9/2) 

In equation 9/2, which relates the speed of sound c to the density p, we would 
now like to eliminate the dV term. For this, we turn to the control volume 
momentum equation (eq: 3/8 p.65): 

fnet ~ “diinVin + dioutkout 

pA - {p + dp)A = -pcAc + {p - dp){c - dV)A{c - dV) 

- dpA = pcA [-C -t (c - dR)] 

— =pdV (9/3) 

c 

We can now combine eqs. (9/2) and (9/3) to obtain: 


c = 



(9/4) 


This can be generalized for any pressure wave, and we state that 



where the s = cst. subscript indicates that entropy s is constant. 


(9/5) 


So now, we see that in a fluid the speed of sound —the speed of pressure 
waves when they travel in a reversible (constant-entropy) manner— is the 
square root of the derivative of pressure with respect to density. This allows 
us to compare it to another term, the bulk modulus of elasticity K, which 
represents the differential amount of pressure that one has to apply on a 
body to see its density increase by a certain percentage: 



(9/6) 


where K is expressed in Pa. 


Inserting eq. 9/6 into eq. 9/5, we obtain a relationship which we had intro¬ 
duced already in chapter 6 (as eq. 6/20 p.l38): 


c = 




Rs=cst. 

p 


(9/7) 


With this equation 9/7 we see that in any fluid, the square of the speed of 
sound increases inversely proportionally to the density and proportionally to 
the modulus of elasticity (the “hardness”) of the fluid. 
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9.2.3 The measure of compressibility 


The most useful measure of compressibility in fluid mechanics is the Mach 
number, which we already defined (with eq. 0/10 p.l5) as: 

[Ma] = - (9/8) 

c 

As we saw in chapter 6 (§6.4.3 p.l38), [Ma]^ is proportional to the ratio of 
net force to elastic force on fluid particles: at high Mach numbers the forces 
related to the compression and expansion of the fluid dominate its dynamics. 

Not all compressible flows feature high Mach numbers. In large-scale atmo¬ 
spheric weather for example, air flows mush slower than the speed of sound 
yet large density changes can occur because of the vertical (static) pressure 
gradient. The same is true in low-speed centrifugal machines because of 
centrifugal pressure gradients. In this chapter however, we choose to focus 
on the simplest of compressible flows: one-directional gas flows with Mach 
numbers approaching or exceeding 1. 


9.3 Thermodynamics of isentropic flow 
for a perfect gas 


9.3.1 Principle 

When a gas is compressed or expanded, its pressure and temperature vary 
together with its density. Before we consider the dynamics of such phenom¬ 
ena, we need a robust model to relate those properties. This is the realm of 
thermodynamics. 

In short, from a macroscopic point of view, a gas behaves as a mechanical 
spring (fig. 9.2): one needs to provide it with energy as work in order to 
compress it, and one recovers energy as work when it expands. Three 
important points must be made regarding these compressions and expansions: 


1. When gases are heated during the movement, the work transfers are in¬ 
creased; conversely they are decreased when the gas is cooled (fig. 9.3). 
Cases where no heat transfer occurs are called adiabatic. 




Figure 9.2 - When they are compressed or expanded, the behavior of fluids can be 
modeled as if they were mechanical springs: they gain energy during compressions, 
and lose energy during expansions. These spring can be thought of as being “fragile”: 
when the movements are sudden, compressions involve more mechanical energy, 
and expansions involve less. 

Figure CC-by-sa Olivier Cleynen 
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p 


infinitesimally slow (reversible) evolutions: 



V 


Figure 9.3 - Pressure-volume diagram of a perfect gas during a compression. When 
the fluid is heated up, more work is required; when it is cooled, less work is required. 
Opposite trends are observed during expansions. 

Figure CC-0 Olivier Cleynen 




Figure 9.4 - When no heat transfer occurs, the a compression or expansion is said 
to be adiabatic. The temperature of fluids can still vary significantly because of the 
compression or expansion. When adiabatic evolutions are infinitely smooth, they 
are termed isentropic. 

Figure CC-0 Olivier Cleynen 
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2. The temperature of gases tends to increase when they are compressed, 
and tends to decrease when they expand, even when no heat transfer 
occurs (fig. 9.4); 

3. A gas behaves as a “fragile” spring: the faster a compression occurs, 
and the more work is required to perform it. Conversely, the faster 
an expansion occurs, and the less work is recovered from it. Cases 
where evolutions occur infinitely slowly are used as a reference for 
comparison: such cases are termed reversible. Cases where the evolu¬ 
tion is both adiabatic and reversible are named isentropic {iso-entropic, 
because they occur at constant entropy). 


9.3.2 Pressure, temperature, density and the speed of sound 

The quantitative relations between pressure, temperature and density of 
perfect gas during isentropic evolutions are classical results in the study of 
thermodynamics. In this section, we wish to re-write them with particular 
focus placed on the speed of sound and the Mach number. The goal is to 
arrive to eqs. 9/20 to 9/22 p.l99, which we need in the upcoming sections to 
describe supersonic flow of gases. 

Let us then focus on the isentropic evolutions of a perfect gas. It is known 
from the study of thermodynamics that in isentropic flow from a point 1 to 
a point 2, the properties of a perfect gas are related one to another by the 
relations: 

(9/9) 

(9/10) 
(9/11) 

for a perfect gas, during an isentropic (constant-entropy) evolution. 



It follows from eq. 9/11 that in an isentropic evolution, a perfect gas behaves 
such that pp~^ = k = cst., which enables us to re-write equation 9/5 for a 
perfect gas as: 



1 = ^Jpp yypy- 



The speed of sound in a perfect gas can therefore be expressed as: 



V 


for any evolution in a perfect gas. 


(9/12) 

(9/13) 


Thus, in a perfect gas, the speed of sound depends only on the local tempera¬ 
ture. 

We now introduce the concept of stagnation or total properties, a measure of 
the total amount of specific energy possessed by a fluid particle at a given 
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instant. They are the properties that the fluid would have if it was brought to 
rest in an isentropic manner: 


ho = h + 

2 

1 1 , 

To = T+ - 

Cp 2 

ho = CpTo 


(9/14) 

(9/15) 

(9/16) 


Using equation 9/13, we can re-write eq. (9/15) as: 


T 


= 1 + 


1 1 


CpT 2 


--V^ = 1 + 


1 1 


I1t2 

r-i 






Eqs. (9/9) to (9/11) can thus be re-written to express the transition from flow 
condition to stagnation condition: 


1 

1 


1 + 


1 + 




(9/17) 

(9/18) 

(9/19) 


These three equations (9/17) to (9/19) express, for any point inside the flow 
of a perfect gas, the local Mach number [Ma] as a function of the ratio of 
the local (current) property (T, p or p) to the local totai property (Tq, po, po)- 
These ratios are always inferior or equal to one, by definition. 

Now, if now we define critical conditions as those that would occur if the 
fluid was exactly at [Ma] = 1, noting them with an asterisk, we can re-write 
eqs. (9/17) to (9/19) comparing stagnation properties to critical properties: 



(9/20) 

(9/21) 

(9/22) 


Here, we compare the local total property (e.g. the total temperature To) to 
the value the property needs to have for the flow to be at sonic speed (e.g. 
the critical temperature T*). This ratio does not depend at all on the fluid 
flow: it is merely a function of the (invariant) characteristics of the gas. 

This is it! It is important to understand that nothing in this section describes 
a particular flow movement. We have merely derived a convenient tool 
expressing flow properties as a function of their stagnation properties and 
the local Mach number, which we shall put to good use further down. 
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9.4 Speed and cross-sectional area 


Let us now go back to fluid dynamics. The last time we considered frictionless 
pipe flow, in chapter 5, we situation was rather simple. From the one¬ 
dimensional mass conservation equation, eq. 5/1: 


piViAi - P 2 V 2 A 2 (9/23) 

it was possible to relate cross-sectional area A and velocity V easily. Now 
that p is also allowed to vary, the situation is decidedly more interesting. 

We wish to let static, pressurized air accelerate as fast as possible: letting V 
increase at a cost of a decease in p according to the rules of an isentropic, 
steady expansion. We know however, that according to mass conservation, 
pVA has to remain constant all along the acceleration, so that only one given 
nozzle geometry {A = f{p,T)) will permit the desired acceleration. How can 
we find this geometry? 

We start from the conservation of mass inside our isentropic acceleration 
nozzle, and differentiating the terms: 


pAV 

Inp -h In A -I- Ink 
dp dA dV 


cst. 

cst. 

0 


for any steady flow. 


Let us set this equation 9/24 aside for a moment. During our isentropic 
expansion, the total energy of the fluid remains constant, thus the total 
specific enthalpy ho is constant. 


h 


y2 



and here also we differentiate, obtaining: 


dh + VdV = 0 


(9/24) 


In last relation, since the process is reversible and without external energy 
input, we can write the first term as dh = dp/p, yielding: 


-dp + VdV = 0 
P 

1 1 

- = -—VdV 
P dp 

for isentropic gas flow. 


(9/25) 


This expression of 1/p is the relation we were after. Now, coming back to 
equation 9/24 and inserting eq. (9/25), we get: 


dp dA dk 


dp 

dk 


dp 


-h 


dA 

~A 


= 0 

= 0 
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In this last relation, since we are in an isentropic flow, we can insert an 
expression for the speed of sound, = 0| derived from equation 9/5, ob¬ 
taining: 


dV 

'Y 

dV 

~ 



dA 
+ — 
A 

dA 
+ — 
A 


0 

0 


This leads us to the following equation relating area change and speed change 
in compressible gas flow: 


dV _ dA 1 

[Ma]' - 1 

for isentropic, one-dimensional gas flow. 


(9/26) 


This equation 9/26 is one of the most stunning equations of fluid dynamics. It 
relates the change in flow cross-section area, dA, to the change in speed dV, 
for lossless flow. Let us examine how these two terms relate one to another, 
as displayed in figure 9.5: 

• When [Ma] < 1, the term l/([Ma]^ - 1) is negative, and dA and dF are 
always of opposite sign: an increase in area amounts to a decrease in 
speed, and a decrease in area amounts to an increase in speed. 

• When [Ma] > 1, dA and dV have the same sign, which —perhaps 
disturbingly— means that an increase in area amounts to an increase in 
speed, and a decrease in area amounts to a decrease in speed. 

Thus, if we want to let a static, compressed fluid accelerate as fast as possible, 
we need a converging nozzle up to [Ma] = 1, and a diverging nozzle onwards. 
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[Ma] < 1 ^ ^ 

subsonic 



Figure 9.5 - Changes in nozzle cross-section area have opposite effects when the 
Mach number [Ma] is lower or greater than 1. 

Note that because temperature changes as the gas compresses and expands, the local 
Mach number varies as the flow changes speed. It is thus possible for the flow to go 
from one regime to another during an expansion or compression. These effects are 
studied further down and described in figure 9.7 p.204. 

Figure CC-0 o.c. 


9.5 Isentropic flow in converging and diverging 
nozzles 


The change of between the subsonic and supersonic flow regimes, by defini¬ 
tion, occurs at Mach one. When the flow is isentropic, equation 9/26 above 
tells us that the [Ma] = 1 case can only occur at a minimum area. 

Of course, there is no guarantee that a pipe flow will be able to become 
supersonic (fig. 9.6). For this, very large pressure ratios between inlet and 
outlet are required. Nevertheless, whichever the pressure ratio, according 
to eq. 9/26, a converging nozzle will only let an expanding fluid accelerate 
to Mach one (critical conditions). Changing the nozzle outlet area itself will 
not permit further acceleration unless an increase in nozzle area is generated 
beyond the critical point. 

This is a surprising behavior. In a converging-only nozzle, increasing the 
inlet pressure increases the outlet velocity only up until [Ma] = 1 is reached 
at the outlet. From there on, further increases in inlet pressure have no effect 
on the outlet pressure, which remains fixed at poutlet = P* ■ The flow then 
reaches a maximum mass flow rate mmax- This condition is called choked 
flow. 

Once the flow is choked, adding a diverging section further down the nozzle 
may enable further acceleration, if the pressure ratio allows for it. Neverthe¬ 
less, the pressure pT = pthroat at the throat remains constant atpx = P*' from 
this point onwards, downstream conditions (in particular the pressure) are 
not propagated upstream anymore: the gas travels faster than the pressure 
waves inside it. Withpx still equal to the critical pressure p*, the mass flow 
through the nozzle cannot exceed mmax- 





















Figure 9.6 - Non-dimensionalized properties of a gas undergoing an isentropic (i.e. 
shock-less) expansion through a converging-diverging nozzle. 

Figure CC-0 o.c. 


We can see, therefore, that it follows from equation 9/26 that there exists 
a maximum mass flow for any duct, which depends only on the stagnation 
properties of the fluid and the throat cross-sectional area. 

The mass flow, at any point in the duct, is quantified as: 
m = pcA [Ma] 

-r-i 

2{y-1) 

(9/27) 

and m will now reach a maximum when the flow is choked, i.e. reaches 
[Ma] = 1 at the throat, where its properties will be p* = pj, T* = Tj, and the 
area A* = Aj: 

,- 

2(r-i) / y 

A general description of the properties of a perfect gas flowing through a 
converging-diverging duct is given in figs. 9.6 and 9.7. 




r + i 


m = A[Ma]po 


RF 


1 -t 


(r-l)[Ma]' 
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Figure 9.7 - The pressure of a gas attempting an isentropic expansion through a 
converging-diverging duct, represented as a function of distance for different values 
of outlet pressure pb* The green curves show the reversible evolutions. 

As long as the back pressure pb is high enough, the flow is subsonic everywhere, 
and it remains reversible everywhere. Whenpr reaches p*, the flow downstream of 
the throat becomes supersonic. In that case, if the back pressure is low enough, a 
fully-reversible expansion will occur. If not, then a shock wave will occur and the 
flow will return to subsonic regime downstream of the shock. 

Figure CC-0 o.c. 


9.6 The perpendicular shock wave 


The movement of an object in a fluid at supersonic speeds (or a supersonic 
fluid flow over a stationary object, which produces the same effects) provokes 
pressure field changes that cannot propagate upstream. Thus, if a direction 
change or a pressure increase is required to flow past an obstacle, fluid 
particles will be subjected to very sudden changes. These localized changes 
are called shock waves. 

The design of bodies moving at supersonic speeds usually focuses on ensuring 
pressure and direction changes that are very gradual, and at the limit (e.g. far 
away from a supersonic aircraft) these are very small. 

The limiting case at the other side of the spectrum is that of the normal shock 
wave, the most loss-inducing (and simple!) type of shock wave. This will 
typically occur in front of a blunt supersonic body. 

A perpendicular shock wave (fig. 9.8) is studied with the following hypothesis: 

• It is a constant energy process; 

• It it a one-dimensional process (i.e. the streamlines are all parallel to 
one another and perpendicular to the shock wave). 

















With these two hypothesis, we can compare the flow conditions upstream 
and downstream of a shock wave occurring in a tube using the continuity 
equation, which we wrote above as eq. 9/23: 


piAiVi — P 2 A 2 V 2 (9/29) 

The energy equation in this case shrinks down to: 

Toi = To 2 (9/30) 

Since, from equation 9/17 p.l99, we have: 



and given that Tqi = T 02 , dividing one by the other yields: 

T 2 _ 1 + ^[MaJi 

Ti i + I^[Ma]2 

In a similar fashion, we can obtain: 

^ _ 1 + )^[Ma]i 

pi 1 + )'[Ma]2 


(9/31) 


(9/32) 


Let us set aside these two equations 9/31 and 9/32 and compare conditions 
upstream and downstream starting with mass conservation, and re-writing it 
as a function of the Mach number: 


Plkl - ^2^2 

Pi _ 1^2 _ \Ma\ 2 Cl 2 

P2 Vi [MaJifli 

[Ma\2^/YR^2 


[Mali ^JYRTl 

[Mall ^P2 /ri\^ 
[Ma]2 pi \T2 / 


(9/33) 


So far, we have not obtained anything spectacular, but now, inserting both 
equations 9/31 and 9/32 into this last equation 9/33, we obtain the staggering 


1 ' 

^ 1 

Pi 

T, 

- 1 - 

[Ma], i 

1 

: [Ma]^ P2 > p, 

T,>T, 

Pi 

Toi 

1 

1 

1 

\ Pl^ P\ 

1 

A = r., 



cross-section area A " " 7 " 1 1 j u 1 

control volume around shock wave 


Figure 9.8 - A normal shock wave occurring in a constant-section duct. All properties 
change except for the total temperature To, which, since the transformation is at 
constant energy, remain constant. 


Figure CC-0 o.c. 
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result: 


[Ma]i _ 1 + r[Ma]i /1 + 

[Ma]2 1 + y[Ma]l \ 1 + ^[MaJj / 


(9/34) 


This equation is not interesting because of its mathematical properties, but 
because of its contents: the ratio [Ma]i/[Ma ]2 depends only on [Ma]i and 
[Ma] 2 ! In other words, if the shock wave is perpendicular, there is only 
one possible post-shock Mach number for any given upstream Mach 
number. 

After some algebra, [Ma ]2 can be expressed as 


[Ma]2 = 


[Ma]? + ^ 

2[Ma]^^-l 


(9/35) 


Values for [Ma ]2 as a function of [Ma]i according to equation 9/35 for par¬ 
ticular values of y can be tabulated or recorded graphically, allowing us to 
predict properties of the flow across a shock wave. We can therefore quantify 
the pressure and temperature increases, and velocity decrease, that occur 
though a perpendicular shock wave. 


9.7 Compressible flow beyond frictionless pipe 
flow 


In this chapter, we have restricted ourselves to one-directional steady flow, 
i.e. flows inside a pipe; within these flows, we have neglected all losses but 
those which occur in perpendicular shock waves. This is enough to describe 
and understand the most important phenomena in compressible flow. 

Of course, in reality most compressible flows are much more complex than 
that. Like other flows, compressible flows have boundary layers, feature 1am- 
inarity and turbulence, heat transfer, unsteadiness, and more. Shock waves 
occur not just perpendicularly to flow direction, but also in oblique patterns, 
reflecting against walls and interacting with boundary layers. Temperature 
variations in compressible flow are often high enough to provoke phase 
changes, which complicates the calculation of thermodynamic properties as 
well as flow dynamics and may result in spectacular visual patterns. 

All those phenomena are beyond the scope of this study, but the reader can 
hopefully climb up to more advanced topics such as high-speed aerodynamics, 
or rocket engine and turbomachine design, where they are prominently 
featured and can surely make for some exciting learning opportunities. 

Video: watch a 140-thousand- 
ton rocket lift away humans on 
their first ever trip to the moon. 

Compressible flow of expand¬ 
ing exhaust gases is what al¬ 
lows the maximum amount of 
thrust to be obtained from a 
given amount of combustion en¬ 
ergy in rocket engines. 

by NASA (public domain) 
https://youtu.be/eZDaM3MHNqg 
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Exercise sheet 9 (Compressible flow) 

last edited April 2, 2018 


These lecture notes are based on textbooks by White [13], Qengel & al.[16], and Munson & al.[18]. 


Except otherwise indicated, we assume that fluids are Newtonian, and that: 

Pwater = 1 000 kg m“^ patm. = 1 bar; Patm. = 1,225 kgm“^ Tatm. = ll,3°C;patm. = 1,5 • 10‘^Nsm“2; 
g = 9,81 ms“^. Air is modeled as aperfect gas (i?air = 287JK“^ kg“'; pair = 1,4; Cpair = 1005Jkg“^ K“^). 


The speed of sound c in a perfect gas is a local property expressed as: 


c = 


(9/12) 


In isentropic, one-dimensional flow, we accept that the mass flow m is quantified as: 


m = A[Ma]po 


Y 

RTo 


1 + 


(y - l)[Ma]' 


-y-l 

2(r-i) 


(9/27) 


This expression reaches a maximum /hmax when the flow is choked: 




r + ij 


2(y-l) 


A*po 


Y 

RTo 


(9/28) 


The properties of air flowing through a converging-diverging nozzle are described in 
fig. 9.9, and fig. 9.10 describes air property changes through a perpendicular shock 

wave. 
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One-dimensional isentropic compressible flow functions for an ideal 
gas with /c = 1.4 
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Figure 9.9 - Properties of air (modeled as a perfect gas) as it expands through a converging- 
diverging nozzle. 

These are numerical values for equations (9/17) to (9/19). In this figure, the perfect gas parameter 
Y is noted k. Data also includes the parameter [Ma]* = y/c* (speed non-dimensionalized relative 
2 Qg to the speed of sound at the throat). 

Figure © Qengel & Cimbala 2010 [16] 

























One-dimensional normal shock functions for an \de^\ gas with k = lA 
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Figure 9.10 - Properties of air (modeled as a perfect gas) as it passes through a perpendicular 
shock wave. 

These are numerical values calculated based on equation 9/35. In this figure the perfect gas 
parameter y is noted k. 

Figure © Qengel & Cimbala 2010 [16] 
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9.1 Flow in a converging nozzle , , 

® ® Qengel & al. [ 16 ] E12-5 

Air from a 4 mm hole in a car tire leaks into the atmosphere. The pressure and tempera¬ 
ture in the tire are 314 kPa and 25 °C. The exterior air conditions are 94kPa and 16 °C. 

If the flow can be approximated as isentropic, what is the air flow rate? 


9.2 Flow in a converging-diverging nozzle . „ 

® ~ O O £xam WS2015-16, see also White [ 13 ] E9.4 

An experiment is conducted where compressed air from a large reservoir is expanded 
and accelerated through a nozzle (figure 9.11). There is no external work transfer to 
or from the air, and the expansion can satisfactorily be modeled as a one-dimensional 
isentropic process (isentropic meaning that the process is fully reversible, i.e. without 
losses to friction, and adiabatic, i.e. without heat transfer). 



Figure 9.11 - A nozzle in which compressed air in a reservoir (left) is expanded as it accelerates 
towards the right. The back pressure pb is controlled by a valve (represented on the right). The 
flow conditions of interest in this exercise are labeled 1, 2 and 3. 

The air (cpair = 1 005Jkg“^ K“9 flows through a converging-diverging nozzle. The en¬ 
trance conditions into the nozzle are Ai = 0,15 m^, pi = 620 kPa, Vi = 180 m s“\ and Ti = 
430 K. 

In the expanding section of the nozzle, we wish to know the conditions at a point 2 of 
cross-sectional area A 2 = Ai = 0,15 m^. 

1 . What is the mass flow through the nozzle? 

2 . If we wish the flow to become sonic (i.e. reach [Ma] = 1 ), what is approximately 
the required throat area Aj? 

3. Which two flow regimes can be generated at section 2 ? What is approximately the 
pressure p 2 required for each? 

The nozzle is now operated so that the flow becomes sonic at the throat, and accelerates 
to supersonic speeds in the expanding section of the nozzle. 

The nozzle cross-sectional area increases further downstream of point 2 , up until it 
reaches an area A 3 = 0,33 m^. Downstream of point 3, the duct is of constant area and 
210 the back pressure pb is set to pb = ps = p 2 - 














4. Describe qualitatively (i.e. without numerical data) the pressure distribution 
throughout the nozzle. 

Finally, we wish to expand the air with a fully-isentropic process, expanding the air from 
the reservoir to point 3 without a shock wave. 

5. Propose and quantify one modification to the experimental set-up that would allow 
an isentropic, supersonic expansion of the flow into the outlet of area A 3 . 


9.3 Capsule re-entry ccooc 

A Soyuz capsule incoming from the international space station is passing through 
Mach 3,5 at a pressure of 0,2 bar and temperature of -50 °C. A bow shock wave is 
generated in front of the bottom surface of the capsule, which is approximately normal 
to the flow for the particles directly in the path of the capsule. 

What are the conditions (temperature, velocity, density and pressure) in between the 
shock wave and the capsule shield? 


9.4 Explosion shock wave 

^ CC-O o.c. 

In an unfortunate coincidence, on the morning following the fluid dynamics exam, the 
red Tesla S of the lecturer explodes on a parking lot (fortunately, no one was hurt). The 
explosion creates an spherical shock wave propagating radially into still air. 

At a given moment, the pressure inside the spherical bubble is 1 380 kPa. The ambient 
atmospheric conditions are 17 °C. 

At what speed does the shock wave propagate, and what is the air speed immediately 
behind the wave? 


9.5 Home-made supersonic tunnel 

^ CC-O o.c. 

A PhD student is thinking about constructing his/her own supersonic wind tunnel. 
To this purpose, s/he designs a converging-diverging tube with variable square cross- 
section. The throat T is 0,1 X 0,1 m; the tunnel expands progressively to an outlet nozzle 
of 0,13 X 0,13 m cross-section. The outlet exits directly in the ambient atmosphere. On 
the inlet side of the tunnel, the student hopes to install a very large concrete reservoir 
(fig. 9.12). 

1 . What is the maximum reservoir pressure for which the flow in the tunnel is 
subsonic everywhere? 

2 . What reservoir pressure would be required to generate a shock-less expansion to 
supersonic flow in the tunnel? 

3. If the outlet temperature was to be ambient (11,3 °C), what would the outlet velocity 
be? 

4. Which mass flow rate would then be passing in the tunnel? 

5. Describe qualitatively (i.e. without numerical data) the pressure distribution 
throughout the nozzle for both flow cases above. 
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Figure 9.12 - A nozzle in which compressed air in a reservoir (left) is expanded as it accelerates 
towards the right. The back (outlet) pressure pb is always atmospheric. 


In order to impress his/her friends, the PhD student wishes to generate a perpendicular 
shock wave inside the tunnel at exactly one desired position (labeled 1 in fig. 9.12 above). 
At this position, the cross-section area is 0,01439 m^. 

6 . Describe qualitatively (i.e. without numerical data) the pressure distribution 
throughout the nozzle for this case, along the two other pressure distributions. 

7. [difficult question] What reservoir pressure would be required to generate the shock 
wave at exactly the desired position? 

8 . What would the outlet velocity approximately be? 
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Answers 


9.1 Since P 2 /P 0 < 0,5283, the flow is choked. Since there is no nozzle expansion, 
the flow cannot accelerate beyond [Ma] = 1 and it will exit the hole with excess 
pressure. With equation 9/28 we obtain m = 9,232 • 10“^ kg s“^ (0,55 kg min“^). 

9.2 1) m = mi = piViAi = 135,6kgs“^; 

2 ) To position point 1 in fig. 9.9 we can either calculate [Ma]i = Vi/ui = 0,433 or 

Ti/To = Ti/ |Ti + j = 0,9639. Either way we readAi/A* « 1,5901, which gets 

us the throat area required for sonic flow: A* = 9,433 • 10“^ m^; 

3) We know the flow is isentropic, so at 2 we can have either p 2 = pi (flow subsonic 
everywhere), or p 2 ~ 1,0446 • 10^ Pa (flow becomes supersonic, p 2 obtained from 
fig. 9.9 after inputing ZI 2 /A*); 

4) This would be the third case in fig. 9.7 p. 204; 

5) The easiest answer is to propose lowering p 3 to 0,0368po- But we could also 
decrease Aj so that A* = A^/3,5 or increase po to l/0,0368p3. 

9.3 All of the answers are listed in fig. 9.10. In 1-D flow, there is only one possible 
downstream solution. Interpolating between [Ma]i = 3 and [Ma]i = 4, we get 
[Ma ]2 = 0,4551, p 2 = 2,883 bar (!!), p 2 = 1,316 kg m“^ T 2 = 750,4 K (477 °C!) and V 2 
a mere 249,9 ms“^ (900kmh“^), down from Vi = 1048ms“^. Re-entering from 
space flight is tough! 

9.5 1) Flow subsonic everywhere: po = 1,103 bar; 

2 ) Flow supersonic in diverging section: po = 7,8247 bar; 

3) We would have [MaJb = 2 and Vb = 676,1 m s“^; 

4) m = 13,98 kg s“^ 

7) Before shock, [MaJi = 1,8; after it, [Ma ]2 = 0,6165. Now, expand the post-shock 
flow subsonically to outlet condition to findpb/po 2 ~ 0,843. From po 2 , go back 
upstream through the shock: poi = ptank = 1,4596 bar; 

8 ) From Tq, find Ti (isentropic expansion), then T 2 (through shock), then Tb (sub¬ 
sonic isentropic expansion): Tb « 487,4 K. Now, Vb = [MaJb -^yRTb « 221,3 ms“k 
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Appendix 



A1 Notation 


= By definition. The = symbol sets the definition of the term on its left 
(which does not depend on previous equations). 

(dot above symbol) Time rate: ' = For example, Q is the rate of heat 
(in watts) representing a heat quantity Q (in joules) every second. 

(bar above symbol) Time-average: A = avg(A) = avg(A-l-A'). The prime 
symbol indicates the instantaneous fluctuation around the average. 

A Indicates a net difference between two values: (AX)a^b = ~ ^a- 

Can be negative. 

italics Physical properties (e.g. mass m, temperature T). 

straight subscripts Points in space or in time (temperature Ta at point A). 
Subscripts “cst” indicate a constant property, “in” indicates “incoming” 
and “out” is “outgoing”. 

Subscript “av.” indicates “average”. 

lowercase symbols Specific values: property per unit mass. For example, 
b = B/m. 

operators Differential d, partial differential d, finite differential 5, total (alt.: 
subtantial) derivative D/Dt (def. eq. 4/7 p.88), exponential expx = e^, 
natural logarithm In x = log^ x. 

vectors Vectors are always written with an arrow. Velocity is V = {u,v,'w), 

alternatively written u; = {u,v,'w). The norm of a vector A (positive or 

—^ ^ 

negative) is |A|, its length (always positive) is ||A||. 

vector calculus Dot product A • B; cross product A A B; gradient VA (def 
eq. 1/8 p.32); divergent V • A (def eq. 2/10 p.50); Laplacian V^A (def 
eq. 4/30 p.94); curl V X A (def eq. A/16 p.218). 

units Units are typed in roman (normal) font and colored gray (1 kg). In 
sentences units are fully-spelled and conjugated (one hundred watts). 
The liter is noted L to increase readability (1L = 10“^ m^). Units in 
equations are those from systeme international (si) unless otherwise 
indicated. 

numbers The decimal separator is a comma, the decimal exponent is pre¬ 
ceded by a dot, integers are written in groups of three (1,234 -10^ = 1 234). 
Numbers are rounded up as late as possible and never in series. Leading 
and trailing zeros are never indicated. 
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A2 Field operators 


Four operators which apply on vector or scalar fields are important in fluid 
mechanics: gradient, divergent, Laplacian and curl. 


A2.1 Gradient 


The mathematical operator gradient (first introduced as eq. 1/8 p. 32) is 
written V . It applies on a scalar field and produces a vector field. It is 
defined as: 


^ ^ d 

dx dy dz 

dA^ dA^ dAt 
V.A = “—i + ~—j + “—k — 
dx dy dz 


¥ \ 


ax 

dA 


dy 
dA 
\ dz 


(A/1) 

(A/2) 


For example, the gradient of a pressure field is the vector field Vp: 



dp ^ dp ^ 

A-j + + 

dx dy 



/ ^ \ 

dp 

F 

dp 
Tz / 


(A/3) 


A2.2 Divergent 


The mathematical operator divergent (first introduced as eq. 2/10 p.50) is 
written V- and is defined as: 


g _ d ^ d ^ d -> 

Tz’"' 

When applied on a vector field, it produces a scalar field: 

d ^ ^ d ^ d ^ ^ 

V • A = — j • A + -x-j-A + -x-k-A 
dx dy dz 


dAy 


dA 


+ 


y dA^ 

+ 


dx dy dz 

When applied on a 2"'’ order tensor field, it produces a vector field: 

/ dAy 


V • Aij = 


dx 

dA- 


XX ^ dA, 


dx 


dy 

l dA 


\ dx 
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dy 

dA, 


dz 


zy 


yz 


dy 


+ 


dz 
dA,, 
dz / 
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( V-A-x \ 

= 

V • Aly 

i 

\ ^ * Aiz j 


(A/4) 


(A/5) 

(A/6) 


(A/7) 


—T 

For example, the divergent of a velocity field is the scalar field X -V: 


V • V 


du dv dw 
dx dy dz 


(A/8) 
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A2.3 Laplacian 


The mathematical operator Laplacian (first introduced as eq. 4/30 p.94) is 
written and defined as: 


= V • V 


(A/9) 


When applied to a scalar field, it is equal to the divergent of the gradient of 
the field, and produces a scalar field: 

V^A = V • VA 

d^A d^A d^A 

{dxy {dyY {dzY 


(A/10) 

(A/ll) 


When applied to a vector field, it is equal to the gradient of the divergent of 
the field, and produces a vector field: 


VM = 


V^A;, \ 
V^A,, 


{dxY 

d'^Ay 

{dxY 
\ (dx)^ 


\ 

V • VA;, \ 

= 

V-VA, 

/ 

1 V • VA, / 

d^Ax 

d'^Ay 

{dyf 

d^A^ 

(dy)^ 

, d^Ax \ 
(5z)2 
d^Ay 

{dzr 
, d^A, 

^ (5z)2 / 


(A/12) 


(A/13) 


For example, the Laplacian of a velocity field is the vector field V W: 



(A/14) 


A2.4 Curl 


The mathematical operator curl (sometimes named rotational) is written Vx . 
It applies to a vector field and produces a vector field. It is defined as: 


Vx = 


V X A = 


i j k 
AAA 

dx dy dz 


I 

_d_ 

dx 

Av 


J 

dy 

A. 
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k 

dz 

A. 
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dAz dAx 
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dAy dAy 


dx 


dy 
(A/16) 


For example, the curl of velocity is the vector field V xV: 


VxV = 


i 

A. 

dx 


j 

_d_ 

dy 


k 

dz 


U V W 


(dw dv\^ ( dw du\^ Idv duX-^ 
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